Spectrometric systems and mechanical methods for improved focus localization of time and space-varying measurements

ABSTRACT

A system and method of dynamically localizing a measurement of parameter characterizing tissue sample with waves produced by spectrometric system at multiple wavelengths and detected at a fixed location of the detector of the system. The parameter is calculated based on impulse response of the sample, reference data representing characteristics of material components of the sample, and path lengths through the sample corresponding to different wavelengths. Dynamic localization is effectuated by considering different portions of a curve representing the determined parameter, and provides for the formation of a spatial map of distribution of the parameter across the sample. Additional measurement of impulse response at multiple detectors facilitates determination of change of the measured parameter across the sample as a function of time.

CROSS-REFERENCE TO RELATED APPLICATIONS

This patent application is a continuation from the U.S. patentapplication Ser. No. 14/397,814 filed on Oct. 29, 2014 and published asU.S. 2015/0109617, which is a national stage entry of InternationalPatent Application No. PCT/US2013/038623 filed on Apr. 29, 2013, which,in turn, claims the benefit of and priority from U.S. provisional patentapplication Ser. Nos. 61/640,682, filed on Apr. 30 2012; 61/789,565,filed on Mar. 15, 2013; 61/794,196, filed on Mar. 15, 2013; and61/794,515, filed on Mar. 15, 2013. The disclosure of each of theabovementioned patent applications is incorporated herein by referencein its entirety.

TECHNICAL FIELD

The present invention relates to spectroscopic systems and methods oftransmitting and receiving two or more concurrent signals and systemsand methods enabling improved measurement focus localization of time-and space-varying spectrometer measurements, with additional uses forvolume and volume variation measurements and challenge based propertymeasurements.

BACKGROUND OF THE INVENTION

The assessment of arterial blood and tissue oxygen saturation has beenshown to be critical for monitoring, diagnosing, and treating acutecardiovascular deficiencies, for example, but not limited to, low bloodoxygen saturation due to chronic obstructive pulmonary disease (COPD),or in extreme cases, exsanguinations. Further, the time line over whichphysiological changes occur is indicative of underlying cardiovasculardeterioration. While techniques exist that allow qualitative monitoringof the level of arterial blood and tissue oxygen saturation, the currentmedical practice appears to lack a tool for quantitative monitoring ofdeteriorating cardiovascular status from the initial cardiovascularevent or acute trauma through arrival at a medical facility.

Spectroscopy was originally the study of the interaction betweenradiation and matter as a function of wavelength (“λ”). Historically,spectroscopy referred to the use of visible light dispersed according toits wavelength, e.g. by a prism.

Later the concept of spectroscopy was expanded greatly to comprise anymeasurement of a quantity as a function of either wavelength orfrequency. Thus, it also can refer to a response to an alternating fieldor varying frequency (“ν”). A further extension of the scope of thedefinition added energy (“E”) as a variable, once the very closerelationship “E”=“hν” for photons was realized (“h” is Planck'sconstant). A plot or measurement of the response of a material orstructure as a function of wavelength—or more commonly frequency—isreferred to as a spectrum.

ISO Standard number 80601-2-61:2011 states, on page 34: “Currenttechnology requires an adequate concentration of haemoglobin, apulsatile change in blood flow, and light transmission through a tissuebed to approximate the in vivo haemoglobin oxygen saturation. PULSEOXIMETER EQUIPMENT is not typically capable of functioning effectivelyduring cardiopulmonary bypass or at extreme low-flow states, and is notat present intended as a means for the measurement of blood flow orblood volume.

“Given these limitations, PULSE OXIMETER EQUIPMENT does not provideprecise measurements of arterial haemoglobin saturation. The presentlymarketed in vivo PULSE OXIMETER EQUIPMENT is not a replacement formeasurement of blood samples by in vitro optical oximeters. The valuesderived from pulse oximetry are not a measurement of blood orsolid-tissue oxygen tension. Pulse oximetry provides no directindication of oxygen delivery to tissue, or of tissue oxygenconsumption.”

Recent theoretical developments have pointed to the potential advantageof a spectroscopic device that uses multiple wavelengths of light toperform oximetry. Spectroscopic devices (including oximeters) usingmultiple well-chosen targeted wavelengths for operation have improvedaccuracy without the use of complex and/or large peripheral systems.There remains a need for a method enabling the selection of multipleoptimal wavelengths for operation of a spectroscopic device configuredto acquire and process data for use in quantitative cardiovascularmeasurements. The targeted selection of multiple wavelengths are alsobeneficial for optimization of the process of accurately and completelydistinguishing between or among different analytes, referred to hereinas analyte contrast.

Furthermore, there remains a need for improved measurement devicesconfigured to operate to measure relative and/or absolute concentrationsof material components in a sample. Conventional models appear to haveignored to incorporate the path length—and especially wavelength-basedvariations in the path length—of an electromagnetic wave traversing themeasured sample. Incorporating this information can increase thesensitivity and accuracy of the measurements, providing a wealth ofevaluation of the property of the material components of interest.Non-limiting examples of these material components are cells, proteins,hemoglobin, glucose, lipids, chromophores, water, pH, and gases such ashyperpolarized gases, carbon dioxide, carbon monoxide, and oxygen. Dueto short transit times of signals through a tissue sample, and the richinformation conveyed in the variations in the impulse and therefore theimpulse response of the sample, the configuration of accurate andfine-grained (fine fidelity) measurement of the impulse response by aspectrometric system is not trivial but, if realized, may producevaluable results.

Given the capability to utilize the fine-grained information for moreaccurate measurements and assessments, there is a need for systems,methods, and apparatus that incorporate the ability to accuratelymeasure very small signal levels and very short transit time therebyproviding the details of the collective differential path lengths of thetransmitted signals.

SUMMARY OF THE INVENTION

Embodiments of the invention include a method for operating aspectrometric device configured to generate electromagnetic waves (EMWs)and juxtaposed with a sample such as to measure at least one property ofat least one material component of the sample using at least one of theEMWs. The method comprises: (i) in response to a first excitationsequence of pulses generated by a unit of the spectrometric device,transmitting a first EMW at a first wavelength by a transmitterassociated with the spectrometric device; (ii) concurrently withtransmitting the first EMW and in response to a second excitationsequence of pulses generated by the unit of the spectrometric device,transmitting a second EMW at a second wavelength by a transmitterassociated with the spectrometric device; and (iii) acquiring said firstand second transmitted EMWs as an aggregated input with a receiverassociated with the spectrometric device. The method additionallyincludes the step of processing the received aggregated input with adata-processing unit associated with the spectrometric device by: (a)correlating the aggregated input with a time-controlled representationof the first excitation sequence of pulses to form a first identifier;(b) correlating the aggregated input with a time-controlledrepresentation of the second excitation sequence of pulses to form asecond identifier; and (c) calculating, based at least on first andsecond identifiers, the at least one property of the at least onematerial components of the sample.

Embodiments of the invention further provide an apparatus adapted tocharacterize first and second material components of a sample, whichapparatus comprises two or more transmitters, each configured toconcurrently transmit a respectively-corresponding electromagnetic wave(EMW) at a respectively-corresponding wavelength in response to anexcitation sequence received by the two or more transmitters; a receiveradapted to receive the EMWs concurrently transmitted EMWs as anaggregated input; and a data-processing unit containing at least onecorrelator and structured to process the aggregated input. In suchapparatus, a distinct correlator is used to quantify each EMW-componentof the aggregated input such as to correlate the aggregated input iswith a time-controlled representation of the excitation sequence togenerate a data sequence and to uniquely quantify the generated datasequence for each EMW-component of the aggregated signal.

Embodiments additionally provide a system for measuring concentrationsof first and second material components in a sample. The system includesa transmitter enabled to transmit a sequence of pulses in response to anexcitation sequence received by the system, a receiver enabled toreceive at least some of the transmitted sequence of pulses; and adevice configured to quantify data representing the received pulses. Thedevice contains a computational unit configured to (a) correlate thereceived pulses with the sequence of pulses; (b) aggregate multiplereceived pulses to quantify a magnitude of data representing pulsesreceived over time; and (c) determine at least one value representing atleast one of the first and second material components from thequantified magnitude.

Embodiments of the invention further provide a method for monitoringcharacteristics of a biological sample with a spectrometric system thatincludes at least one transmitter and at least one receiver operablyjuxtaposed with the sample. The method comprises receiving, with adata-processing unit of the system, an output from the at least onereceiver caused by emission (by the transmitter), of at least oneelectromagnetic wave at at least one corresponding wavelength, where theoutput represents an optical property of the tissue. The method furtherincludes receiving, with the data-processing unit of the system,reference data representing empirically-defined spectrally-dependentcharacteristic of the tissue. The method additionally includesdetermining an impulse response of the sample, the impulse responsebeing associated with emission of the at least one electromagnetic waveby the at least one transmitter; and determining, as a function of time,a wavelength-dependent material parameter characterizing the samplebased on the determined impulse response and the reference data. Thedetermining, as a function of time, a wavelength-dependent materialparameter characterizing the sample may include determining awavelength-dependent material parameter characterizing the sample basedon the averaged wavelength-dependent path length that has been estimatedwith the use of the determined impulse response.

In a related embodiment, the method additionally includes sampling acurve representing the determined wavelength-dependent materialparameter as a function of time to obtain values of the materialparameter at sampled points in time, and forming a weighting function bycalculating variance values among the values of said material parameterat sample points in time, the weighting function being defined by saidvariance values as a function of time. Furthermore, the method mayinclude determining, as a function of time, a modifiedwavelength-dependent material parameter characterizing the sample basedon (i) the computed impulse response that has been weighted by theweighting function and (ii) the reference data.

Embodiments of the invention additionally provide an optical-pathadjustment apparatus for modifying an optical path between a source oflight and a receiver of light. Such apparatus includes: (i) the sourceof light with at least one transmitter enabled to emit a correspondingbeam light at at least one corresponding wavelength; (ii) the receiverof light in optical communication with the at least one transmitter; and(iii) a device including (iii-a) a first plurality of selectivelyreorientable reflectors disposed such that, when each of the reflectorsfrom the first plurality is in a corresponding reference orientation,the first plurality defines a first surface; (iii-b) a second pluralityof selectively reorientable reflectors disposed such that, when each ofthe reflectors from the second plurality is in a corresponding referenceorientation, the second plurality defines a second surface that istransverse to the first surface and substantially parallel to thecorresponding beam of light. The first and second pluralities are inoptical communication with the at least one transmitter such as toenable the apparatus to scan the corresponding beam of light which, whena reflector of the first plurality and a reflector of the secondplurality are reoriented with respect to the corresponding referenceorientations, has interacted with the reflector of the first pluralityafter having interacted with the reflector of the second plurality.

Embodiments of the invention also provide an optical-data-collectionapparatus. The apparatus includes an optical detector and a lightirradiation system enabled to deliver a beam of light from a startingpoint to the optical detector such as to scan an object that has beenpositioned at an intermediate point along a direction of propagation ofthe beam between the starting point and the optical detector. The lightirradiation system contains at least one light transmitter defining thestarting point, and a plurality of spatially reorientable reflectorsdisposed such that, when at least three of said reflectors are inrespectively corresponding reference states of orientation, the at leastthree of said reflector define a plane that is substantially parallel tothe beam of light. In a specific embodiment, the light irradiationsystem enables controlling a distance traversed by the delivered beamthat has scanned the object by reflecting the beam with a reflector fromthe plurality to redirect the beam towards the optical detector.

BRIEF DESCRIPTION OF THE DRAWINGS

FIG. 1A is a diagram schematically illustrating an embodiment of avertical cavity surface emitting laser (VCSEL) based oximeter.

FIG. 1B is a diagram schematically illustrating an alternate embodimentof a VCSEL-based oximeter.

FIG. 1C is an improved measurement of a cardiac stroke volume that canbe enabled based on an embodiment of the invention.

FIG. 2 is a graph illustrating the spectral energy distribution of threelight emitting diodes (LEDs).

FIG. 3A is a graph illustrating a choice of eight operationalwavelengths in a pulse oximeter of the related art employing LEDs.

FIG. 3B is a graph illustrating a smaller number (four) of operationalwavelengths for oximeters of FIG. 1A and FIG. 1B that still provide foradequate operation. The reduction of the number of the operationalwavelengths is achieved with the use of an algorithm of the presentinvention.

FIG. 3C is a graph illustrating an optimal choice of spectrallyindependent bandwidths of a light-based spectroscopic device, determinedwith an algorithm of the invention. In addition, the graph shows the“penalty” curve (the “cost” curve, or the “figure of merit” curve)representing operational cost of employing such spectroscopic device

FIG. 3D is a graph illustrating the spectral distribution of theextinction coefficient characterizing Bilirubin (a possible confoundinganalyte) in chloroform.

FIG. 4 is a graph illustrating an optimal choice of spectrallyindependent bandwidths of a light-based spectrographic measurementdevice, measuring glucose levels, determined with an algorithm of theinvention

FIG. 5 is a flow chart illustrating an optimization algorithm defining asimulated annealing solution.

FIG. 6 is a flow chart illustrating a sequence of validation of asimulated annealing solution to the optimization algorithm incorporatingthe sequence of steps illustrated in FIG. 5.

FIG. 7 is a diagram illustrating an arm with two detailedcross-sectional views.

FIG. 8A is a graph illustrating linear and nonlinear transfer functions.

FIG. 8B is a set of graphs illustrating the ability to receive andrecover the excitation sequence from additive noise.

FIG. 8C is a set of graphs illustrating the ability of a system of theinvention to receive and recover the impulse response from additivenoise.

FIG. 8D is a diagram relating the impulse response graph from FIG. 8C toradiation paths through an example cross-section of a measurementsample.

FIG. 8E is a modification of the diagram of FIG. 8D relating changes inthe impulse response graph from FIG. 8C to the changing diameter of aheterogeneous object in the example cross section.

FIG. 9A is a diagram showing components of a system adapted to measureone or more properties of one or more material components.

FIG. 9B is a diagram illustrating components of a system adapted toincrease signal-to-noise ratio associated with measurement data obtainedusing a transmitter and receiver.

FIG. 9C is a diagram illustrating components of a system adapted totransmit and receive multiple excitation signal sequences to increasesignal-to-noise ratio and/or decrease motion artifacts of associatedmeasurements.

FIG. 10A is a diagram illustrating an embodiment of apparatus configuredto precisely control the distance from the light-scattering materialcomponent to a receiver.

FIG. 10B is a diagram illustrating an addition to the embodiment shownin FIG. 10A that enables the precise control of the distance from thelight-scattering material component to a receiver in two dimensions.

FIG. 10C is a diagram showing the basic layout of the system of FIG. 10Aand the related geometrical parameters.

FIG. 10D shows the embodiment of FIG. 10A in a different operationalconfiguration.

FIG. 10E shows the embodiment of FIG. 10A in another operationalconfiguration.

FIG. 10F shows the embodiment of FIG. 10A with an auxiliary barrierintroduced in a path of light.

FIG. 10G shows the embodiment of FIG. 10A in a top plan view.

FIG. 10H shows a top plan view of an embodiment related to that of FIG.10A, in which the micro-mirror array is configured with a withrow-to-row offset.

FIG. 10I shows a top plan view of another embodiment of the inventionhaving a micro-mirror array with hexagonal mirrors arranged in ahoneycomb fashion.

FIG. 10J shows a top plan view of another embodiment of the inventionhaving a micro-mirror array with rectangular mirrors arranged in adiamond pattern.

FIG. 10K shows a side plan view of another embodiment of the inventionas specifically applied to disk drives.

FIG. 11 is a diagram illustrating components of a spectrometric systemutilizing a variable delay to recover the impulse response of a sampleto an interrogating input.

FIG. 12 is a diagram illustrating an analog and digital input mixercircuit suitable for use in the system shown in FIG. 11.

FIG. 13 is a graph illustrating the form of an impulse response of anoptically measured sample using an embodiment of the invention.

FIG. 14 is a graph illustrating the variation of measurement output oftwo wavelength sources plotted with respect to the delay value of FIG.11.

FIG. 15A is a diagram illustrating an embodiment of a linear feedbackshift register (LFSR) for use as a source of an excitation sequence.

FIG. 15B is diagram illustrating electronic circuitry associated withthe LFSR embodiment of FIG. 15A.

FIG. 16A is a diagram showing an alternative implementation of anexcitation sequence generator.

FIG. 16B is a diagram showing an alternative implementation of anexcitation sequence generator.

FIG. 16C is a diagram showing another alternative implementation of anexcitation sequence generator.

FIG. 17A is a diagram showing a portion of a single bit digitalcorrelator for use in a receiver of an embodiment of the invention.

FIG. 17B is a diagram of electronic circuitry associated with thedigital correlator of FIG. 17A.

FIG. 17C is a flow chart representing an algorithm for determination ofthe gain setting of the input to the digital correlator of an embodimentof the invention.

FIG. 17D is a diagram representing an N-bit digital correlator for usewith a receiver of an embodiment of the invention.

FIG. 17E is a diagram of electronic circuitry associated with the N-bitdigital correlator of FIG. 17D.

FIG. 18A is a diagram of a system employing an analog delay line as partof the correlator.

FIG. 18B is a diagram representing electronic circuitry of the analogdelay line of FIG. 18A.

FIG. 19 is a diagram depicting an embodiment of a non-programmableprocessing element for producing numeric measured property values.

FIG. 20A is a diagram schematically illustrating an example resettableanalog integrating amplifier.

FIG. 20B is a diagram schematically illustrating a traveling-wave-likeparallel implementation to produce key values of the impulse response inoutput registers.

FIG. 20C is a timing diagram illustrating the control signal sequenceconstraints for the system operation of 20B.

FIG. 21 is a diagram of a system for determining the concentration ofmaterial components based on the received signal stream and timeaccurate representations of the various excitation sequences.

FIG. 22 is a flowchart showing the process for using a challenge tomeasure a property of a material component.

FIG. 23 is a diagram illustrating a common high quality technique forintroducing a measured amount of carbon monoxide to a human subject.

FIG. 24A is a diagram showing the elements used to calculate the impulseresponse of material properties of a sample from the measured receivedamplitude of transmitted excitation sequence trains of pulses fromdiffering wavelengths and locations.

FIG. 24B is a diagram illustrating one method of the invention forintroducing material component, path length, and wavelength basedadjustments to the calculations of the impulse response materialproperties of a sample from the measured received amplitude oftransmitted excitation sequence trains of pulses from differingwavelengths and locations.

FIG. 25 is a set of graphs illustrating example pressure profiles fromdifferent measurement locations responsive to heart beat variation.

FIG. 26 is a graph illustrating the basis of a method for determiningstroke volume and blood flow volume.

FIG. 27 is a flow diagram illustrating the steps of a technique forimproving the sensitivity of measurement of time changing materialproperties in a heterogeneous structure.

FIG. 28 is a set of graphs showing an example impulse response withoutthe presence of a fluorescing dye and with the presence of a fluorescingdye in the measurement sample.

DETAILED DESCRIPTION

References throughout this specification to “one embodiment,” “anembodiment,” “a related embodiment,” or similar language mean that aparticular feature, structure, or characteristic described in connectionwith the referred to “embodiment” is included in at least one embodimentof the present invention. No portion of this disclosure, taken on itsown and/or in connection with a figure, is intended to provide acomplete description of all features of the invention.

In the drawings, like numbers represent the same or similar elementswherever possible. No single drawing is intended to support a completedescription of all features of the invention. In other words, a givendrawing is generally descriptive of only some, and generally not all,features of the invention and may, for the purposes of simplification ofthe drawing, contain not all elements of a particular view or allfeatures that can be presented. The invention may possibly be practicedwithout one or more of the specific features of the invention. Althougha particular detail of an embodiment may not be necessarily shown ineach and every drawing describing such embodiment, the presence of thisdetail in the drawing may be implied unless the context of thedescription requires otherwise. The described single features,structures, or characteristics of the invention may be combined in anysuitable manner in one or more further embodiments.

Moreover, if the schematic flow chart diagram is included thatillustrates the processing flow, the order in which processing steps orparticular methods occur may or may not strictly adhere to the order ofthe corresponding steps shown.

The invention as recited in claims or clauses appended to thisdisclosure is intended to be assessed in light of the disclosure as awhole. While the invention may be described in reference to examples ofan oximeter device, it is understood that, in general the scope of theinvention covers, generally, a device generating and optionallydetecting electromagnetic radiation.

The assessment of arterial blood and tissue oxygen saturation has beenshown to be critical for monitoring, diagnosing, and treating acutecardiovascular deficiencies, such as chronic obstructive pulmonarydisease (COPD), or in extreme cases exsanguination, for example.Further, the time line over which physiological changes occur isindicative of underlying cardiovascular deterioration. While techniquesexist that allow for qualitative monitoring of the level of arterialblood and tissue oxygen saturation, the current medical practice appearsto lack a tool for quantitative monitoring of deterioratingcardiovascular status from the initial cardiovascular event or physicaltrauma through arrival at a medical facility. Recent theoreticaldevelopments have pointed to the potential advantage of multi-wavelengthoximeters, as a way to improve accuracy without extremely large devices.In addition to shortcomings of the commercially availablemulti-wavelength techniques described above, the determination ofwavelengths at which the commercial oximeters operate is somewhatarbitrary and, as proven by practical use, results in the mixing ofspectral data acquired from the patient. As a result, the dataprocessing and retrieval of vital information is unnecessarilycomplicated. Currently used methodologies are therefore quantitativelysuboptimal.

Pulse oximetry is generally used to continuously monitor the arterialblood oxygen saturation of the patient. The “pulse” comes from thetime-varying amount of arterial blood in the tissue during the cardiaccycle. Typical pulse oximetry sensors employ a photodetector and a lightsource producing light at two or more wavelengths to measure the lightthat scatters and/or is transmitted through blood-perfused tissues.Wavelength selection traditionally emphasizes sensitivity to changes inarterial oxygen saturation (SaO₂, SpO₂), with at least one of theemitted wavelengths chosen to fall within the spectral region where theabsorption coefficient of oxygenated hemoglobin (i.e., O₂Hb), ismarkedly different from that of deoxygenated hemoglobin (i.e., HHb). Oneexample of a choice commonly practiced for the multiple sources of lightfor pulse oximetry is to use two light-emitting diodes (LEDs), onegenerating light with a spectrum centered at about 660 nm and the secondLED generating light with a spectrum centered at about 900 nm.

It is also possible to measure other analytes in the arterial blood,such as dyshemoglobins (i.e., MetHb and COHb), with at least oneadditional wavelength for each analyte.

Oximeters may employ vertical cavity surface emitting lasers (VCSELs).The beam of light emanating from a traditional edge-emitting (also knownas in-plane) laser device is usually elliptical in cross section and isoften astigmatic. The angle of divergence of a beam from anedge-emitting laser is significant (on the order of 30 degrees by 10degrees) thereby reducing the efficiency of coupling into an opticalfiber. In contrast, a VCSEL device uses a shaped aperture to provide amore constrained light path, which results in a beam of circular crosssection and of sufficient diameter to make the beam generally lessdivergent, typically at about a 10 degree angle. Operationalcharacteristics of VCSELs are not discussed here in any more detail.

FIG. 1A illustrates schematically a measurement probe module 110 inwhich the laser diodes are an integral part of the probe. The probemodule 110 contains one or more VCSELs 120, optionally spatiallyconfigured in a predetermined pattern or array. Each of the presentVCSELs 120 generates a beam of light substantially centered at aselected wavelength. The probe module 110 also contains a light detector130 configured to receive and detect light emitted by the laser diode(s)120. The light detector 130 may include at least one of any suitabledetector component, such as a photomultiplier tube (PMT), an avalanchephotodiode, and a PIN diode, for example, and in a specific case mayinclude an array of such detectors. The probe module 110 also contains aconnector (not shown) adapted to operably cooperate the probe module 110with a cable (not shown) that interconnects the probe with a monitor(not shown). The probe module 110 is shown to be juxtaposed to thepatient's perfused tissue 140 (shown as a finger, for illustrationpurposes only). The probe module 110 is adapted to operate in thetransmission mode (with the array of laser diodes 120 on one side of thefinger 140 and the light detector 130 on the other side of finger 140).Alternatively, the probe may be adapted to operate in the reflectance orbackscatter mode (where the elements 120 and 130 are positioned on thesame side of the perfused tissue 140). The array of laser diodes 120 mayinclude a diffuser 150 held in place by a mount 160. The diffusion ofthe light output is to ensure that light portion 166 generated by thearray of laser diode light emitting devices 120 is not spatiallyconcentrated and cannot cause harmful effects if mistakenly redirectedto the patient's eyes. The change in the path length of the lightportions 166 as they pass through the arteriolar bed 170 is a figure ofmerit used for determination of the concentration of blood analytes. Theclose spacing of the laser diodes in the array 120 results in theplurality of light beams produced by the laser diodes traversingsubstantially a common path through the arteriolar bed 170, whichimproves the accuracy of the measurements.

Referring to FIG. 1B, the probe module 110 is juxtaposed to thepatient's tissue 170 (shown as an upper arm 180, as an example only) andmultiple light detectors 130 are arranged in a predetermined geometry(e.g., array (square or rectangular), triangular, and pentagon) alongthe same surface of the upper arm 180. Optical paths 184, 185, and 186determined by propagation of light from the light source 120 to thelight detectors 130, 131, and 132 respectively, may traversesubstantially different depths of the tissue 170 and, possibly, thearteriolar bed. The light detectors can generally be positioned atpredetermined distances one from another. In this case, a pulse wavepropagating through the tissue 170 is separately detectable by thespatially separated light detectors 130, 131, and 132. The timedifferences in the pulse-based analyte composition data associated withlight received by these detectors is used to determine variouscharacteristics of the status of the patient such as, for example,cardiac stroke volume and beat-to-beat variation of cardiac strokevolume. Body locations at which the probe is placed may differ,depending on the application.

It is appreciated that sources of monochromatic light are generallydesirable for use as emitters in an oximeter to avoid overlap of thespectral distributions among the emitters and to provide precisecoordination of the spectral bandwidth of a given emitter with respectto the hemoglobin extinction spectral distribution curves. Inparticular, the full-width half-maximum (FWHM) value of an LED bandwidthtypically exceeds at least several tens of nanometers. In contrast, theFWHM of a VCSEL's bandwidth is approximately 1 nm. In addition, the LEDsare typically manufactured without tight quality control of the centralwavelength of the corresponding bandwidth because conventionalapplications of these light sources are in display technologies, not inquantitative measurement systems requires an operation at predeterminedwavelength(s).

As another example, the operational variations in central wavelengths,bandwidths, and shapes of spectral curves among LEDs that in otherapplications may be presumed to be substantially identical, inapplications such as those described below may significantly differ. Tothis end, FIG. 2 shows an example of spectrum 210 of an idealized LED.The curve has a symmetrical bandwidth centered at 600 nm. A second curve215 corresponds to measurements of a nominally 600-nm-centered spectrumof an LED as manufactured. The imperfections of the LED-fabricationprocess result in a red shift of the central wavelength of the spectralband of the curve 215 as compared to that of the idealized curve 210,and in addition, a visibly present spectral output above 700 nm. Thetail end 217 of the spectral emission curve 215 can have a measurableeffect on the received signal from the perfused tissue. The magnitudeand extent of the tail 217 portion of the spectral emission curve 215 isdifficult to predict, measure, and account for as compared to the signaloutput at the nominal central wavelength. Yet another spectral curve 220corresponding to a third LED, the spectrum of which is nominallycentered at a 600 nm wavelength illustrates asymmetry of the spectraloutput with respect to the corresponding central wavelength. Incontradistinction to LEDs, VCSELs can be manufactured with predictableand repeatable spectral outputs, both in terms of central wavelengthsand bandwidths.

Another complicating factor is that the LEDs are mounted in the probemodule, and are juxtaposed to the patient's skin. The LEDs are thereforesubject to significant temperature fluctuations during the operation ofthe probe module, which may cause changes in wavelength output by theLEDs. This effect can introduce a dynamic variation across the LED arrayand cause a measurable source of error. Finally, another source of erroris the so-called “venous prefiltering,” wherein the spectral outputs ofthe LEDs are unevenly and unpredictably attenuated across the span ofgenerated wavelengths by the venous and non-pulsatile arterialcomponents of the blood. Such attenuation of light is a function of theoxygen saturation of the blood and wavelength of the light, varies fromsubject to subject, is temporal in nature, and varies within a givenpatient and across patients. The arterial blood flow is highly variablein the extremities of a patient, where the pulse oximetry readings aretaken. The difference in oxygen saturation between arterial and venouscomponents of the blood can be from as little as less than one percentto greater than twenty-five percent. With greater spectral bandwidth ofthe light source comes greater potential error from this source oferror. Sources of errors are inherent in the LEDs as well as in themethod of placing the LEDs on the patient's appendage to performreadings.

Accordingly, selection of the multiple operating wavelengths, at which asource of light (or a combination of sources of light) of a pulseoximeter should operate, becomes of critical importance. To this end,FIG. 3A presents a graph showing spectral distribution curves for fourblood analytes (oxyhemoglobin, OxyHb; deoxyhemoglobin, DeoxvHb;carboxyhemoglobin, COHb; and methemoglobin, MethHb) and schematicindication of eight wavelengths (610 nm, 620 nm, 630 m, 660 nm, 700 nm,730 nm, 805 nm, and 905 nm) and eight spectral distributions (labeled,collectively, as 300) that respectively correspond to the eight LEDsources employed in a commercially-available hospital oximetry system.It is appreciated that the significant spectral overlap of the LEDs'light outputs (presumably centered at 610 nm, 620 nm, 630 m, 660 nm, 700nm, 730 nm, 805 nm, and 905 nm) results in spectral mixing of theacquired data and at best, merely contribute to signal averaging toreduce background noise. In fact, as established by the methods andalgorithms discussed below, light beams at only four wavelengths chosenfrom the above-mentioned eight wavelengths are sufficient to acquire thepractically-significant majority of data to identify and characterizemultiple hemoglobin species. FIG. 3B illustrates such selection of thefour spectral distributions 320 from the set 300 of FIG. 3A. Using dataacquired with the use of several different light sources operating atdifferent wavelengths allows the user to obtain a solution representingmultiple biochemical species.

While specific examples describing embodiments of the invention arepresented below with applications to optical spectroscopy and, morespecifically, oximetry, quantitative oximetry, or pulse oximetry ofblood samples, embodiments of the invention can be used with other media(for example, gas or fluid) to identify and/or measure differentanalytes (for example, molecules or proteins). Although the particularimplementations of the algorithm of the invention are described usingmatrix equations, the use of continuous equations for determination ofone or more analytes is also within the scope of the invention.Similarly, VCSEL(s) could be replaced with other sources of light suchas, for example, edge emitting semiconductor lasers, or light emittingdiodes. Moreover, embodiments of the invention can be practiced with theuse of different waves, not necessarily optical waves, such as radiowaves or acoustic waves, for example.

When the sample under test includes biological tissue and/or blood, theproperties of the sample determined with light at wavelengths that havebeen selected according to an embodiment of the invention are various.These properties include, for example, concentration of a functional ordysfunctional hemoglobin; glucose; a lipid; a protein; a chromophore; agas (such as percent of oxygen or carbon dioxide); water, and pH level.Other properties that may be measured include normal or abnormal cellcount and normal or abnormal protein levels. The scope of the presentinvention is not limited to biological tissue or the limited propertieslisted. The scope of the invention also includes measuring relativelevels and quantitative amounts of material components in general, aspresent in living as well as non-living samples.

An embodiment of the invention providing for derivative processing mayalso provide for example, one or more of blood oxygen saturation level,blood glucose level, blood protein level, blood lipid level, total bodywater, central blood volume, concentration of fluorescent biomoleculesor dyes, respiratory tidal volume, cardiac stroke volume, andbeat-to-beat variation of cardiac stroke volume. The sample-probingwaves can be varied by including time-based variations of position,angle, intensity, phase, and wavelength of the wave-emitting source. Thesample may also be modified, for example, by introducing a calibratedamount of carbon monoxide (CO) to provide quantitative measurements ofone or more properties of the sample. The scope of the present inventionalso includes recording the measurements over time, such that the timevarying nature of the concentrations can be determined. Propertiesderived from this collection of data include heart rate, heart ratevariability, stoke volume, stoke volume variability, and similarproperties.

Implementations of an algorithm according to the present inventionfacilitate not only optimal optical detection of various species ofhemoglobin that are found in chemical equilibrium with oxygen-saturatedhemoglobin (oxyhemoglobin), such as deoxyhemoglobin (the desaturatedform of hemoglobin) as well as dysfunctional methemoglobin,sulfhemoglobin, and carboxyhemoglobin (the latter being related toexposure to high environmental levels of carbon monoxide), but alsoenablement of determination of the relative percentage of each speciesin blood at any given time. To this end, a set of curves representingspectral distributions of the extinction coefficients (further referredto as absorption curves) for the multiple species of hemoglobin such asthose illustrated in FIG. 3A. Understandably, to resolve multiplespecies from the optically-acquired spectral data one has to solve asystem of equations. As an example, in the case of four hemoglobinspecies, at least four independent equations are required. To derivethese equations, at least four sets of coefficients are needed toincorporate at least four linearly independent optical measurements ofthe spectral characteristics of a blood sample, and therefore at leastfour different optical wavelengths are required. The opticalwavelengths, at which meaningful hemoglobin measurements can beperformed, are in the range of about 450 nm to about 1000 nm. FIG. 3Cillustrates examples of the absorption curves for the species ofhemoglobin (specifically, curves of spectral dependence of extinctioncoefficients) and identifies optical wavelengths (as shown, 611 nm, 650nm, 673 nm, and 857 nm) that are optimal in terms of resolving the fouridentified hemoglobin species with the lowest noise-to-signal ratio(NSR) or, alternatively, the highest signal-to-noise ratio (SNR). It isappreciated that a choice of different absorption curves and associatedwavelengths of operation for the spectrometric device will correspond toa different NSR.

Analysis of the absorption curves of four hemoglobin species (e.g.,oxyhemoglobin, deoxyhemoglobin, methemoglobin, and sulfhemoglobin) showsthat any given species may define a high absorption figure of merit atsome narrow range of frequencies, and a much lower absorption figure ofmerit at a different narrow range of frequencies. In order to measureand solve for the relative or absolute quantity of all four hemoglobinspecies simultaneously, ideally one has to select four linearlyindependent equations with four unknown wavelength-variables. One methodto optimize the selection of wavelengths is to make these equations aslinearly independent of one another as possible in order to yield thelowest NSR for all four measurements. Therefore, the four opticalwavelengths should be chosen such that the respectively correspondingabsorption curves have amplitudes, at these wavelengths, that differfrom one another as much as possible. The wavelengths should be chosensuch that each interaction of electromagnetic radiation (EMR) orelectromagnetic waves (EMWs) with the analyte is captured. If allcomponents are transparent to the EMR or EMWs at a particularwavelength, then that EMR or EMWs have very little interaction with thesample. If any of the analytes are substantially opaque at a givenwavelength, the signal level received by a detector will be difficult tomeasure. Improving the combination of linear independence and theinteraction of the energy with the analyte is captured by improving thepropagation of variance in the system of equations. This is examinedbelow in detail (and, in particular, in reference to Eqs. (29-39)).

According to Ward Cheney and David Kincaid (Numerical Mathematics &Computing, 7th Ed., Brooks Cole; Apr. 27, 2012), for matrix operationsit is possible to compute a “condition number” that provides a gauge ofthe transfer of error in the input to the output of the operation. Asystem with a condition number of one is said to be well conditioned,higher calculated condition numbers indicate progressivelyill-conditioned systems in solving the system Ax=b. If the linear systemis sensitive to perturbations in the elements of A, or to perturbationsof the components of b, then this fact is reflected in A having a largecondition number.

In FIG. 3A, the condition number of the eight-wavelength solution wascalculated to be 30.5; in FIG. 3B, the condition number was 33.1; inFIG. 3C, with optimally chosen wavelengths, the condition number wascalculated to be 19.8; the larger the condition number, the moreill-conditioned the system. Note that in FIG. 3B, the choice ofwavelengths was constrained to the choice of wavelengths included inFIG. 3A. While the condition number is higher in FIG. 3B than in FIG.3A, other figures of merit, including the propagation of error, are muchimproved. The condition number is used here to illustrate one potentialfigure of merit. Note that without the constraint of wavelength choice,the condition number of FIG. 3C is improved in comparison with that ofeither FIG. 3A or FIG. 3B.

According to the idea of the invention, selection of the operatingwavelengths for VCSELS for use in a pulse oximeter or opticalspectrometer is carried out by forming a merit (or cost) function basedon a cost-function approach that is not limited to, for example, thecondition of the mixing matrix or propagation of variance. Some of theconstraint parameters of the cost-function may includewavelength-specific penalty coefficients according to (i)manufacturability of the VCSEL or other components; (ii) bandwidth orspectral extent of the light source output; (iii) a central wavelengthabout which the bandwidth is substantially centered; (iv) transitionareas that result in rapidly changing accuracy with minor changes inoperating center wavelength; (v) differential path length or the averagepath that light of different frequencies traverses in the medium <δ>;(vi) optical path length error derived from the statistical distributionof path lengths associated with, for example, physiologically relevantchanges in the medium; and (vii) analyte absorption equations. Accordingto an algorithm of the invention, the devised overall cost-functionapplies all of these constraints (and additional constraints that may berequired) at the same time in order to optimize the process of definingthe operational wavelengths for the pulse oximeter or, more generally, aspectroscope for measuring spectrally dependent properties of materialcomponents.

The “Manufacturability” constraint parameter defines how easily, orlikely, a VCSEL that operates at a chosen wavelength may bemanufactured. In calculating the overall cost-function, a weighting orscaling (as discussed below, for example in reference to Eq. (40)) isperformed if the chosen wavelength falls outside of themanufacturability constraining limits.

The “Bandwidth” constraint parameter facilitates the simulation ofhypothetical VCSELs or LEDs or even a white light source to account forvariance in the desired bandwidth and to allow for the optional input ofactual light source profiles in-order-to test effectiveness inseparating sharp transition regions.

The “Central wavelength” constraint parameter facilitates testing theeffects of manufacturing tolerances on the performance of a given lightsource. If a particular type of spectral distributions has a largevariance in central wavelengths, this parameter can be included as aconstraint to minimize the error in the solution as a function ofmanufactured center frequency versus the desired central frequency.

The “Transition areas” constraint parameter takes into account therealization that the selection of a wavelength may be such that, in thevicinity of that wavelength, an absorption curve associated with theoperation of the oximeter is rapidly changing (thereby introducingvariability into the processing of acquired data).

The “Differential Path Lengths” constraint parameter accounts for thewavelength dependency of optical scattering coefficients, meaning thatthe optical path length will be different for each light source. Intraditional oximetry, the path lengths are assumed equal and as such arenot part of the consideration. However, in cases characterized by lowoxygen saturation the difference in optical path lengths can result in apoor estimation of oxygen saturation. Accordingly, the algorithm of theinvention takes into account a difference in path lengths correspondingto optical data acquisition at different wavelengths.

The “Interaction” constraint parameter accounts for the wavelengthdependency of differentially excessive transparency and differentiallyexcessive opaqueness of a sample or a material component of the samplesuch that there is either minimal interaction of an EM wave with thesample (transparency), or minimal level of signal received by adetection unit (opaqueness).

A method to minimize error in the solution utilizes the proposedalgorithm. Such method includes simulated annealing, gradient descent(also known as steepest descent), and linear programming. In addition oralternatively, constraints such as full forward and inverse models usinga Monte-Carlo approach are optionally included into or, alternatively,removed from a method depending on the desired system optimization,thereby to create a comprehensive understanding of the measurement ofphysiological variables. For example, constraints or penaltiescorresponding to a particular analyte need not be included in an overallcost-function if that analyte is not present in the medium.

Accordingly, as shown in FIG. 3C, a curve labeled as “penalty”represents a particular subset of a cost function, a scaling factorincorporated into the algorithmic process according to the invention tooptimize the determination of the wavelength choice. This “penalty”function represents, for example, a consideration that lasers (and, inparticular, VCSELS) at some wavelengths are easier to manufacture thanlasers at other wavelengths. The use of the “penalty” function, asdiscussed below, effectively modifies the algorithm by shifting theselection of optimal laser wavelengths to minimize the penalty functionor penalty figure of merit. If the optimization algorithm is run withoutinclusion of the penalty curve, the optimization process arrives atoptimal wavelengths that may be quite different from those chosenwithout consideration of the penalty curve. Because of, for example, newlaser manufacturing approaches, the shape of the penalty curve maychange, and results of a new optimization may be determined byexercising the algorithm with the updated/changed penalty curve. If avalid solution requires the use of a laser source the manufacture ofwhich is particularly expensive, the algorithm selects the correspondingwavelength and a “cost” threshold determines that the solution is notfeasible provided a given design budget. Additional considerations canbe incorporated into the statement of the problem, including otherelements beyond manufacturing constraints on the optical sources, forexample, relative optical responsiveness at varying wavelengths of, ormanufacturing constraints on, the selected photodetectors that comprisethe optical receiver portion of the system.

The wavelength selection process used to separate two distinct elementsof blood (species of hemoglobin, also referred to herein aschromophores) is generally based on an empirical observation thatmeasurements at one or more wavelengths are required to spectrallyseparate one or more chromophores. This observation serves as a basisfor ensuring the minimal spectral mixing of the measured optical signalsrepresenting absorption of probing light by the one or morechromophores. As the number of chromophores of interest increases, sodoes the minimum number of wavelengths required to solve the system ofwavelength-dependent equations. It has previously been documented in theopen scientific literature to be computationally intractable to definemore than a couple of chromophores on a limited subset of thewavelengths of interest, making the solution suboptimal by definition.However, according to the idea of the invention, the algorithm forselection of wavelengths for an oximeter does provide an optimizedsolution (corresponding to laser sources operating at optimally-definedwavelengths or, alternatively, laser sources available commercially andoperating at wavelengths that are close to those defined by anoptimization algorithm) based at least on minimizing the pre-determinedfigure of merit. An example of the pre-determined figure of merit isprovided by propagation of variance from the measurements to thesolution. The algorithm may be used to locate near-optimal solutions, orconstrained solutions, that use particular wavelengths fixed as part ofthe solution. This approach may be advantageously used when a given setof VCSELs at fixed wavelengths are commercially available, and newwavelengths are available only at a considerable incremental cost. Thealgorithm can additionally be used to identify an optimal set ofwavelengths, which can then drive the process of specially designing andmanufacturing VCSELs that operate at central wavelengths correspondingto the wavelengths from the optimal set.

When considering the choice of wavelengths at which the spectrometricdevice, used in a particular application, should preferably operate itis necessary to consider the material components of the sample to bemeasured. Some of these material components—referred to herein aspotentially confounding components—are the material components of thesample the presence of which will affect the spectrum associated withthe measurement(s) of the sample, the presence of which will affect atleast one of the levels of the signal propagated through the sample, andthe level of interaction of the signal with the sample. Suchconsideration requires the use of propagation of variance as opposed tothe condition number. The curve 360 of FIG. 3D illustrates the opticalextinction coefficients of Bilirubin in chloroform as they spectrallyvary over the range from about 240 to about 700 nanometers. (Bilirubinis the yellow breakdown product of normal heme catabolism. Heme is foundin hemoglobin, a principal component of red blood cells. Bilirubin isexcreted in bile and urine.)

As can be empirically determined by the knee of the extinction curve 360Bilirubin, as a material component of the measured medium, has greaterextinction coefficients below 530 nm than above 530 nm. Thischaracteristic illustrates some of the tradeoffs that are available withpotentially confounding analytes. If the concentrations of the analytesare high enough to effect the measurements of the analytes of interest,then the error introduced by the range of potential concentration of theconfounding analyte can be introduced as a spectrally-dependent cost ofthe measurement. This will lead the algorithms for determination ofoptimal wavelength(s) of operation of the spectroscopic device to choosewavelengths that are located in the spectral regions in which a givenconfounding analyte has little or no effect on the measurement.Alternatively, an additional wavelength (to account for a confoundingmaterial component of the sample) can be added to the algorithm, and theconcentration of the confounding material component can be extractedfrom the other analyte determinations.

These two choices (providing a cost curve for interference, or addingthe additional wavelength) are not mutually exclusive. It has beenempirically found that if both are used as inputs to the algorithm whenthe additional wavelength is not needed, two of the wavelengthsproduced, as an output by the algorithm will be identical. Suchoptimization is a result of using propagation of variance as the figureof merit. As was shown in the discussion of FIGS. 3A, 3B and 3C, asmaller number of wavelengths well-chosen, based on the localoptimization of a figure of merit representing operational cost ofemploying the device, improves the performance of the device.

Optical noninvasive sensing methods facilitate detection of additionalmedically important, physiological and biochemical variables besideshemoglobin species. As one example, FIG. 4 offers a graph illustratingabsorption curves for blood glucose, protein, and lipid in the near-IRrange of 1350 nm to 1850 nm. The central wavelengths of spectraldistributions 400 shown in this graph are longer than those in FIGS. 3A,3B, and 3C because the differentiating extinction characteristics of thematerial components (that is, protein, lipid, and glucose) of interestexhibit themselves at longer wavelengths. An absorption curve for wateris provided in FIG. 4 for comparison. It is expected, therefore, thatgiven appropriate operating wavelengths, detection and measurements ofglucose, lipids normal and abnormal proteins, as well as naturalanalytes, as well as manmade materials (for example certain chemicaltags, markers, or dyes) in blood can be effectuated. In addition to thematerial components measured as shown in FIG. 4, other materialcomponents that can be measured include normal or abnormal cell counts.Finally, other concentrations such as percent of a chromophore, percentof a gas (such as oxygen or carbon dioxide), percent water, and pH levelcan also be measured. The above-mentioned material components areimportant as clinical parameters and in some cases as disease markers(for example, diabetes).

The solid vertical lines, marked with wavelength readings (such as 1371nm, 1597 nm, 1679 nm, and 1707 nm) represent the central wavelengthcorresponding to four operating spectral distributions of light selectedfor detection of the above-mentioned species with the use of anoptimization algorithm of the present invention. As the FWHM value for aVCSEL is on the order of one nm, the width of these vertical lines inthe graph approximately represents the spectral distribution of VCSELs'light. The bell-shaped curves 400 (the peaks of which are substantiallycentered on the solid vertical lines), on the other hand, illustrate, incomparison, the broadband spectral distribution of LEDs that would becentered at the same wavelengths. The FWHM value for each LED is on theorder of 50 to 100 nm or at least 50 times that of a VCSEL.

The use of these wavelengths, or a set of wavelengths determined tomatch the constraints described herein, that is an optimal or nearoptimal separation of the major material components, can be used as anon-invasive blood glucose level detector. The techniques described formultiple detectors 130 may be used.

Non-Scattering Absorbing Medium

As the light photons pass through a homogenous, lossy, nonscatteringmedium they are attenuated according to the attenuation coefficientμ_(a)(λ) of the bulk medium (BM) in relation to the traversed distanceδ, the irradiance of light changes in accordance with theBouguer-Beer-Lambert exponential law, the probability of a photon beingabsorbed over a distance is Δδ is Δδμ_(a)(λ), and the mean free pathbetween attenuating interactions l_(t)(λ) is given by l_(t)(λ_(j))=μ_(a)⁻¹(λ_(j)). As individual attenuating components or chromophores of themedium, denoted as C_(i), are considered instead of the properties of abulk medium, the properties of the bulk medium can be expressed as thesummation of all N independent attenuating chromophores:I(λ_(j))=I ₀(λ_(j))exp[−Σ_(i=1) ^(N)μ_(a) ^(C) ^(i) (λ_(j))δ^(C) ^(i)]  (1).

Often μ_(a)(λ) is replaced with the molar extinction or molarattenuation coefficient ϵ(λ_(j)) [M⁻¹ m⁻¹] to more easily allowformulation of the transmission of light as a function of concentrationof the individual attenuating components [C_(j)] that comprise thenonscattering medium. For this substitution, for a given chromophore i,the attenuation A (in units of optical density, OD) and transmissioncoefficients T are related as

$\begin{matrix}{\mspace{20mu}{{A\left( \lambda_{j} \right)} = {{\log_{10}\left( \frac{I_{0}\left( \lambda_{j} \right)}{I\left( \lambda_{j} \right)} \right)} = {{{\epsilon\left( {\lambda_{j},C_{j}} \right)}\left\lbrack C_{i} \right\rbrack}{\delta^{C_{i}}.}}}}} & (2) \\{{A\left( \lambda_{j} \right)} = {{{- \log_{10}}T} = {{- {\log_{10}\left( {\exp\left\lbrack {- {\sum\limits_{i = 1}^{N}{{\mu_{a}^{C_{i}}\left( \lambda_{j} \right)}\delta^{C_{i}}}}} \right\rbrack} \right)}} = {\frac{{\mu_{a}^{C_{i}}\left( \lambda_{j} \right)}\delta^{C_{i}}}{\ln(10)}.}}}} & (3)\end{matrix}$

Accordingly,μ_(a) ^(C) ^(i) (λ_(j))=ln(10)δ^(C) ^(i) ϵ(λ_(j) ,C _(j))[C _(i)]  (4)andA(λ_(j))=Σ_(i=1) ^(N)δ^(C) ^(i) ϵ(λ_(j) ,C _(j))[C _(i)]  (5).

If the optical density coefficient, A, is measured at multiplewavelengths, j=1 . . . M it is expressed in vector and matrix notationsasA=Σ _(j=1) ^(M)Σ_(i=1) ^(N)δ^(C) ^(i) ϵ(λ_(j) ,C _(j))[C _(i)]  (6)andA=ϵ[C]δ  (7).

Equation (7) can be re-written to account for bias G and noise N of themeasurements, to result inA=ϵ[C]δ+G+N  (8).

Given a known set of molar extinction coefficients, the measured valuesof optical density, and assuming that the non-scattering medium ishomogeneous, the solution to Eq. (7) is provided by[C]=ϵ ⁻¹ A(δ^(C) ^(BM) )⁻¹  (9).Scattering Non-Absorbing Medium.

In a non-absorbing medium, light intensity is reduced along a ballisticpath according to the scattering coefficient μ_(s)(Δj) in relation tothe distance δ traversed by the photons. The probability of a photonpassing a distance Δδ without being scattered is Δδμ_(s)(λ_(j)) and,accordingly,I(λ_(j))=I ₀(λ_(j))exp[−Σ_(i=1) ^(M)μ_(s) ^(C) ^(i) (λ_(j))δ^(C) ^(i)]  (10).

The light scattering may be anisotropic, and the directionality of suchscattering is described by the unitless factor g. With the use of apolar angle ϕ (0≤ϕ≤π) and azimuthal angle ψ (0≤ψ≤2π), g can becalculated as

$\begin{matrix}{g = {\frac{\int_{4\pi}{{p(\phi)}{\cos(\phi)}{\sin(\phi)}d\;\phi\; d\;\psi}}{\int_{4\pi}{{p(\phi)}{\sin(\phi)}d\;\phi\; d\;\psi}}.}} & (11)\end{matrix}$

where p(ϕ), which is the probability that a photon is scattered at anangle ϕ defined between the direction of incident light (the incidentphoton's unit vector) r and the scattered photon's unit vector r′, isapproximated by the Henyev-Greenstein function

$\begin{matrix}{{p(\phi)} = \;{\frac{1 - g^{2}}{\left( {1 + g^{2} - {2g\;{\cos(\phi)}}} \right)^{3/2}} = \;{\sum\limits_{q = 0}^{\infty}{\left( {{2q} + 1} \right)g^{q}{{P_{q}\left( {\cos(\phi)} \right)}.}}}}} & (12)\end{matrix}$

The scattering data can be acquired with the use of, for example, anintegrating sphere and then correlated with the results of inverseMonte-Carlo simulations. When g=1, the light is considered to becompletely forward-scattered, and when g=−1, the light is considered tobe completely backward-scattered.

Scattering Absorbing Medium

To account for both scattering and absorbing properties of the medium,the total coefficient of interaction between a photon and the medium canbe expressed asμ_(t)(λ_(j))=μ_(a)(λ_(j))+μ_(s)(λ_(j))  (13).

The unitless albedo factor a is defined as

$\begin{matrix}{{a = \frac{\mu_{a}\left( \lambda_{j} \right)}{\mu_{t}\left( \lambda_{j} \right)}},} & (14)\end{matrix}$and the probability of scattering in such a medium is expressed as

$\begin{matrix}{{p(\phi)} = {{a\;\frac{1 - g^{2}}{\left( {1 + g^{2} - {2\; g\;{\cos(\phi)}}} \right)^{3/2}}} = {a\;{\sum\limits_{q = 0}^{\infty}{\left( {{2q} + 1} \right)g^{q}{{P_{q}\left( {\cos(\phi)} \right)}.}}}}}} & (15)\end{matrix}$

Anisotropy of light scattering in an absorbing medium is accounted forby introducing the reduced scattering coefficientμ′_(s)(λ_(j))=μ_(s)(λ_(j))(1−g)  (16)and the reduced total coefficient of interactionμ′_(t)(λ_(j))=μ_(a)(λ_(j))+μ′_(s)(λ_(j))  (17).

It is conventionally recognized that, as g≥0, the likelihood of forwardscattering increases, thereby decreasing the apparent attenuation oflight along the path of travel. This conventional view does not takeinto account the information about a path that an individual photon oflight has taken to reach a given point and only accounts for the totalflux of light, i.e., <δ>≠δ. Such a conventional approach, therefore,causes errors in the estimates of concentration. In particular, as<δ>»δ, the concentration will be overestimated. Accordingly, the pathlength δ may be multiplied by a scalar constant to account forscattering and for the effective path of light <δ> that is not accountedfor by the geometrical path of light δ.

Scattering medium often defines the so-called diffusion regime where itis assumed that a steady-state fluence rate, Φ, of light propagatingaway from a continuous-wave isotropic point light source of wavelengthλ_(j) can be modeled as a function of radius, Φ(r), in an infinitemedium using the generic diffusion equation in which the effectiveattenuation coefficient, is defined for the diffusion approximation:

$\begin{matrix}{{\Phi(r)} = {\frac{3{\mu_{s}^{\prime}\left( \lambda_{j} \right)}}{4\pi\; r}{{\exp\left( {{\mu_{eff}\left( \lambda_{j} \right)}r} \right)}.}}} & (18)\end{matrix}$

Given these differences between the conventional model and the approachproposed here, it can be seen that information about the path length foran EM wave taken between transmitter and receiver is beneficial towardimproving the accuracy of the derived properties measurements.

Oximetry: Considerations

The major chromophores of a blood sample include the functionalhemoglobins oxyhemoglobin (O₂Hb) and deoxyhemoglobin (Hb), as well asthe dysfunctional hemoglobins (i.e., dyshemoglobins or nonfunctionalhemoglobins), carboxyhemoglobin (COHb), methemoglobin (MetHb),Glycosylated hemoglobin (GHb/Hb λ_(1c)), and sulfhemoglobin (SulfHb).The total concentration of [Hb] is represented by [tHb], which is thesummation of all the fractions according to[tHb]=[O₂Hb]+[Hb]+[COHb]+[MetHb]+[SulfHb]  (21).

The summation of functional hemoglobin is given by[pHb]=[O₂Hb]+[Hb]  (22).

The oxygen saturation is determined as

$\begin{matrix}{{S_{O_{2}} = \frac{\left\lbrack {O_{2}{Hb}} \right\rbrack}{\lbrack{pHb}\rbrack}},} & (23)\end{matrix}$and the oxyhemoglobin fraction is determined as

$\begin{matrix}{F_{O_{2}{Hb}} = {\frac{\left\lbrack {O_{2}{Hb}} \right\rbrack}{\lbrack{tHb}\rbrack}.}} & (24)\end{matrix}$Two-Wavelength Isosbestic Oximetry

If it is assumed that O₂Hb and Hb are the only two functional hemoglobinspecies present in significant concentrations, then only two wavelengthsare required to estimate S_(O2). The two wavelengths λ₁ and λ₂ can bechosen such that for λ₁ the corresponding ϵ (λ₁,O₂Hb) is maximallydifferent from ϵ (λ₁,Hb), and such that λ₂ is an isosbestic point, i.e.that ϵ (λ₂,O₂Hb) exactly equals ϵ(λ₂Hb). With that, the system of twolinear equations requiring a solution is

$\begin{matrix}{\begin{bmatrix}{A\left( \lambda_{1} \right)} \\{A\left( \lambda_{2} \right)}\end{bmatrix} = {{{\begin{bmatrix}{\epsilon\left( {\lambda_{1},{O_{2}{Hb}}} \right)} & {\epsilon\left( {\lambda_{1},{Hb}} \right)} \\{\epsilon\left( {\lambda_{2},{O_{2}{Hb}}} \right)} & {\epsilon\left( {\lambda_{2},{Hb}} \right)}\end{bmatrix}\begin{bmatrix}\left\lbrack {O_{2},{Hb}} \right\rbrack \\\lbrack{Hb}\rbrack\end{bmatrix}}\left\lbrack \delta^{BM} \right\rbrack}.}} & (25)\end{matrix}$

It is assumed that the extinction coefficients are known and theattenuation (optical density) values have been determined empirically.Accordingly, the ratio of the optical density values at the two chosenwavelengths is

$\begin{matrix}{{\frac{A\left( \lambda_{1} \right)}{A\left( \lambda_{2} \right)} = \frac{{\epsilon\mspace{11mu}{\left( {\lambda_{1},{O_{2}{Hb}}} \right)\left\lbrack {O_{2}{Hb}} \right\rbrack}} - {\epsilon\mspace{11mu}{\left( {\lambda_{1},{Hb}} \right)\left\lbrack {O_{2}{Hb}} \right\rbrack}} + {\epsilon\mspace{11mu}{\left( {\lambda_{1},{Hb}} \right)\lbrack{pHb}\rbrack}}}{\epsilon\mspace{11mu}{\left( {\lambda_{2},{O_{2}{Hb}}} \right)\lbrack{pHb}\rbrack}}},} & (26)\end{matrix}$and the oxygen saturation ratio is determined, therefore, from

$\begin{matrix}{S_{O_{2}} = {\frac{\left\lbrack {O_{2}{Hb}} \right\rbrack}{\lbrack{pHb}\rbrack} = {\frac{{\frac{A\left( \lambda_{1} \right)}{A\left( \lambda_{2} \right)}\epsilon\mspace{11mu}\left( {\lambda_{2},{O_{2}{Hb}}} \right)} - {\epsilon\mspace{11mu}\left( {\lambda_{1},{Hb}} \right)}}{{\epsilon\mspace{11mu}\left( {\lambda_{1},{O_{2}{Hb}}} \right)} - {\epsilon\mspace{11mu}\left( {\lambda_{1},{Hb}} \right)}}.}}} & (27)\end{matrix}$Generalized Two-Wavelength Oximetry

The derivations provided above demonstrate a solution of the system ofequations when an isosbestic point is chosen for one of the λ's.Generally, one of the primary reasons behind using an isosbestic pointis the simplification of such a solution to one that does not require alarge amount of computational effort. Generally, however, if the abovelinear equations (25) are simply viewed as an unmixing matrix of twosuperimposed signals with known, unique mixing coefficients ϵ, any λ₁and λ₂ can be selected to estimate the oxygen saturation value—notnecessarily the isosbestic wavelength points. Indeed, for twoindependent wavelengths λ₁ and λ₂ (or two non-overlapping spectralbandwidths respectively centered at λ₁ and λ₂), the Eq. (7) can bere-written as

$\begin{matrix}{\begin{bmatrix}{A\left( \lambda_{1} \right)} \\{A\left( \lambda_{2} \right)}\end{bmatrix} = {{{\begin{bmatrix}1.0 & 0.0 \\0.0 & 1.0\end{bmatrix}\begin{bmatrix}\left\lbrack {O_{2}{Hb}} \right\rbrack \\\lbrack{Hb}\rbrack\end{bmatrix}}\left\lbrack \delta^{BM} \right\rbrack}.}} & (28)\end{matrix}$

In this case, the coefficient matrix defines mixing ϵ factors that aremutually orthogonal and unitary. The corresponding noise amplificationand propagation of variance are given as:

For noise amplification:[g(λ₁),g(λ₂)]=diag[(ϵ^(H)Ψ⁻¹ϵ)⁻¹](diag[ϵ^(H)Ψϵ])  (29a),and[g(λ₁),g(λ₂)]=[1,1]  (29b).

For propagation of variance:[σ(λ₁),σ(λ₂)]=√{square root over(Σ[(ϵ^(H)ϵ)⁻¹(ϵ^(H))var(Noise)])}  (30a),and[σ(λ₁),σ(λ₂)]=[1,1]  (30b).

In contradistinction, for two spectral bandwidths that do overlap afigure of merit needs to be defined, which takes into account suchspectral mixing/overlap, in order to optimally select λ₁ and λ₂ andamplitude. According to an embodiment of the invention, the choice oftwo wavelengths for a spectrally overlapping case includes thedetermination of noise amplification and/or propagation of variancetechniques. For example, with numerical coefficients chosen that are notlinearly independent and therefore have a condition number greater thanone and that result in mixing of the information, the solution becomesmore complicated because A(λ₁) now comprises signal components from O₂Hband Hb. Similarly A(λ₂) now includes signal components from O₂Hb and Hb;whereas A(λ₁) should contain signal components solely from O₂Hb, andA(λ₂) should contain signal components solely from Hb. This contaminatedor impure situation is thereby described as a mixing matrix ornon-diagonal matrix. In the case of spectrally overlapping bandwidthscorresponding to

$\begin{matrix}{{\begin{bmatrix}{A\left( \lambda_{1} \right)} \\{A\left( \lambda_{2} \right)}\end{bmatrix} = {{\begin{bmatrix}0.5 & 0.03 \\0.2 & 0.4\end{bmatrix}\begin{bmatrix}\left\lbrack {O_{2}{Hb}} \right\rbrack \\\lbrack{Hb}\rbrack\end{bmatrix}}\left\lbrack \delta^{BM} \right\rbrack}},} & (31)\end{matrix}$the noise amplification is observed to increase as the mixing vectorsbecome less orthogonal.

The propagation of variance equation should now account for thedecreased orthogonality of the mixing vectors and the decreasedϵ-sensitivity, not otherwise predicted with noise amplification. In thepropagation of variance, any decrease in ϵ sensitivity causes anindependent decrease in SNR, because there is less interaction betweenlight and the analyte for a given source-to-detector distance associatedwith the data acquired at each wavelength. However, ϵ, the interactionpenalty, is not observed with noise amplification, and is a veryimportant distinction between the two techniques. Because of thisdifference, noise amplification is the less preferred technique comparedto the propagation of variance technique.

Here, for the above example, the interaction is demonstrated:

For noise amplification:[g(λ₁),g(λ₂)]=diag[(ϵ^(H)Ψ⁻¹ϵ)⁻¹](diag[ϵ^(H)Ψϵ])  (32a),and[g(λ₁),g(λ₂)]=[1.1135,1.1135]  (32b).

For propagation of variance:[σ(λ₁),σ(λ₂)]=√{square root over(Σ[(ϵ^(H)ϵ)⁻¹(ϵ^(H))var(Noise)])}  (33a),and[σ(λ₁),σ(λ₂)]=[2.0676,2.7759](33b).

The above formulations do not account for the absolutewavelength-dependent ϵ sensitivity and therefore a more general form ofnoise amplification and propagation of variance parameter may be needed;one skilled in the art can understand the additional constraints thatapply for analytes. Unlike the noise amplification approach, the use ofthe propagation of variance technique for estimation of the oxygensaturation level S_(O) ₂ depends on both the orthogonality of the ϵsensitivities and the relative sensitivity of ϵ at each of thewavelengths. The optimal solution occurs when the signals aresubstantially not mixed spectrally (which corresponds to equation (28)),resulting in two independent equations and the relative ϵ sensitivities,each of which is equal to 1. This is important because noiseamplification and propagation of variance techniques provide an unbiasedestimate of the optimal λj's that should be selected for a given set ofchromophores in the experiment. One skilled in the art will observe thatboth of these approaches may be further adapted to the specificapplication.

Multiwavelength Oximetry

According to embodiments of the invention, the above-discussed methodsare further extended to a case of multiwavelength oximetry (employingmore than two operating wavelengths). Another approach may include afull-matrix solution of the equation such as equation (7) by analogywith that discussed in reference to equations (31, 32a, 32b, 33a, and33b). It is appreciated that, with the increase of the number ofoperational wavelengths, the overall variance is increased as well. Inthe example shown below, five operational wavelengths λ₁ through λ₅ areused for determination of only two species (oxyhemoglobin anddeoxyhemoglobin):

$\begin{matrix}{\begin{bmatrix}{A\left( \lambda_{1} \right)} \\{A\left( \lambda_{2} \right)} \\{A\left( \lambda_{3} \right)} \\{A\left( \lambda_{4} \right)} \\{A\left( \lambda_{5} \right)}\end{bmatrix} = {{{\begin{bmatrix}0.5 & 0.03 \\0.2 & 0.4 \\0.01 & 0.03 \\0.18 & 0.23 \\0.21 & 0.13\end{bmatrix}\begin{bmatrix}\left\lbrack {O_{2}{Hb}} \right\rbrack \\\lbrack{Hb}\rbrack\end{bmatrix}}\left\lbrack \delta^{BM} \right\rbrack}.}} & (34)\end{matrix}$

In this example, accordingly, for noise amplification:[g(λ₁),g(λ₂)]=diag[(ϵ^(H)Ψ⁻¹ϵ)⁻¹](diag[ϵ^(H)Ψϵ])  (35a)and[g(λ₁),g(λ₂)]=[1.2155,1.2155]  (35b).

For propagation of variance:[σ(λ₁),ø(λ₂)]=√{square root over(Σ[(ϵ^(H)ϵ)⁻¹(ϵ^(H))var(Noise)])}  (36a)and[σ(λ₁),σ(λ₂)]=[1.8787,2.5257]  (36b).

It is important that the wavelengths selected in multi-wavelengthoximetry be optimized to improve the orthogonality and analytesensitivity of the mixing equation. Mere addition of supplementarywavelengths can have the effect of decreasing the SNR in the resultingestimate of S_(O2). Further, according to the invention, selection ofadditional (supplementary) wavelengths is made dependent on the numberof chromophores (elements of blood) to be measured. Specifically, themultiple wavelengths selected for determination of two chromophores(such as O₂Hb and Hb, for example) may be different from the same numberof multiple wavelengths selected for a five-chromophore oximetrymeasurement. The result of selecting multiple wavelengths according tothe idea of the invention allows for the detection of dyshemoglobins,for example, and for improved sensitivity of the pulse oximeter orspectrometer device without adversely affecting the noise figure in theestimates of concentration of the functional hemoglobins. In onespecific example:

$\begin{matrix}{{\begin{bmatrix}{A\left( \lambda_{1} \right)} \\{A\left( \lambda_{2} \right)} \\{A\left( \lambda_{3} \right)} \\{A\left( \lambda_{4} \right)} \\{A\left( \lambda_{5} \right)}\end{bmatrix} = {{\begin{bmatrix}{\epsilon\left( {\lambda_{1},C_{1}} \right)} \\{\epsilon\left( {\lambda_{2},C_{2}} \right)} \\{\epsilon\left( {\lambda_{3},C_{3}} \right)} \\{\epsilon\left( {\lambda_{4},C_{4}} \right)} \\{\epsilon\left( {\lambda_{5},C_{5}} \right)}\end{bmatrix}\begin{bmatrix}\left\lbrack {O_{2}{Hb}} \right\rbrack \\\lbrack{Hb}\rbrack \\\left\lbrack {{CO}{Hb}} \right\rbrack \\\lbrack{MetHb}\rbrack \\\lbrack{rHb}\rbrack\end{bmatrix}}\begin{bmatrix}\delta^{C_{1}} \\\delta^{C_{2}} \\\delta^{C_{3}} \\\delta^{C_{4}} \\\delta^{C_{5}}\end{bmatrix}}},} & \left( {37a} \right)\end{matrix}$

where equation 37a is derived from equation 7 for the specific exampleusing five exemplary material components (O₂Hb, Hb, COHb, MetHb, andrHb) of interest. Exemplary extinction coefficients are providedyielding:

$\begin{matrix}{{\begin{bmatrix}{A\left( \lambda_{1} \right)} \\{A\left( \lambda_{2} \right)} \\{A\left( \lambda_{3} \right)} \\{A\left( \lambda_{4} \right)} \\{A\left( \lambda_{5} \right)}\end{bmatrix} = {{\begin{bmatrix}0.50 & 0.03 & 0.01 & 0.01 & 0.01 \\0.20 & 0.4 & 0.01 & 0.01 & 0.01 \\0.01 & 0.04 & 0.40 & 0.07 & 0.01 \\0.18 & 0.23 & 0.01 & 0.60 & 0.01 \\0.21 & 0.13 & 0.01 & 0.07 & 0.40\end{bmatrix}\begin{bmatrix}\left\lbrack {O_{2}{Hb}} \right\rbrack \\\lbrack{Hb}\rbrack \\\left\lbrack {{CO}{Hb}} \right\rbrack \\\lbrack{MetHb}\rbrack \\\lbrack{rHb}\rbrack\end{bmatrix}}\left\lbrack \delta^{BM} \right\rbrack}},} & \left( {37b} \right)\end{matrix}$

Similarly, as shown in equation 9, it is also desirable to be able toderive the concentration of the material components of interest givenattenuations for each waveform transmitted. Using matrix operations onequation 37a, it is shown that the concentrations of each materialcomponent can be derived from the molar extinction coefficientsϵ(λ_(i),C_(i)), attenuations A(λ_(i)) of the waveforms, and thetraversed distances δ^(Ci). This form of the equation is shown inequation 37c for using the distance associated with bulk media.

$\begin{matrix}{{{{\begin{bmatrix}{\epsilon\left( {\lambda_{1},C_{1}} \right)} \\{\epsilon\left( {\lambda_{2},C_{2}} \right)} \\{\epsilon\left( {\lambda_{3},C_{3}} \right)} \\{\epsilon\left( {\lambda_{4},C_{4}} \right)} \\{\epsilon\left( {\lambda_{5},C_{5}} \right)}\end{bmatrix}\begin{bmatrix}{A\left( \lambda_{1} \right)} \\{A\left( \lambda_{2} \right)} \\{A\left( \lambda_{3} \right)} \\{A\left( \lambda_{4} \right)} \\{A\left( \lambda_{5} \right)}\end{bmatrix}}\left\lbrack \delta^{BM} \right\rbrack}^{- 1} = \begin{bmatrix}\left\lbrack {O_{2}{Hb}} \right\rbrack \\\lbrack{Hb}\rbrack \\\left\lbrack {{CO}{Hb}} \right\rbrack \\\lbrack{MetHb}\rbrack \\\lbrack{rHb}\rbrack\end{bmatrix}},} & \left( {37c} \right)\end{matrix}$

equation 37d shows the same equation (37c) using the inverted exemplaryextinction coefficients shown previously in equation 37b.

$\begin{matrix}{{{{{{\begin{bmatrix}2.08 & {- 0.12} & {- 0.05} & {- 0.02} & {- 0.05} \\{- 1.02} & 2.60 & {- 0.04} & {- 0.02} & {- 0.04} \\0.11 & {- 0.08} & 2.51 & {- 0.29} & {- 0.06} \\{- 0.22} & {- 0.95} & {- 0.01} & 1.69 & {- 0.01} \\{- 0.72} & {- 0.61} & {- 0.02} & 0.27 & 2.54\end{bmatrix}\begin{bmatrix}{A\left( \lambda_{1} \right)} \\{A\left( \lambda_{2} \right)} \\{A\left( \lambda_{3} \right)} \\{A\left( \lambda_{4} \right)} \\{A\left( \lambda_{5} \right)}\end{bmatrix}}\left\lbrack \delta^{BM} \right\rbrack}^{- 1} = {\quad\quad}}\quad}\quad}{\quad{\begin{bmatrix}\left\lbrack {O_{2}{Hb}} \right\rbrack \\\lbrack{Hb}\rbrack \\\left\lbrack {{CO}{Hb}} \right\rbrack \\\lbrack{MetHb}\rbrack \\\lbrack{rHb}\rbrack\end{bmatrix},}}} & \left( {37d} \right)\end{matrix}$

Referring to equation 37b, noise should also be taken intoconsideration. For noise amplification:g(λ₁),g(λ₂)]=diag[(ϵ^(H)Ψ⁻¹ϵ)⁻¹](diag[ϵ^(H)Ψϵ])  (38a),and[g(λ₁),g(λ₂)]=[1.3546,1.3440]  (38b).

For propagation of variance:[σ(λ₁),σ(λ₂)]=√{square root over(Σ[(ϵ^(H)ϵ)⁻¹(ϵ^(H))var(Noise)])}  (39a)and[σ(λ₁),σ(λ₂)]=[2.0936,2.7927]  (39b).Weight Function or Parameter

It is appreciated that, given the derivation of equation (7), additionalwavelength-specific penalties or weighting factors can be incorporatedvia a two-dimensional (2D) or vector (1D) matrix W representing acost-function that incorporates, for example, noise penalties. When thespectral channels of the oximeter (i.e., the operational wavelengths andassociated spectral bandwidths) are independent (non-overlapping), W isthe identity matrix W₁. Alternatively the W=W_(sep) matrix can includepenalties associated with the wavelength separation (which relates theseparation to the increased difference in δ or the differential pathlength factor <δ> as the operational wavelengths are spaced furtherapart). In yet another implementation, an alternative or additionalW=W_(cost) weighting matrix can include penalties for the complexity orcost of manufacture of the source of light operating at a particularwavelength. Further, all or at least some of such weighting factors canbe used in series to make the optimization of wavelength selection morepractical. Accordingly, to incorporate a weighting function, the generalmatrix equation (7) can be re-written, for example, from[g(λ₁)W(λ₁),g(λ₂)W(λ₂)] and [σ(λ₁)W(λ₁),σ(λ₂)W(λ₂)]  (40a),to,[g(λ₁)W _(I) W _(cost)(λ₁),g(λ₂)W _(I) W _(cost)(λ₂)]  (40b)and[σ(λ₁)W _(I) W _(cost)(λ₁),σ(λ₂)W _(I) W _(cost)(λ₂)]  (40c).Pulse Oximetry

According to one implementation of the invention, the oximetry dataprocessing method that takes into account the “pulse” nature of thepulse oximetry data-acquisition (such as, for example, (i) the presenceand/or parameters of the arterial pulse wave and (ii) the temporalpoint-spread function or impulse response associated with the tissuebeing measured, including the time-ordered sampling of the effect of thetransfer function of the tissue on light passing through the object)determines the temporal dependence of light attenuation by the tissue.For example, in the case when the two chosen operationalwavelengths/bandwidths do not overlap, the equation (7) may expanded toexplicitly include the tissue and venous component factors:

$\begin{matrix}{\begin{bmatrix}{A\left( {\lambda_{1},t} \right)} \\{A\left( {\lambda_{2},t} \right)}\end{bmatrix} = {{\begin{bmatrix}{\epsilon\mspace{11mu}\left( {\lambda_{1},{O_{2}{Hb}}} \right)} & {\epsilon\mspace{11mu}\left( {\lambda_{1},{Hb}} \right)} & {\epsilon\mspace{11mu}\left( {\lambda_{2},{tissue}} \right)} \\{\epsilon\mspace{11mu}\left( {\lambda_{2},{O_{2}{Hb}}} \right)} & {\epsilon\mspace{11mu}\left( {\lambda_{2},{Hb}} \right)} & {\epsilon\mspace{11mu}\left( {\lambda_{2},{tissue}} \right)}\end{bmatrix}\begin{bmatrix}\left\lbrack {O_{2}{{Hb}(t)}} \right\rbrack \\\left\lbrack {{Hb}(t)} \right\rbrack \\\left\lbrack {{tissue}{\mspace{11mu}\;}(t)} \right\rbrack\end{bmatrix}}{\quad{\begin{bmatrix}{\delta^{O_{2}{Hb}}(t)} \\{\delta^{Hb}(t)} \\{\delta^{tissue}(t)}\end{bmatrix}.}}}} & (41)\end{matrix}$

In this example, and further considering the time-derivative of theequation (41) expressed as

$\begin{matrix}{{\frac{d\left\lbrack {A(t)} \right\rbrack}{dt} = {\frac{d}{dt}\begin{bmatrix}{A\left( {\lambda_{1},t} \right)} \\{A\left( {\lambda_{2},t} \right)}\end{bmatrix}}},} & (42)\end{matrix}$the following assumptions are made:

-   -   (i) for a given chromophore, a portion of the absorption        spectrum (of such a chromophore) that remains unchanged with        time is set to zero by the above time derivative;    -   (ii) the effective distance traveled by light is substantially        the same for all chromophores and equals Δ<δ>, or the change in        distance caused by the arterial pulse, and    -   (iii) only two chromophores are being considered (which is an        optional assumption specific only to this example).        The cancellation of the change in path length due to the        arterial pulse-wave from the equation is achieved by, for        example, defining a ratio of the equations (39) corresponding to        two different wavelengths:

$\begin{matrix}{{\frac{R_{\lambda_{1}}}{R_{\lambda_{2}}} = {\frac{\frac{d\left\lbrack {A\left( {\lambda_{1},t} \right)} \right\rbrack}{dt}}{\frac{d\left\lbrack {A\left( {\lambda_{2},t} \right)} \right\rbrack}{dt}} = {\frac{{\epsilon\mspace{11mu}\left( {\lambda_{1},{O_{2}{Hb}}} \right){\Delta\left\lbrack {O_{2}{Hb}} \right\rbrack}} + {\epsilon\mspace{11mu}\left( {\lambda_{1},{Hb}} \right){e\lbrack{Hb}\rbrack}}}{{\epsilon\mspace{11mu}\left( {\lambda_{2},{O_{2}{Hb}}} \right){\Delta\left\lbrack {O_{2}{Hb}} \right\rbrack}} + {\epsilon\mspace{11mu}\left( {\lambda_{2},{Hb}} \right){e\lbrack{Hb}\rbrack}}}\frac{\Delta\left\langle \delta \right\rangle}{\Delta\left\langle \delta \right\rangle}}}},} & \left( {43a} \right) \\{\mspace{79mu}{{\Delta\lbrack{pHb}\rbrack} = {{{\Delta\left\lbrack {O_{2}{Hb}} \right\rbrack} + {{\Delta\lbrack{Hb}\rbrack}{\Delta\lbrack{pHb}\rbrack}}} = {{\Delta\left\lbrack {O_{2}{Hb}} \right\rbrack} + {{\Delta\lbrack{Hb}\rbrack}.}}}}} & \left( {43b} \right)\end{matrix}$

The solution for the oxygen saturation level for pulse oximetry can bedetermined under these assumptions as

$\begin{matrix}{S_{p,O_{2}} = {\frac{\Delta\left\lbrack {O_{2}{Hb}} \right\rbrack}{\Delta\lbrack{pHb}\rbrack} = {\frac{{\frac{R_{\lambda_{1}}}{R_{\lambda_{2}}}\epsilon\;\left( {\lambda_{2},{Hb}} \right)} - {\epsilon\mspace{11mu}\left( {\lambda_{1},{Hb}} \right)}}{{\epsilon\mspace{11mu}\left( {\lambda_{1},{O_{2}{Hb}}} \right)} - {\epsilon\mspace{11mu}\left( {\lambda_{1},{Hb}} \right)} + {\frac{R_{\lambda_{1}}}{R_{\lambda_{2}}}\left( {{\epsilon\mspace{11mu}\left( {\lambda_{2},{Hb}} \right)} - {\epsilon\left( {\lambda_{2},{O_{2}{Hb}}} \right)}} \right)}}.}}} & (44)\end{matrix}$

Given the similarity between the pulse oximetry solution (44) and thestandard oximetry solution (27), one skilled in the art can appreciatethat it is feasible that a wavelength selected to be optimal forstandard oximetry may also be optimal for pulse oximetry.

Qualitative and Quantitative Pulse Oximetry

The major differences between qualitative pulse oximetry andquantitative oximetry are: A) pulse oximetry measures the ratio ofoxyhemoglobin to deoxyhemoglobin by first assuming that the ratio of thewavelength-specific optical path lengths is unity, or alternately that aratio factor is well known, thus removing unknown variables from thesystem of equations; B) pulse oximetry assumes that thewavelength-specific time-varying arterial pulse measured by the pulseoximeter is caused by increased path length of light through thearteries; C) pulse oximetry assumes that the ratio of time-varyingsignals at different wavelengths is related linearly, or by a low-orderpolynomial, to invasive blood gas measurements (i.e., by medicallaboratory measurement of the oxygen saturation of a sample of blooddrawn from the patient) to calibrate the device to percent saturation orto the ratio of oxyhemoglobin to the combined amount of oxyhemoglobinand deoxyhemoglobin [i.e.,oxyhemoglobin/(oxyhemoglobin+deoxyhemoglobin)]; D) pulse oximetry relieson the invalid Beer-Lambert law; E) quantitative oximetry assumes thatlight diffuses down a concentration gradient from a source, rather thanthe Beer-Lambert law's known-incorrect assumption of a direct opticalpath without scattering between the optical source and the photodetector; F) quantitative oximetry attempts to account for the scatterin tissue which alters the optical path length from the geometric, orshortest path length assumed by the Beer-Lambert law, to the averagepath length which is greater than the geometric path; and G)quantitative oximetry does not rely on the arterial pulse and insteadmeasures the absorption at specific wavelengths that directly relate tothe concentration of oxyhemoglobin and deoxyhemoglobin.

In tissue, the path that light follows from the source to the detectoris not directly along the minimal geometric path. Instead, the light ishighly scattered, resulting in a net path length longer than thegeometric path length. In addition, the net path length is wavelengthspecific. In algorithms of the related art, see for example, Ayogi et al(Int'l Anesthesia Res. Soc., v. 105, No. 6. December 2007), the opticalpath length factors for the wavelengths are assumed to be similar,allowing them to be divided out. In addition, the differential pathlength factor is specific to anatomy, e.g. finger versus the scalp orpinnae of the ear, while normal anatomical variation, e.g. fingerthickness, makes the differential path length factor highly inaccurate.In reality requiring the path length to be the same for all wavelengthsand for the placement on patients and anatomy of patients to be exactlythe same is unlikely; therefore, measuring the actual optical pathlength for each wavelength is advantageous.

Examples of Algorithms for Solving Matrix Equations

Simulated Annealing

In one implementation, the technique referred to as simulated annealingis used to determine the optimum selection of wavelengths. Simulatedannealing (SA) is a general meta-heuristic for locating a goodapproximation to the global optimum of a given function in a largesearch space. For the wavelength choice task, simulated annealing ismore tractable than exhaustive enumeration. The name of the method comesfrom annealing in metallurgy, a technique involving heating andcontrolled cooling of a material to increase the size of its crystalsand reduce their defects; both are attributes of the material thatdepend on its thermodynamic free energy. Heating and cooling thematerial affects both the temperature and the thermodynamic free energy.As illustrated in the flow diagram of FIG. 5, in this solution space,temperature is an analog for both the probable range of change selectedin ‘permute the solution’ 510 between successive explorations of thesolution space, and the likelihood P that that a “worse” solution willbe chosen at ‘evaluate choice’ P(E_(OLD), E_(NEW), Temp) 515. The analogfor the thermodynamic free energy of metallurgy is the figure of meritE_(WAVESET) used to evaluate the choice of wavelengths. Cooling isimplemented in the Simulated Annealing algorithm as a slow decrease inthe probability of accepting worse solutions as the algorithm exploresthe solution space. Accepting worse solutions is a fundamental propertyof this heuristic because it avoids the traps of locally optimalsolutions that are worse than the global optimum.

With reference to FIG. 5, a simulated annealing optimization algorithminitializes various data structures 505, including a starting point forthe solution (selection of wavelengths), the equation to calculate theworthiness of the solution E_(WAVESET), the present value of thegoodness E_(OLD) (referred to here as the energy of the solution, whereenergy is desired to be minimized), and a setting of the current“temperature” (Temp). Additional initialization may take place withoutloss of generality. The chosen solution is permuted (the selection ofwavelengths is updated) at random from the possible solutions. The rangeof permutations is governed by the use of the “temperature” of thecomputation, said “temperature” to be reduced over the course ofsuccessive runs such that the range of possible solutions to be examinedis updated in a stepwise fashion (for example, decreased) over time, asdiscussed below. As each permuted solution is evaluated, 515, if theobjective function, E_(WAVESET), yields a better number (which occurs ifthe “energy” of the permuted solution state is less than the previousstate), then with probability P(E_(OLD), E_(NEW), Temp) the presentstate is updated to reflect the permuted state, step 520. Proceedingfrom this step, the evaluation to determine if permutations at thisenergy level have completed is performed. A number of techniques may beused for this evaluation, including but not limited to a countassociated with a particular “temperature” step, or an evaluation ofchange over the last selected number of evaluation runs. If there are nomore perturbations to be performed at this “temperature,” the“temperature” is updated 530. If there are no more “temperature”selections to be performed, i.e., if the computation is complete 535,then the result is returned 540 and the algorithm finishes 545. If theevaluation 535 so determines, computation will continue at 510,producing an additional permutation at the new “temperature.”

The properties of ‘evaluate choice’ P(E_(OLD), E_(NEW), Temp) 515 aresuch that as the temperature (Temp) decreases, the likelihood of thelower of E_(OLD) and E_(NEW) will be chosen. As described above theobjective function, E_(WAVESET)(λ₁, . . . , λ_(N)) is based on thewavelengths selected (λ₁, . . . , λ_(N)). Processing the wavelengthchoices through the spectrometric graphs determines the conversionmatrix to be used from received signal to material composition. Thematrix thus formed is evaluated for propagation of variance, and thisresult is weighted by the various cost functions related to the chosenwavelengths and the relation between the wavelengths. This operationproduces the figure of merit for the particular choice of wavelengths.See Equations 40a-40c.

With reference to FIG. 6, it is advantageous to verify that the statechosen by the simulated annealing at step 545 is the optimal choice. Ithas been found that the confidence in the optimality of the chosen statecan be increased by performing multiple trials of the sequenceillustrated in FIG. 5 (steps 500-545). The validation of the simulatedannealing solution includes choosing a “temperature cooling” profile atstep 655 and running multiple tasks for this “temperature profile.”These tasks are distributed in step 660 and the results are collected atstep 665. If results collected based on multiple runs are substantiallythe same, as determined at step 670, then the confidence is high thatthe optimal choice has been returned, 680, and the algorithm ends atstep 690. If the results from multiple runs differ, based on comparisoncarried out at step 670, then “select slower cooling” 675 engages aslower cooling profile than had been used previously at “DistributeTasks” 660. The operation of evaluating the tasks described in referenceto steps 660-670 is repeated until the cohort of solutions convergesbased on comparison performed at step 670.

The choice of wavelength(s) and the weighting functions (or constraintparameters) representing partial or full operational cost of employingthe device and used in the evaluation may include finite or discretevariables and continuous variables. The random perturbation of thewavelengths in the solution state can be augmented by a number oftechniques, including those known in the art but not limited to gradientdecent (i.e., if a change improved the result last time, make a furtherchange in the same direction) and statistical sampling.

Additional Algorithms

Convex optimization, as described by Stephen Boyd and LievenVandenberghe, (in Convex Optimization, Cambridge University Press; Mar.8, 2004) and linear programming as described by R. Fletcher (inPractical Methods of Optimization, Wiley, 1987) are two additionalclasses of techniques that can be used to find improved solutions to thequantified problem as stated above.

Structures, Systems, and Methods

A combination of structures, systems, and methods (collectively referredto as techniques herein) can be used to improve accuracy of themeasurements in the presence of noise, to increase sensitivity of themeasurement system (for example, with respect to the measurement of verysmall component percentages), to discriminate between different materialcomponents, and to provide and improve the ability to discern changes inthe material makeup of a measured sample over time and location. Thevariation in the measurements to be made can come from many sources, anda combination of the techniques is needed to address the combination ofthe variation sources identified below.

Causes of Measurement Signal Variation and Degradation

Causes of signal variation or degradation of accuracy during ameasurement include heterogeneity of the sample, the presence ofconfounding materials, electronic noise, non-linearity of theelectronics and optics, ambient radiation, randomness, which includesvariations in signal path, quantization at translation points, quantumeffects of very low signal levels, fluorescence within the sample, andnon-linear component interaction. Each of these is discussed below,although this list should not be considered exhaustive.

Heterogeneity of the Sample.

Heterogeneity of the sample being measured causes variations inmeasurement results depending on the paths taken by radiation throughthe sample. If the sample is not perfectly homogenous, then thevariations in composition can be considered a noise element of themeasurement. Spatial localization of the measurements, eithermechanically or using fine-grained temporal characterization of thesample response, can be used to limit the noise caused by suchinhomogeneity. Averaging the results over numerous measurementspartially alleviates this problem, but it is often precise informationon the variation over time and space that is useful diagnostically. Forthis reason, several of the techniques described below are beneficiallyemployed.

Referring to FIG. 7, an arm 700 is schematically illustrated, with twocross-sectional views 710 (taken along the plane J-J′) and 720 (takenalong the plane I-I′), where 710 is a notional cross section taken closeto the wrist, and 720 is a notional cross section in the vicinity of anelbow. FIG. 7 is intended to illustrate the heterogeneous makeup of thesample used when an embodiment of a spectrometric device according tothe invention is juxtaposed with a human. Different parts of the samplewill modify the signal differently as it passed through the sample. Thebone structure 730 is also shown in the two cross-sectional views 710and 720. The cross sectional views include muscle tissue 740, not shownin the main view for purposes of clarity. The major vein system 750 andarterial system 770 are schematically illustrated to show the variationof the tissue throughout the sample, and to show that the size of theveins and arteries varies, as does their proximity to each other. Forthis reason, measurements that are configured to localize the responseto the tissue of interest provide a better result. Also pictured in thecross sectional views are notional locations of the tendons.

Confounding Materials.

While the concentration of particular material components in a tissuesample may be of interest, it will also be important to assess thepotential presence of additional components of the sample that mayconfound the measurements of the material components of interest. Aparticular confounding component may vary in concentration ordistribution. Potentially confounding materials are spectrometricallyactive, that is, their effects on the propagation of a signal atdifferent wavelengths vary. The description of bilirubin offered inreference to FIG. 3D above provides an example of a potentiallyconfounding material and several appropriate treatments of a potentiallyconfounding material.

Electronic Noise.

An electronic system, that is designed to operate at temperatures aboveabsolute zero, is subject to generation of thermal noise in each of thesystem's components. The characteristics of additional types of noise(such as Shot noise or 1/f noise, for example) depend mostly on devicetype and/or manufacturing quality and semiconductor defects, such asconductance fluctuations. In measurement systems, the noise is an erroror undesired random disturbance of a useful information signal,introduced before or after the transmitter, detector, and decoder. Thenoise is a summation of unwanted or disturbing energy from natural, andsometimes, man-made sources. If noise is caused by a number ofuncorrelated sources, the central limit theorem (also known as the “Lawof large numbers”) states that the noise exhibits a statistically“normal” or Gaussian, distribution. This means that a number oftechniques (discussed below) can be used to extract very accuratemeasurements from below the noise floor. The quality of such extractionis dramatically improved if different causes of noise do not correlatewith one another or if they are correlated, they correlate in anidentifiable fashion.

Non-Linearity of the Electronics and Optics.

When non-linearities are present in the measurement system, and whendifferent measurements are carried out under different non-linearityconditions, the algorithm based on the matrix approach as discussedabove, which is based on linear transformations, may produce erroneousresults. Referring to FIG. 8A, in which an input signal level iscompared to 4 different system output signal levels, linear signal 810represents a linear output the value of which is two times its inputvalue. Similarly, linear output 815 shows a one-to-one relationship withthe input signal, and linear output 820 shows a one-half input levelsignal. The non-linear data series has a low but increasing gain region825, a high gain region 826, and a high output but very low incrementalgain region 827.

If the signal from the EMW at one wavelength encounters a stageoperating in a regime similar to 825, while the EMW at a secondwavelength is operating in a regime similar to 826, the relativedifferences in signal level at the second wavelength will be exaggeratedwhen compared to that at the first wavelength.

Linearization before combination is important for accurate results.Calibration and compensation circuits can be used as well as calibratedrange selection circuits to avoid the inaccuracies that may be caused bythe non-linearities of the electronics and the optics.

Ambient Radiation.

There is often ambient radiation illuminating the samples beingmeasured. Techniques to distinguish the effects of the ambient radiationfrom the measuring probe must be used. In addition, if the ambientradiation overloads the input of the receiver, this task becomes moredifficult.

SNR Improvement Techniques

The signal accuracy improvement techniques include the previouslydescribed improved choice of wavelength used for spectrometricmeasurements. The signal accuracy improvement techniques also includethe techniques described below, including, but not limited to pulsingthe transmitted radiation, providing variety in the on and offsequences, performing correlation between the transmitted radiation andthe received signal from the sample under measurement, repetition of themeasurement, using a variety of techniques to localize the portion ofthe sample being measured, using additional receivers, using thevariations in the impulse response produced by a variety of techniques,selecting particular portions of the impulse responses according toalgorithms, and using physical structures to modify the path from thetransmitter to the sample and from the sample to the receiver.

The design of many communications and measurement systems must take intoaccount the fact that the transmitted signal (in case of communicationssystems) or measured values (in case of many other systems) may bedegraded due to distance propagation, noise in the environment, and manyother factors. The “signal” that emerges at the output end of acommunication path or at the detector of a measurement system can,however, be improved (that is, the SNR of the detected data can beenhanced) with the use of various statistical processing approaches. Therelative strength enhancement of the received signal is referred to as“signal processing gain,” or simply “processing gain,” and cansignificantly improve the SNR of the end-to-end system.

In one embodiment, the sequence of “ones” and “zeros” of the excitationsequence is generated by an appropriate hardware component of the systemoperating according to an algorithm on the transmit-side of the systemand, as such, is defined by and/or in the system. On the receive-side ofthe system, the received signal can be converted into a stream ofelectrical “ones” and “zeros” and then passed through a “digitalcross-correlator,” which conceptually moves (i.e., a sliding windowrelative to time) the ones and zeros past a matching, pre-storedpattern. When the two signal patterns are delivered to thepost-processing unit and line up or are found to conform to one anotherin the post-processing unit, the “digital correlator” responds that ithas found a “match” between the transmitted pattern and its pre-storedpattern, and generates a “match signal.” The amplitude of the “matchsignal” is a measure of the “processing gain” from using the correlationapproach; in general, the longer the pattern of ones and zeros, thestronger will be the amplitude of the “match signal.” If the transmittedpattern is different from the pattern stored in the cross-correlator, no“match” will be reported by the system.

The advantage of this approach is that the received signal can beafforded to be substantially below the “noise floor” of the signal path.With reference to FIG. 8B, it can be seen that an extended excitationsequence or signal can be constructed from an edge-rich concatenation ofzero-values and two-values such that a sequence with transitions atvarying intervals is produced. This excitation sequence can be buried orhidden in noise. With reference to FIG. 8B, the amplitude of theexcitation sequence is attenuated such that the noise to signal ratio(NSR) reaches the level of 34:1 on the three graphs of the second rowfrom the top, 68:1 on the graphs of the third row from the top, and136:1 at the fourth row from the top. The four graphs in the middlecolumn of graphs of FIG. 8B illustrate the burst of the excitationsignal added with a zero-mean 4σ Gaussian noise source. The four graphsin the rightmost column show the result of performing a crosscorrelation of the excitation sequence signal with the noisy signal. Thepresence of the 10,000-bit excitation sequence signal is recovered fromthe noise even if the transmitted signal is buried in the backgroundnoise, at levels of 1/34th, 1/68th, or even 1/136th of the noise level.The longer the excitation sequence, the higher will be the potential SNRof the recovered signal; a properly chosen 2000-bit excitation sequencewill have a higher SNR than a 1000-bit sequence, and so on. Notably, thecross-correlator may optionally yield a measure of elapsed time as well;if there is a reference time “tick,” the correlator will indicate notonly how strong the recovered signal is, but when the match occurred intime. In the following paragraphs, we describe in more detail how thiscross-correlation approach can be employed. For the oximetry applicationdescribed, where the body-worn oximeter is to be tiny and powered by asmall coin cell, the cross-correlator will beneficially be implementedas custom-designed integrated circuit.

With reference to FIG. 8C, a series of graphs of a signal amplitude overtime are shown. Graph 885 represents the impulse response of a samplebeing measured by a spectrometric system. This impulse response is usedin the rest of the illustration in FIG. 8C. Graph 886 represents theidealized single chip input pulse to the system. The term “chip time”refers to the duration between opportunities for a discretetime-quantized excitation sequence to change from one level to another.The term “chip” in this context refers to the signal during the timebetween two successive change opportunities. The horizontal axis ofgraph 886 has tick marks at chip boundaries. Specifically, because thisillustration shows a spectrometric system characterized by chip durationof five time units, the tick marks are numbered at times equal to 0, 5,10, 15, 20, 25, and 30 and the period in between each tick mark isreferred to as a chip. Therefore, a single chip pulse 886 is a signalwhich has an amplitude at a low level, then at a chip boundary makes atransition to a high amplitude, and at the next chip boundary makes atransition to a low amplitude.

The graphs of FIG. 8C are provided for the purposes of illustrationonly; amplitudes are normalized where appropriate to a high level of oneand a low level of zero. Given the system impulse response 885, theresponse of the system to a single chip pulse input 886 is representedby the graph 887. This pulse response 887 can be used as a visualreference for the waveforms represented by the graphs below, 890, 894,and 895. An implementation of the invention produces a series of chipsin the form of an excitation sequence 888. Notably, the excitationsequence 888 as shown can also be represented by the sequence of binarydigits: {1, 0, 0, 0, 0, 0, 1, 1, 0, 1, 1, 0, 1, 0, 1, 0, 0, 0, 0, 1}.Graph 889 represents the response, or the output, of a system withimpulse response 885 and excitation sequence (input) 888. This signal isno longer a binary signal, but best thought of as an analog signal, or acontinuously varying signal. There are varieties of circuits andmathematical operations that can be used to take the excitation sequenceas a baseline, and find the correlation over a time range of thatbaseline with the output 889. One such method is to take the dot productof the two signals to provide a correlation value at time offset zero,and the dot product of the output 889 and the excitation sequence 888time delayed by Δt to find the value for the time Δt. Graph 890represents the result that is the correlation of the input of the systemwith the output of the system. The noiseless correlation 890 issubstantially identical to the pulse response 887. The graphs 887 and890 appear different because the single chip pulse used as input to 887starts at time 5 rather than time 0.

When noise is introduced into a measurement system, the noise becomespart of the output. When measuring the effect of the controlled input tothe system, the effect of the noise must be reduced. Several techniquesare illustrated in the present disclosure to allow a measurement to bemore accurate—even in the presence of noise. Graph 892 represents theexcitation sequence 888 on the local scale. Gaussian noise at +25 dBcompared to the signal level has been added to the output signal with anominally maximum envelope represented by the graph line 891 and anominally minimum envelope represented by the graph line 893. Graph 894represents the result of the correlation operation on the noisy outputbased on an input excitation sequence of length 2¹⁴-1, while graph 895represents result of the correlation operation on the noisy output basedon an input excitation sequence of length 2¹⁸-1.

With reference to FIG. 8D, graph 885 from FIG. 8C is shown in thecontext of an explanation of the waveform developed over time as theimpulse response of a sample configured as shown. The lower part of thediagram includes a sample volume 830 juxtaposed with a transmitter 832and a receiver 835, located a distance away. When an energy impulse isinjected into the sample at time zero, there is, at time zero, noresulting output detectable at the receiver. Radiation which travelsballistically from the transmitter 832 to the receiver 835 will travelthe ballistic or straight-line path 841. Until that portion of the inputenergy arrives at the output, as is indicated on the first interval ofthe graph 840, there is zero received amplitude. At time point 841, EMradiation traversing the sample along the shortest path 841 activatesthe receiver and causes the spike in received energy (indicated by aspike of graph 885 at about 4 units measured along the abscissa). Someradiation follows paths through the sample close to the surface,following paths that are short because they remain substantially inregion 842. This radiation is responsible for the rise over time in thegraph region 842. Depending on the characteristics of the samplematerial, with regard to diffusion, absorption, reflection, andanisotropic features of these effects, more or less of the radiationwill follow short paths compared to long paths. The longer pathspenetrate farther into the sample and take more time for the EMradiation to reach the receiver; and statistically, over time, much ofthe remaining radiation is not detectable by the receiver as it eitheris absorbed, or exits the sample undetected by a receiver. Thesuccessive depth of regions shown by 843, 844, 845, and the bulk of thesample 846 illustrate the probable penetration of the energy that isreceived in the time of graph 885 corresponding to the regions. It wouldbe misleading to say that the correspondence is exact; it is not.However, different delays to correspond to different regimes of paths,and if the regime can be characterized, the results of thesecharacterizations can be used to inform later measurements. Embodimentsof the invention include systems that advantageously use thefine-grained ability to segregate the measurements of the delay cohortsof paths, and to perform differential characterization of these delaycohorts.

While the general use of the term ‘impulse response’ is the response ofa system to a unit sized burst of energy injected in zero time, theembodiments herein can more accurately be described by their response toa chip-long pulse. This chip-long rectangular-shaped pulse may be longerthan the duration of detectable levels of a typically conceptualizedimpulse response; in that case, the pulse is similar to a time-limitedstep function, or more simply, a pulse. Alternately, the chip-time maybe considerably shorter than the detectable duration of an impulseresponse and the chip pulse then more closely resembles an impulse. Oneskilled in the art will recognize that the shape and duration of theinput pulse is convolved in time with the actual impulse response andthat what is received and quantified is the pulse response. Depending onthe various noise reduction techniques embodied in the various systemsdescribed, the pulse response, which can be measured with significantlyimproved accuracy and resolution, can be usefully de-convolved toprovide even finer-grained information. Although the use herein does notstrictly fit mathematical interpretation of the term, when used herein,‘impulse response’ encompasses the extracted response, of the systembeing measured, to a single chip-wide pulse of input. The data collectedis quantized and limited in time; however, the term is used herein as itcaptures the combination of impulse, ‘chip’, step, and pulse response.Impulse response is additionally used describing the calculated resultof combining received impulse responses from multiple λ_(i)'s.

As used herein, concentration is an amount of a particular materialcomponent or analyte. This may be a percentage, or when volume, forinstance, is determined, an absolute value may be calculated.

An additional reason for employing lasers (VCSELs) rather than LEDs isthat the VCSELs can be pulsed much faster than the LEDs, which is ofvalue when long excitation sequences are to be transmitted, and providesfor a more accurate measure of path length.

Examples of Systems for Measuring Properties of Material Components

Referring to FIG. 9A, a system to measure one or more properties of oneor more material components includes a transmitter 900 that forms aprobe 915 (for example, in the form of light at one or morewavelengths). The probe is directed to interact with material component920, as a result of which the properties or characteristics of the probeare changed to reflect such interaction and to form an altered probe;the specific nature of the change depends on the properties of theanalyte. The presence of the altered probe 925 is detected by thereceiver 930. The data representing the altered probe 925 and acquiredwith the receiver 930 is quantified by the data quantifier 935, in waysdescribed below, and the quantified data is then processed by thecomputational unit 940 to produce an ascertained property 945.

With reference to FIG. 9B, the system may be augmented to increase thesignal-to-noise ratio from measurements obtained using a transmitter 900and receiver 930 by having the transmitter driven with an excitationsequence generator. In this case, each of the components of light (ofthe probe) at different wavelengths is modulated by a corresponding,unique excitation sequence. In this embodiment, the receiver 930 feedsthe signal through a data quantifier 935 to a computational unit 940composed of at least the parts of a correlation circuit 955 and anaggregating circuit 960. In this embodiment, each correlator has asinput to its operation both the received signal and the unaltereddefined excitation pulse sequence 950 through a direct connection 965.The correlator can beneficially reference the excitation sequence toperform the cross correlation of the received signal with the probe.Examples of the implementation of such correlation procedures aredescribed below. With reference to FIG. 9C, an additional system isshown, adapted to transmit and receive multiple excitation sequences tofacilitate correlation thereby increasing signal-to-noise ratio and/ordecreasing motion artifacts. The transmitted signal is directed tointeract with the material component that may be a turbid medium or amolecular dye 970. The probe includes light at one or more wavelengths,light portions at each wavelength modulated by a corresponding uniqueexcitation sequence. Multiple correlators can be used to provide forconcurrent processing of the multiple wavelengths of the altered probe975.

Controlling the Location of Sampling Point(s) Using Micro-Mirrors

In making measurements with light sources, it is often desirable tocontrol within tight tolerances the location where the light makescontact with the sample to be measured. One application that makes useof such measurements is the read/write mechanism on an optical diskdrive such as digital video disc (DVD) or blu-ray disc (BD). The currentstate of the art positions the laser using a worm drive with servocontrol. The advantage of such a system is the position of the laser iscontinuously variable with accuracy limited by the screw, the motor, andcontrol system. However, the disadvantage is that it requires asignificant amount of hardware, time delay, and energy (to move themass) to accurately position the laser with respect to the samplelocation.

Another example application is Universal Product Code (UPC) scanners. Inone implementation, a laser light is scanned across the UPC using asingle rotating mirror. In another implementation, a laser light isscanned across the UPC using a hologram. Unlike the invention describedherein, the incident angle of the light to the UPC is not controlled,nor is the distance from the light source to the sample or from thesample to the receiver.

The advantage of the invention described herein versus a worm drive isthat it requires less hardware, little time, and lower energy toposition the laser beam at the sample location to be measured. Relatedbenefits are increased reliability due to simplified mechanical partsand the reduced size of the apparatus. Unlike UPC scanners, theinvention described herein simultaneously controls the distance and theincident angle to the sample. Though controlling the distance andincident angle is not necessary for UPC scanners, it is necessary forother sampling applications.

With reference to FIG. 10A, the use of micro-mirrors 1080, possibly ascontained in a micro-electro-mechanical system (MEMS) but not limited tosuch, each with controllable manipulation between at least three (tiltedleft 1081, non-tilted 1082, and tilted right 1083) positions allows theoptical path and thereby the interaction of the probe 1015 with thematerial component 1020 to be altered during operation. Though themirrors are shown in an array pattern in FIG. 10A, this is onlyexemplary, as other patterns are anticipated. The transmitter 1000 maybe any light source such as a laser, VCSEL, or LED, but is not limitedto those. The receiver 1030 may be a photodiode, PIN diode (a specifictype of photodiode), photomultiplier tube, or charge coupled device(CCD)), but is not limited to those. The transmitter emits a waveparallel to the face of the micro-mirror array 1080. One or moremicro-mirrors at precise locations in the micro-mirror array 1080 aremoved from a non-tilted 1082 position to a tilted 1081 position thatdirects the probe 1015 toward the material component 1020 to bemeasured. In one embodiment, the altered probe 1025 is directed to thereceiver 1030 by at least one micro-mirror tilted 1081 in one directionand at least one micro-mirror tilted 1083 in the opposite direction.

FIG. 10B illustrates an extension to the embodiment illustrated in FIG.10A. In this embodiment, there are three micro-mirror arrays: 1080 a.1080 b, and 1080 c. Micro-mirror arrays 1080 b and 1080 c are 1×Nmirrors in dimension. In addition, the two arrays are perpendicular tothe plane of micro-mirror array 1080 a. The transmitter 1000 transmitslight in parallel to micro-mirror array 1080 b where it is reflected offtilted mirror 1081 b (1 of N possible micro-mirrors in the array)directing the light across one row of micro-mirror array 1080 a andparallel to micro-mirror array 1080 a. At that point, the apparatusworks as shown in FIG. 10A, with the light interacting with the materialcomponent and being reflected off of micro-mirror 1083. Unlike FIG. 10A,the light is then reflected off micro-mirror 1081 c in the micro-mirrorarray 1080 c. The light reflected by micro-mirror 1081 c is thenreceived by receiver 1030. With this invention, each row on themicro-mirror array 1080 a may be used to sample a two dimensional areaof the material component. Micro-mirror array 1080 a is shown with eachrow offset from the previous one (i.e., staggered rows). The staggeringof the rows is only exemplary. The rows could also be one above theother and therefore not staggered.

FIG. 10C shows the efficacy of this invention for precisely controllingthe distance from a point 1040 where the light is reflected from thematerial component 1020 to the receiver 1030. For the example shown, twomirrors are used to control the distance and the array isone-dimensional (i.e., there is only one row of mirrors). The firstmirror is tilted as shown by 1081 and the second mirror is tilted asshown by 1083. The distance (i.e., d) from the point 1040 to thereceiver 1030 is shown by equation 45. Using a simplifying assumptionthat the material component 1020 and the micro-mirror array areparallel, then when the mirrors to the right of the two tilted mirrorsare tilted, and the prior two tilted mirrors are put in non-tiltedpositions, then the distance is that shown in equation 46. Thus, thedistance is decreased by the distance between two adjacent mirrors. Theequation is generalized to equation 47 where N is the number of mirrorsto the right of the positions shown in the figure. Note that thissimplified equation is dependent on the material component beingparallel to the micro-mirror array and the angle α (and therefore θ)being constant. One of the advantages of micro-mirrors is that the angleof tilt is repeatable and consistent from mirror to mirror.d=x ₂ −x ₁ +h  (45)d=x ₂ −x ₁ +h−x ₃  (46)d=x ₂ −x ₁ +h−Nx ₃  (47)

The advantage of having an array of precisely spaced micro-mirrors isthe capability of selecting with fine resolution the distance to thereceiver from the “point” where the probe is scattered by the materialcomponent. By optimizing this distance, the desired point to be measuredcan be located and the resulting measurements improved. In addition, thematerial component can be scanned (by changing which one or moremicro-mirrors is tilted) to quickly locate areas of the materialcomponent with the properties desired.

Other advantages include compensation for variability in the y direction(as shown in FIG. 10C). If y increases then Nx₃ can be increased to keepthe total distance “d” constant.

The paragraphs above have made use of a simplifying assumption, thatthere is a single point where the probe 1015 reflects from the materialcomponent 1020, producing the altered probe 1025. However, many materialcomponents of interest do not simply reflect the light, but in addition,absorb, refract, and scatter the light. Therefore, there is not a singleray coming from the light interacting with the material component but amultitude of rays scattered by the material component as it interactswith the light. This understanding leads to other useful alternativesfor operating this device. FIG. 10D shows a single tilted mirror 1081with receiver 1030 receiving light reflected from one or more of thenon-tilted mirrors 1082. FIG. 10E shows the complementary operation withtransmitted light reflecting off one or more non-tilted mirrors 1082,interacting with the material component, and being scattered back to atilted mirror 1083 and reflected to the receiver 1030.

Some material components will diffuse the light. Some of this light willexit the sample at parts of the micro-mirror array that are not intendedto be active. To reduce the amount of diffused light reaching thereceiver though non-tissue paths a barrier can be introduced. FIG. 10Fshows how such a barrier 1040 may be used. Even if the barrier is incontact with the material component, the light may still travel to thereceiver as it travels through the material component and exits on theright hand side of the barrier. An alternative barrier technique isshown in FIG. 10G. In this case, a column of micro-mirrors is used as abarrier. By tilting one or more micro-mirrors 1083 in the oppositedirection to those micro-mirrors 1081 reflecting the light, then lesslight is allowed to pass from the left of the micro-mirrors 1083 to thereceiver 1030. However, as previously shown, the 1083 micro-mirrorsreflect the light to the right toward the receiver.

FIG. 10G shows the second dimension of operation of a two dimensionalmicro-mirror array. The tilted micro-mirrors 1081 are used to selectportions of the transmitted light (probe 1015). In this manner, withmirrors smaller in dimensions than the width of the transmitted wave oflight, a wave smaller than that generated by the transmitter isreflected for a probe. The micro-mirrors 1083 tilted in the oppositedirection are used to trap light that passes by the micro-micro-mirrors1081. Another possible use of the light trapping mechanism is toscavenge energy from the light that will not otherwise serve a usefulfunction.

Finer granularity of distances than those between micro-mirrors can beachieved by using staggered rows of micro-mirrors, such as shown in FIG.10H. In this case, the stagger distance is less than the distance fromthe centerline of one mirror in a column to the next mirror in the samecolumn. FIG. 10H also introduces multiple transmitters 1000 and multiplereceivers 1030. One use of multiple transmitters is to allow multiplewavelengths to be transmitted, and if desired simultaneously. One of theadvantages of this invention is that each wavelength can, if desired,have a unique distance to the receiver. Multiple receivers may be usedto facilitate the selection of multiple distances from the sample to thereceiver. Finally, each receiver can also be dedicated to a particularwavelength, or range of wavelengths, of light. This can be useful to aidin separating signals and in selecting more than one type of receiver,as in ones particularly suited to the wavelengths being transmitted.

FIG. 10I shows a micro-mirror pattern using mirrors with a hexagonalshape. This pattern helps maximize the amount of reflected light byminimizing the amount of light passing between the tilted 1081 mirrors.

FIG. 10J shows a micro-mirror array using a diagonal pattern.Micro-mirrors in production today often have the pivot axis across thediagonal.

FIG. 10K illustrates the application of the invention for use in readingor writing an optical disk 1090 (also showing the hole 1092 for thespindle). The tracks 1091 for this example are aligned with a single rowof micro-mirrors. Using the micro-mirror arrays illustrated in FIG. 10Ballows for all rows in micro-mirror array 1080 a to be used as necessaryto align the light rays with the tracks on the optical disk 1090. Most,if not all, optical disks currently manufactured use spiral tracks.Several techniques may be used with this invention to facilitate thetracking of spiral tracks. One such technique is a worm drive with servocontrol, as typically used in optical disk drives. However, in thisapplication much less range of motion is required than in typicaloptical disk drive laser assemblies. Another technique is to enabletracking through the micro-mirrors by adjusting the angle of tilt forthe mirrors both on the transmit and on the receive side. An additionaladvantage of this invention is with known track-to-track distances, oncea track is located, then the mirror positions required to find any othertrack is known and can be changed rapidly.

In reference to all FIG. 10A-K, each describes an implementation wherethe angle of the tilted mirrors is constant. There are obviousadvantages of allowing multiple angles and even continuously variableangles. For instance, continuously variable angles may be used forsampling continuous locations rather than the discrete locations formirrors with a single angle.

The use of micro-mirrors in the path of the probe or altered probe mayreduce the received signal. This technique can be made more useful inpractice by utilizing the various encoding techniques andpost-processing techniques described herein.

Sequence Generator and Excitation Sequences

Linear Feedback Shift Register.

An excitation sequence generator may be implemented as a programmablelinear feedback shift register (LFSR); different settings on theinternal registers of the LFSR will produce different excitationsequences. There are several ways of designing the LFSR to generatedifferent patterns, which patterns may be quite long without repetition;and the transistor-level implementation is not complex and will not be avery high-transistor-count circuit. In one embodiment, with reference toFIG. 15A, a shift register 1500 produces a shifted output 1510. Thisoutput 1510 is “fed back” to the syndrome logic 1515. The syndromelogic, based on the configuration register 1520 provides the output 1510to configured elements of the shift register, wherein the output isexclusive-ored with the bit being shifted in to modify that bit value sothat it may be stored in the next element of the shift register whentriggered by the clock pulse 1505. Optionally, the configurationregister 1520 can be written with a new value, said value being suppliedfrom a value source, e.g. but not limited to a processor (not shown)data bus. Toggling the write control signal 1530 will cause theconfiguration register 1520 to update to the value supplied in the valuebus 1525.

Incorporated in this LFSR design is a “start pulse” output signal line1540, which may be optionally used in the receiver and correlatorcircuitry to determine the sequence temporal boundaries. Thecommunication of sequence boundaries can optionally be used forautomatic gain control in the receiver as described below. It isrecognized that the sequence start signal 1540 can be produced in anumber of ways. In the pictured embodiment, a pattern matching circuit1535 is used.

With reference to FIG. 15B, the details of an exemplary LFSR are shownwith the shift register elements 1555 shown as digital latches,receiving an input and propagating that input to a storage element inthe latch and providing the contents of the storage location at theoutput. The input is propagated to the internal storage location on therising edge of the digital signal designated clock 1505. The leftmostelement in the shift register takes its input from the leftmost elementof the configuration register 1565. The logic is arranged in groupscorresponding to the logic needed to update each element of the shiftregister. The collection of one each of a configuration register element1565, an AND logic gate 1570, an exclusive-or (XOR) logic gate 1545, anda shift register storage latch 1555 comprises the syndrome logic 1515used to update the shift register each time the input clock 1505toggles. Each successive (traversing from left to right) shift registerelement takes its input from the output of an XOR logic gate 1545. Thetwo inputs of the XOR gate 1545 are the previous element of the shiftregister 1555 and the output of an AND logic gate 1570. The inputs tothis AND gate 1570 are the value of the element of the configurationregister corresponding to this bit group and the output of the shiftregister 1510. It is appreciated that this logic, when configured withan appropriate configuration value will provide a deterministic sequenceof logic ones and zeros in a repeating pattern. With reference to FIG.15B, the center of the figure contains four ellipsis patterns, i.e., ‘ .. . ’ these patterns are used by those familiar with the convention toindicate the presence of an indefinite number of repetitions of theelements on either side of the ellipsis patterns. This convention ofillustrating a larger structure with a smaller number of repeatedpatterns is used in several of the illustrations including at leastFIGS. 15B, 17B, 17E, 18B, and 19.

Noise Based Generation.

With reference to FIG. 16, additional embodiments of excitation sequencegeneration are illustrated. FIG. 16A illustrates the use of a noisegenerator 1600, such as but not limited to a microphone, a back biaseddiode junction, a receiver of atmospheric noise, or other time varyingprocess. The output of the noise generator 1600 is fed to a low passfilter 1610. The low pass filter provides an average value reference asone input of a comparator circuit 1620. The other input of thecomparator circuit 1620 comes directly from the noise generator 1600. Aswill be appreciated, the instantaneous value of the noise generator maybe at times above the longer-term average provided by the low passfilter 1610, at which time the output 1630 of the comparator 1620 willbe driven to the extreme of a logic one. At other times, theinstantaneous value of the noise generator will less than thelonger-term average, and the output of the comparator 1620 will bedriven toward the other extreme, a logic zero. Using this circuit ananalog time varying signal can be converted to a binary signal.

Storage Based Generation.

A technique that may be used for the production of a defined excitationsequence is a predetermination of a sequence and the storage of theelements of that sequence in an addressable storage component 1640, asillustrated in FIG. 16B. The excitation sequence can be produced bydriving the address lines of the memory device with a counter 1650.Successive memory locations will be presented to the output of thememory, creating the sequence of values comprising the excitationsequence 1630. As a variation, when more than one sequence is needed, asin the embodiments which provide for the concurrent driving of multiplelasers, multiple bit wide memories can be used, with each bit indexbeing used to drive a separate laser. When calibration sequences differfrom measurement sequences, having access to multiple sequences canpresent advantages. In this case, the counters may be used concatenatedwith configuration registers to cycle through one group of memoryaddresses, then with a change to the configuration value, to cyclethrough another group of memory locations.

With reference to FIG. 16C, a processor 1690, connected through a memoryinterface unit 1670 to a memory unit 1660, fetching instructions fromthe memory and decoding 1680 and executing those instructions mayexecute a stored program or algorithm, and by writing to a memorylocation which is able to output a logic signal 1630, produce theexcitation sequence under program control.

The excitation sequence repeats in time over and over again, withoptional brief durations between each sequence when no bit transitionsare transmitted. Another way to generate this quiet time is to include afixed set of zeros at the end of the defined excitation sequence. Theoutput from the excitation sequence generator is fed to the laser drivercircuitry, which in turn modulates the excitation sequence onto thelaser. The laser driver and the light source itself is an analog portionof this embodiment of the measurement device. Either an on-off keyingmodulation can be used in one embodiment, or in another embodiment thelaser light is modulated from bright to dim and back, but does notcompletely extinguish as long as each instance of the excitationsequence is being transmitted.

Examples of Correlators Useful for Spectrometry

Binary Digital Correlator for Pulse Oximetry.

With reference to FIG. 17A, techniques for producing a digitalcorrelator for use in spectrometry are described. These techniques canalso be used in the additional applications described herein.Correlation embodiments that improve SNR significantly mayadvantageously be implemented in the form of an integrated circuit,thereby allowing the correlator to function as part of a mobile,body-worn unit. In the case of pulse oximetry, which is a relativemeasure of blood oxygen saturation, not an absolute measure, it is notrequired to measure the path length of the light through the tissue.

One embodiment of a correlator for spectrometry contains a preponderanceof digital circuitry using a binary received signal to detectcorrelation.

On the receive side of the circuitry, the laser light, after passingthrough the tissue, impinges on a solid-state photodetector 930 such asa P-I-N (PIN) diode, an avalanche photodiode, or a photomultiplier tube,any of which will convert the optical signal into an analog,time-varying electrical signal. This signal is then fed into an analogamplifier stage 1710, whose amplification factor is modified in aclosed-loop fashion as follows: the start pulse 1540 from the excitationsequence generator informs the controller 1700 when the pulse train isto begin; and if by the end of the excitation sequence, as indicated bythe next start pulse 1540, the correlator 1730 has not generated acorrelation “match” 1736, then for the next receipt of the pulse train,the amplification factor will be increased. This increase will continuein repetitive steps until either the correlator shows a correlationmatch, or the amplification factor is at its maximum. A thresholddetector 1720 on the output of the automatic amplification stage 1710 isneeded to convert the analog electronic signals from the photodetectorinto a binary pulse train, but again this is well-known circuitry.

The correlator digital circuitry 1720 is clocked with the same clockpulse train 1505 that drives the excitation sequence generator.

With reference to FIG. 17B, the elements of the digital correlatorconsist of latches 1750 (grouped in registers) to hold thepre-determined bit pattern of the excitation sequence; a linear train offlip-flops 1740 (latches) through which the pulse train is propagated,clock cycle by clock cycle 1505; a group of XNOR gates 1745 which detectif the propagating bits in the received pulse train match the ones orzeros in the static excitation sequence holding register (one-bits mustmatch one-bits, and zero-bits must match zero-bits, across the entireexcitation sequence, for a full match to be identified). The output ofeach of the XNOR gates is fed to a binary adder 1765, which generates asum value on each clock pulse. When a full match occurs, the binaryadder will generate a large binary output value, which defines that amatch occurred, and when in time the match occurred. On average, onlyhalf of the bits will match the pattern during each clock cycle.However, during the clock cycle corresponding to the alignment, that is,when a match occurs across all excitation sequence bits simultaneously,the output of the binary adder 1765 will be at a maximum value, theexpected magnitude of which is well known and scaled by sample losses.

As used herein, an “XNOR gate” is a two-input one-output circuit withthe characteristic that the output will be a one when the inputs match;otherwise, it will be a zero. As used herein, the term “binary adder” isa circuit with N-inputs and ┌log 2(N+1)┐ outputs such that the outputstaken together represent a binary number that is the count of the inputsthat are one. An example of the binary adder is a circuit with 10inputs, I₀, through I₉, and 4 (┌log 2(10+1)┐) outputs, O₀, O₁, O₂, andO₃. When the inputs, I₀ . . . , I₉ are, for example {0, 1, 1, 1, 0, 0,0, 0, 0, 0}, the outputs, O₀ . . . , O₃ would be {1, 1, 0, 0}representing the value of three. The same output would be produced withinput patterns of, for example {0, 0, 0, 1, 1, 1, 0, 0, 0, 0}, {0, 0, 0,0, 0, 0, 1, 1, 1}, or {0, 1, 0, 1, 0, 0, 0, 1, 0, 0}.

With reference to FIG. 17A, the input to the correlator is a comparisonbetween the output of the amplifier 1710 and a gain reference voltage1715. The comparator 1720 compares the two values and, similar to thenoise resolving circuit of FIG. 16A, produces a binary output 1735 tofeed into the correlator.

With reference to FIG. 17B, the register elements 1750 that hold thereference of the excitation sequence to be compared against, may beupdated with new values. The new values are presented by the N-bit databus 1760, and will be stored in the excitation sequence register bits1750 when the write pulse 1755 toggles. In this manner, a circuit may beused with multiple reference sequences.

A carefully designed digital correlator, fed from an excitation sequencegenerator, can detect a laser-generated optical signal deeply buried innoise and degraded due to absorption and scatter of the laser light. Thereceive-and-process system has a front-end gain control stage whose gainis controlled in steps until the correlator sees the excitation sequencebit stream or until the AGC amplifier reaches an upper limit. Thecontroller 1700 in FIG. 17A will follow the process of FIG. 17C todetermine the range over which the signal can be detected. Without lossof generality, but rather to facilitate understanding of the algorithmused, the following variables are used: G is the gain of the amplifieror the voltage of the reference signal used as input to the comparator1720. C is the count of matching correlation elements. It is the outputof the binary adder 1735. C_(max) is the largest value of C observedbetween two start pulses; therefore, it is the level of highestcorrelation during one excitation sequence. T is the threshold value ofC that is chosen to represent a high likelihood that the transmittedsignal has been detected. G_(min) is the smallest value of G for which avalue of C_(max)>T has been observed. G_(max) is the largest value of Gfor which a value of C_(max)>T has been observed. V is the numericalvalue of the measured output and is produced by a function ƒ, forexample but not limited to the average, of the two measured values,G_(min) and G_(max). Following this notation, the algorithm starts 1770initializing G 1772 to its lowest value and proceeding to initializingG_(min), G_(max), and C_(max) 1774 to maximum, minimum, and minimumvalues respectively. Once initialized, the algorithm waits for a startpulse 1776. After a start pulse has occurred, the algorithm observes thevalue of C, recording the maximum value of the value of C observedbefore the next start pulse 1780 is observed.

After the start pulse 1780 is observed, C_(max) is compared to thethreshold T 1782. If C_(max)>T, then G_(min) and G_(max) are updatedwith G_(min):=min(G_(min), G) and G_(max):=max(G_(max), G) 1784. G isthen increased by a step 1786, and then tested against the last valuefor G. 1788. If G is not greater than the last value of G, the processcontinues from 1774. If the last value of G has been evaluated, then thealgorithm calculates V as a function of G_(min) and G_(max) 1790, andends 1792.

The value V of the detected signal can be used to determine the value ofthe parameters to be measured, by, for instance but not limited to, amatrix multiplication of the measured values and a reference matrix. Thematrix multiplication can be performed by means well understood using anelectronic processor. Additionally, with reference to FIG. 19, a set ofmeasured values can be fed as binary-coded input 1900 values into anarray of multiplier circuits 1920, with outputs 1930, which feed toadders 1940, one for each output 1950 corresponding to a calculatedproperty. Using structures of this form or similar compositions, a fixedcircuit can be composed in a non-programmable structure to calculate thenumerical value of a property to be measured. The removal of separatenon-linearities can also be performed, by at least one of table lookup,range rescaling, limit alerts, and gain adjustment techniques asdescribed related to FIG. 17A at element 1900. This conversion elementis usefully included when the range of use of the inputs exercisenon-linearities in the measurement path.

In order to quantify the results, the attenuation of the received signalshould be initially measured. Such measurement includes transmitting anexcitation sequence, receiving the attenuated sequence, and correlatingthe received values in time. A device, such as that shown in FIG. 11, isused for such a purpose. Once the attenuation of the received signal isderived, then equations 37c and 37d, for instance, can be used tocalculate the concentration of a property in the material component. Theapparatus shown in FIG. 19 is an appropriate circuit for performing suchmathematical operations assuming the measure times the weightingcoefficient and the summation are stateful operations (i.e., values fromprior time slices are available to perform calculations). The “weights”in these calculations refer to wave properties such as the optical paths(or traversed distances) and extinction coefficients. Other fixedcircuits, such as systolic arrays would also be appropriate to performsuch calculations. Depending on real-time constraints and speed ofprocessing, programmable processors such as those used fordigital-signal-processing (DSP) or even general-purpose processors mayalso perform the necessary calculations.

In addition to equations 37c and 37d, it is understood that theintroduction of bias and noise (see equation 8) to the calculations istypically required to attain quantifiable (versus relative orqualitative) results. Calibration operations, for example, are one meansof deriving the values to be used for bias and noise.

The method to correlate a received signal with a transmitted signalprovides for determining the delay in time between the two and therebythe delay from transmitting to receiving is a known quantity.Additionally, connections from the transmitters to the receiver, such asthe connection 965 in FIGS. 9B and 9C, may be used for synchronization.As a result, the delta in time between signals received from two or moretransmitters is also a known quantity and the signals can thereby becorrelated in time.

N-Bit Digital Correlator.

While very applicable to a pulse oximeter (and providing for compact andlow power implementation, and thus so advantageous in use), theembodiment described above with a binary comparison has given up much ofthe processing SNR gain made possible by the use of long excitationsequences. This is due to the single point of discrimination betweenones and zeros at the comparator 1720. For many spectrometryapplications, providing additional signal discrimination is useful. Withreference to alternative FIG. 17D, similar to FIG. 17A with twosignificant changes, an N-bit digital correlator is shown. Thecomparator 1720 has been replaced with an N-bit analog to digitalconverter (ADC) 1721. The amplifier now optionally takes an input forgain from the gain reference block 1716.

With reference to FIG. 17E, the N-bit wide output of the ADC 1721 isused as input 1735 to the first 2 times N-bit wide delay line register1742. The input is fed to N of the bits of the register 1742; the input1735 is also fed to a negation block 1736. The output of the negationblock 1736 feeds the other N bits of the delay line register. Based onthe clock pulses 1505, the values in one delay line register 1742 feedthe next register such that a value from one point in time will passfrom one register to the next and on to the end with a step made duringeach clock pulse.

Responsive to the excitation sequence register bits 1750, N bits of theoutput of the delay line registers are chosen by each multiplexor 1747.The excitation sequence register bit set to one will select thetime-delayed value of the ADC, a zero bit in the excitation sequenceregister bit will select the negated time-delayed value. In any giventime cycle, the outputs of all of the multiplexers will be addedtogether at 1767 to create the output for the controller. If theexcitation sequence stored is M bits long, then the sum of the outputs1747 will be log₂(M)+N bits wide if no bits from the least significantend of the numeric bus are dropped.

It is appreciated that, without loss of generality, the negation block1736 can be placed in the data path after each delay line register 1742,and the delay line register can then have a collapsed width of N-bitsrather than 2-N bits. Additionally other rearrangements of the topologyof operations that accomplish the same goals can be applied. Thisprinciple is demonstrated below in the analog embodiments and in thehybrid embodiment.

Analog Delay Line Correlator-A.

If it is desired to carry out quantitative spectrometry, a measure ofthe optical path length of the radiation in the tissue is needed. Since,in, for example, making measurements of samples of the human body thepath length can be as short as 2-3 cm, and since light will travel in aturbid medium at perhaps 0.2-0.3 times the speed of light in a vacuum(or on the order of 200-500 ps/inch), very fine time resolution, in therange of 5-10 ps, is needed to resolve the optical path length withsufficient accuracy to be useful (see the discussion elsewhere in thisdocument as to why a measure of optical path length yields anabsolute-reading oximetry value rather than a relative value). Severalembodiments of a receive-and-correlation processing chain have beenidentified, as described herein.

With reference to FIG. 18A and FIG. 18B, the embodiments are based onthe concept of a traveling-wave amplifier, a concept that has been knownand practiced at least since the 1950s for various high-poweramplifiers, e.g., satellite communications systems, but has not beenapplied in either the context (spectrometry) nor in the structure(varying the inputs to the various sub-elements in two distinct ways) asdescribed herein. The design of the traveling wave amplifier is hereinaltered considerably, to such an extent that the design is only “akin”to a traveling wave amplifier, but is not actually one. In a classicaltraveling wave amplifier, the analog input signal is propagated, or“travels,” along a transmission line that is “tapped” at regularintervals, with the “tap” feeding the inputs of a successive chain ofdiscrete amplifiers. As the analog signal propagates along the “feed”transmission line, the voltage variations of the input wave are used todrive the input of each successive amplifier, which in turn injects a“slug” of current into a second “output” transmission line. Thus, thesignal in the output transmission line is amplified by the individualamplifiers through the superposition principle, such that at the end ofthe output transmission line a strongly amplified copy of the inputsignal emerges. This is the explanation of a conventional traveling waveamplifier. Regarding the present invention, the design changes to thetraditional amplifier are described below, which novel modificationsmake the concept work as an “analog correlator,” in contradistinction tothe “digital correlator” described previously. The need for this designis driven by a requirement for a high level of timing accuracy withrespect to the occurrence of the excitation sequence pattern at thereceive chain, in comparison to when the excitation sequence wastransmitted, that will yield a measure of optical path length in thetissue.

A correlator based on analog signal processing can be developed as amodification to, and extension of, a traveling wave amplifier. Withreference to FIG. 18A, the optical pulse train is converted by anamplifier 1810 to a time-varying analog signal 1812 representing theoptical pulse train, which is then propagated down an input transmissionline 1825, which in turn is “tapped” at regular intervals, with thetime-varying voltages at each “tap” fed into one input of an analogamplifier 1830. The analog amplifier in combination with the excitationsequence register bit 1750 and the multiplier 1830 implements the samecircuit function as in the digital correlator, i.e., XNOR 1745 ormultiplexor 1750. The time-varying output current from each analogmultiplier 1855 is fed into a current summing junction 1860, and theoutput of the current summing junction 1840 will only have a maximumnonzero value if the stored and propagating excitation sequence codesline up one-for-one at some point as the excitation sequence patternpropagates down the input transmission line. In the embodiment shown inFIG. 18A, the controller 1800 provides input to the gain control 1815 ofthe amplifier 1810. The second input to the multiplier is the excitationsequence register bit 1750.

With reference to FIG. 18B, the incoming analog excitation sequencepulse train amplified 1810 and allowed to propagate down the inputtransmission line 1812 (as described above). As each pulse in theexcitation sequence passes an amplifier 1850 and multiplier 1830, thatamplifier generates a slug of current that is fed by the multiplieroutput into a current summing junction, optionally through a resistor1855, or into one of a set of current summing junctions. This choice isa practical question of the particular technology used forimplementation. If the summing junction has a large dynamic range, thenfewer summing junctions will be needed and fewer cascade steps will berequired. The summing junction sums all of the current slugs from everyamplifier along the chain. Each amplifier/multiplier along the inputtransmission line in effect replicates in analog form the XNOR gatestructure described for the digital correlator, such that matching“ones” and “zeros” between the stored pattern bits and the propagatingexcitation sequence bits statistically line up, zero-to-zero, andone-to-one. It is appreciated that the amplifier and multiplier showncan be combined in a variety of technologies using techniques known inthe field.

Conversely, a one in the stored pattern aligned with a zero in thepropagating wavefront, or a zero in the stored pattern aligned with aone in the propagating wavefront, will produce uncorrelated outputamplitude from a given amplifier, contributing positive or negativecontributions to the summing junction. On an instantaneous basis, theoutput from the summing junction will be higher or lower depending onhow many zero-bits line up between the stored reference and thepropagating waveform, and how many one-bits line up in the same manner.Since this is an analog circuit, the maximum output from the summingjunction will only occur when the propagating code bits line up with thestored code bits, at which time the output from the summing junctionwill be maximum.

Optical path length and hence propagation time is very important in theabsolute reading oximeter, since optical path length enables thequantitative calculation of oxygen saturation. To achieve the timeresolution needed for path lengths of only a few cm, time resolutions inthe 10-20 ps range need to be achieved. In one embodiment, to accomplishthis resolution, we use a reference length of electrical transmissionline on the oximeter circuit board (such as, but not limited to, alength of stripline or microstrip), which contains both the lasers atone end, and the photodetector and analysis circuitry on the other end,designed with a known propagation delay; this is the “gold standard”delay. We then compare the gold standard delay with the measured delayfrom decoding the excitation sequence at the receiver end.

In the first embodiment of this reference approach, an electrical “startpulse” 1540 from FIG. 15A is generated by the same circuitry that drivesthe laser, and propagates down the reference delay line to theprocessing circuitry; the rising edge of the start pulse at theprocessing circuitry is used as a time reference against which theoutput of the analog correlator circuitry is compared. Because thepropagation delay along the known length of the reference delay line ismade short by controlling the material from which the circuit board ismanufactured, the rising edge of the start pulse creates a “time zero”point in the analog correlator circuitry, such that the output of theanalog correlator is with reference to the rising edge of the startpulse; and since the time duration of the delay line is known, the timedifference between the arrival of the start pulse and the output of thecorrelator can be added to the known duration of the delay line. Thisthen yields an accurate measurement of the actual optical path delay,and hence optical path length.

The analog correlator circuitry must employ very “fast” transistors, inthe sense that they must exhibit very fast output rise times whenexposed to very fast edge rates of an incoming signal. If thetransistors cannot maintain the fast input rise times (and in fact ifthe rise times of the excitation sequence pulses are not rapid as well),then highly accurate time resolution of the system will be impossible toachieve. Fortunately, very fast rise time transistors are available fromseveral modern integrated circuit technologies.

The analog correlator receives a “start pulse” which is transmitted downa reference delay line on the oximeter circuit board. Since the timedelay of the reference delay line can be created at printed circuitboard fabrication time and subsequently measured accurately withlaboratory equipment, the delay is well characterized and can bearranged to be shorter (or longer) than the approximate optical pathdelay through the tissue. The start pulse is used to initiate theoperation of the analog correlator, and then the identification timefrom the correlator can be added to the known delay time of thereference delay line. With the total propagation time of the opticalpath known, the physical path length can be determined accurately, whichin turn can be employed in the oximetry equations to yield absoluteoxygen saturation values. The correlator implementation must employ veryfast transistors to preserve the fast excitation sequence pulse risetimes. It is understood that several different implementations arefeasible, one example of which is the operation of a very fast “clockedcounter”, which is “started” by the arrival of the reference signal fromthe delay line, and “stopped” by the pulse output of the correlator,such that the number of “counts” in the counter, when it is stopped,yields a direct time measure of the propagation delay of the probe lightthrough the medium being measured.

As a beneficial enhancement of the analog correlator, rather thantransmitting the start pulse, a duplicate of the analog correlator isfeed directly from the excitation sequence generator. The two analogcorrelator circuits are operated in parallel. One of the correlatorcircuits processes the analog optical-to-electrical pulse train asdescribed above. The second correlator processes the electricalexcitation sequence pulse train that has arrived down the referencedelay line. The excitation sequence holding registers in bothcorrelators are loaded with the same excitation sequence pulse train.Thus, if a match occurs in one correlator, a match will occur in theother correlator as well. However, the two correlators will detect thematch at slightly different times, because the propagation delay in thecircuit board material can be designed to be shorter than thepropagation of the optical signal in tissue. Thus, the measured delay inthe reference line can be added to the difference in detection timesbetween the two analog correlators, to create an estimate of thepropagation delay, and hence the optical path length, in the tissue.

If it turns out that the optical path lengths at different wavelengthsare substantially different from one another, a “training” cycle may beused for each wavelength. That is, before beginning an actualmeasurement cycle, it may be necessary for each wavelength to be “pathlength-calibrated” by passing the start pulse down the reference delayline while each individual VCSEL is transmitting its unique excitationsequence code; when all VCSELs have done so, and the individual pathlengths have been measured and stored in the circuitry, then all VCSELscan begin transmitting for the purpose of making the material componentspecies measurements.

Hybrid Correlator.

With reference to FIG. 11, an embodiment is shown for determining theamplitude of the one or more wavelengths used to enable the comparisonof the transmission and absorption of one or more single-wavelengthlight sources as the source illuminates a measurement subject. FIG. 11presents the block structure of a mixed signal mechanism that can beused to collect this measurement. By adjusting the variable delay line1160, the mechanism shown is able to take a sequence of measurements;each measurement corresponds to a particular delay of light through themeasurement subject.

A coded data stream (zeros and ones) from the excitation sequencegenerator 1100 is used by the driver 1110 to modulate a light source1120, in the preferred embodiment a VCSEL. The VCSEL is coupled to themeasurement subject 1130 (the light shines on the skin), and a sensor1140 (e.g. but not limited to a photo detector) detects the light andits output is amplified 1150, producing a time varying analog signal.The correlation function 955 is performed by multiplying this timevarying signal 1170 by a time delayed version of the coded data stream.This multiplied signal is summed 1180 and the result captured at theoutput 1190 and processed in combination with other values produced byother time delays and the result of correlation with the excitationsequences used to modulate optional additional wavelengths of light.

The coded data stream takes two paths to the output; one is a simpletime delay through the variable delay line 1160, and the other is anoptical path through the material to be measured. The two paths arecombined at the multiplier 1170. The multiplier 1170 is configured suchthat the 1's of the excitation sequence from the delay line 1160 causethe multiplier 1170 to pass the analog signal from the amplifier 1150 tothe summing amplifier 1180. Alternately, 0's from the delay line causethe multiplier to invert the sign of the analog signal from theamplifier 1150 before passing that signal to the summing amplifier 1180.The result of this multiplication will be that the output 1190 of thesumming amplifier 1180 will have the correlated parts of the analogsignal amplified, and the uncorrelated parts of the analog signal willaverage over time toward zero.

Referring to FIG. 12, a representative circuit to perform thismultiplication is illustrated. A Gilbert cell multiplier using acombination of FETs 1230, 1235 and bipolar transistors 1250 for thedigital/analog mixing is shown, with a differential analog input1210,1215, a digital input 1220 and the compliment of the digital input1225, working together to produce, at the output 1260, 1265, a piecewiseversion or inverted version of the analog input.

The multiplier is configured so that I's from the delay line will causethe multiplier to pass the analog signal from the amplifier to thesumming amplifier, and 0's from the delay line will invert the sign ofthe analog signal from the amplifier before passing that signal to thesumming amplifier.

It is useful, for comparison purposes, to produce a ratio of thecorrelated response of the system to a reference. Two references may beprovided by versions of this invention. Using the switch to pass thesignal or to ground the input of the VCSEL Driver circuit, when theswitch is grounding the input, measurements may be taken while the VCSELis driven at a constant brightness. Also usefully, the VCSEL Driver canbe driven with a different excitation sequence input than the one usedat the multiplier. This approach will drive the VCSEL with a timevarying signal, but one that is not correlated to the detection circuit.This in turn will produce an output at the summing amplifiercorresponding to noise reduced by the processing gain of the excitationsequence correlation.

When multiple values for the delay are used, one per measurement period,multiple measurements can be taken. The set of specific delay valuesused, combined with the corresponding measured output of the summingamplifier, provide a set of points that can be used to characterize themeasurement subject. The impulse response of the optical leg of themeasurement system can be roughly described as answering the question,“When a bundle of photons are released into the sample at time zero, howlong does it take for them to reach the sensor?” While each photon takesa particular path through the sample, and emerges at only one time,there is a probability density function that can be used to predict theaggregate energy of the collection of photons at each of the likelytransit times through the sample. Referring to FIG. 13, an example ofthis impulse response is shown. There is a ballistic transport delayassociated with the geometric path length from time zero to the time thefirst photons appear at the sensor. No photons can reach the sensorbefore this minimum transit delay. There is a spike 1330 at the timeassociated with straight-line transit from emitter to sensor. The heightof this spike relative to the rest of the energy observed ischaracteristic of the material sample and the coupling to the material.If the material is a homogeneous diffusor, the observed impulse appearsas pictured, with paths distributed about a mean. For this example, agamma-like distribution was chosen to represent the tissue response.Generally, however, the tissue response may be represented by adifferent distribution.

While plotting the impulse response is not a necessary step inevaluating the measurement subject, a plot of the values collected canshow the correspondence between those values with the characteristicimpulse response (i.e., point spread function) of the measurementsubject.

As shown in FIG. 14, if the impulse response of the material is similarto that shown in FIG. 13, the measured points corresponding to thevarious delay line settings will lie upon the impulse response.

Based on the multiple relatively orthogonal coding sequences describedelsewhere, the sensor output can be processed simultaneously by multiplemixers to produce the measurement value for each input wavelength. Beingable to run these measurements simultaneously can advantageouslyeliminate motion-based inaccuracies when comparing multiple wavelengths.

In a heterogeneous measurement sample, the impulse response may not be asimple Gamma function. The longer the delay, the deeper the penetrationof the photons into a sample. Because the makeup of the sample may varywith depth, obtaining the relative absorption of multiple wavelengthsover a range of depths provides an advantage in characterizing themeasurement sample.

Multiple mixers/correlators, each fed with the time delayed excitationsequence input of a reference from a VCSEL driver input, can be runsimultaneously to obtain ratios between wavelength extinctioncoefficients, thereby providing for material sample characterization. Asshown in FIG. 14, where two impulse responses are shown, notionally fortwo different wavelengths, A₁ is the signal level measured forWavelength 1 at delay A, and A₂ is the signal level measured forWavelength 2 at delay A. In a similar fashion, B₁ and B₂ correspond todelay B, and C₁ and C₂ correspond to delay C. The two curves composed ofthe varied delays are not simply scaled versions of each other, butrather

$\begin{matrix}{\frac{A_{1}}{B_{1}} \neq \frac{A_{2}}{B_{2}} \neq {\frac{A_{3}}{B_{3}}.}} & (48)\end{matrix}$

The result of this is that the various component concentrations can becalculated to vary by path length.

Multiple mixers, each fed with a time delayed excitation sequence usedas the input of a different VCSEL driver, can be run simultaneously toobtain ratios between wavelength extinction coefficients, with the delayused in the correlation providing multiple ratios, each corresponding toa notional “depth” in the measurement sample. This provides a set ofselectable sites for measuring the variation of multiple properties.

Analog Delay Correlator-B.

With reference to FIG. 20A, a resettable integrating amplifier 2000,alternately called a summing amplifier 1180 is shown as a schematicrepresentation of an example implementation. The function of thiscomponent is to, with equal weight over time, collect the analog signallevel of its input 2001, and provide the integrated value over time asoutput 2008. This integrated value isOutput₂₀₀₈=∫_(T) ₀ ^(T) ^(RESET) Input₂₀₀₁ dT  (46).

The circuit is created using an operational amplifier 2003, connected topower 2006 and ground 2005 sources to provide energy for operation.Operational amplifiers have two inputs, an inverting input 2002 and anon-inverting input 2004. The input 2001 is connected to the invertinginput 2002 of the operational amplifier 2003. The output of theoperational amplifier 2007 feeds both the output 2008 of the integratingamplifier subsystem system 2000 as well as terminals on the feedbackcapacitor 2009 and the source/drain of reset field effect transistor(FET) 2011. The other terminal of the feedback capacitor 2009 and thedrain/source of reset field effect transistor (FET) 2011 are connectedto the inverting input 2002 of the operational amplifier 2003. When theRESET pin 2010 is not asserted, the feedback capacitor isolates theoperational amplifier input from the output and stores an integratingcharge, providing the integrating function over time. When the resetinput 2010 is asserted, the FET 2011 conducts allowing the feedbackcapacitor 2009 to quickly discharge. The output 2007 of the operationalamplifier 2003 and both terminals of the feedback capacitor, take on theinstantaneous voltage level of the input terminal 2001. The integratingamplifier subsystem 2000 is used in the schematic in FIG. 20B as part ofthe analog delay correlator-B.

In this circuit, differentiated from the solution described related toFIGS. 18A and 18B, rather than having a series of delays 1812 in theanalog received signal 1810, instead, in this implementation, the inputsignal 2016 is replicated with balanced connections that provide thesame delay 2015 and is fed simultaneously to each buffer 2020 such thata timing-identical copy of the input 2016 is presented to eachmultiplier 2025. Said multipliers may be implemented as the exampleGilbert mixer of FIG. 12 with buffer 2020 output as the analog input1210 in mixer circuit 2025. Other mixers and multipliers known in theart may be used. The digital input of the multiplier 2025 is derivedfrom the reference signal which is a precisely time imaged version ofthe signal produced by the excitation sequence series of pulses. Thissignal is the input at buffer 2030, and is distributed to each of amultiplicity of multipliers through controllable time delay elements2040. In the figure, these time delay elements are shown connected incascade form. The inputs to the time delay elements can alternately beconnected in parallel or a combination of parallel and cascade. Theadvantage of using the cascaded delay for the excitation sequence isthat the signal is a binary full swing signal that can be regenerated ateach delay stage without injecting any additional noise into the signalderived at the output of the multiplier 2026. The excitation sequenceseries of pulses serving as reference for each group of elements is timedelayed such that it corresponds to a desired point in the impulseresponse for each group. The result of this series or collection ofmeasurements is what is represented in the impulse response graphs shownin FIGS. 8C, D and E, FIG. 13. FIG. 14, FIGS. 24A and B, and FIG. 28.The production of this value is performed for a single time point alongthe full impulse response curve by the circuit illustrated in FIG. 11,based on the setting of the variable delay line.

Each multiplier output 2026 is the input to the resettable integratingamplifier 2000. The output of the integrating amplifier 2050 providesthe input to the analog to digital converter (A/D) 2055. The reset input2010 and the sample input indicated by the clock triangle on the analogto digital converter 2055 are operated to sample the extractedcorrelated signal level, convert the value to digital form 2060, andlatch the value in the output register 2070, responsive to a clocksignal feed to the data latch 2070. The value latched in register 2070is a surrogate for the inverse of the attenuation, that is, it is thereceived signal strength. The value latched in the register may be readand averaged over time to provide more precision of measurement, as thenoise from the quantification, the noise for the system, and if theexcitation sequence is sufficiently long, the correlation noisecontinues to be reduced. For this reason, an implementation of thepresent invention uses an A/D with fewer bits of precision than areachieved in the end, and very long excitation sequences, sampling thereceived signal strength multiple times during the course of theexcitation sequence, thereby increasing the available SNR beyond therange of the A/D. The full correlator subsystem 2075 can be incorporatedinto the larger system shown in FIG. 21.

Referring to FIG. 20C, the relative timing of the control signals thatare used to compose the apparatus of FIG. 20B into a system aredescribed. The control signals include the integrator reset input 2010,the A/D sample clock on latches 2055, and the data register latch clockof 2070. These are shown in the relative timing diagram of FIG. 20C. Thetiming traces from top to bottom are the A/D sample clock 2056 and thedata register latch clock 2071, and the excitation sequence repetitionperiod are shown. Using time as the x-axis, the positive going edge ofthe control signal is designated as the time that the event istriggered. The width of the triggering pulse is dependent on theparticular technology chosen. Following design practice that includevalidation of set-up and hold times is required and assumed. The “S”curves, for example, the one going from the pulse on timing arc 2056 tothe pulse on timing arc 2071, are used to indicate that the trigger asthe beginning of the arrow must be complete before the next event, thetrigger at the end of the arrow. For example, the A/D sample 2056 mustbe complete before the integrating amplifier 2011 is reset, or the dataat the output of the A/D is latched 2071. The data at the output of theA/D must be latched 2071, before the integrating amplifier 2011 isreset. In one embodiment, the sequence of triggering multiple channelsof attenuation measurement are concurrent, and are triggered from thetransmitter indicating the start and end of excitation sequences ofknown and complete statistics 2080, that is, that the signal is balancedsuch that the ‘chip’ to ‘chip’ interference in the impulse responseaccumulation is canceled out. Time period 2081 encompasses the sequenceneeded to capture data from the correlating subsystem. This sequence isrepeated during the period indicated by 2083. In the intervening timerepresented by the ellipsis of 2082, this sequence can be repeated manytimes to improve the accuracy of the measurements. The full duration ofthe measurement, including the repetitions of 2082, are advantageouslyencompassed by a set of excitation sequences that do not repeat over theperiod 2084 plus the period 2085.

System Composed of Multiple Correlators

With reference to FIG. 21, in one embodiment, a collection of thecorrelators 2075 is arrayed sharing the received signal 2016. Eachsub-system however has the excitation sequence reference correspondingto the propagation amplitude that it will quantify. The inputs 2030 eachare unique for the correlating subsystem 2075. As described in FIG. 19,the output of the correlator may have residual non-linearities. Element1900 in this figure provides a digital mapping function to restorerelative linearity to the collection of A(λ_(i)) to the matrix system2100. The result is an output of the concentration levels of thematerial components 2120 for which the system has been designed.

Multi-Dimensional Measurements

The techniques disclosed herein are not strictly limited to single-pointor one-dimensional measurements. Rather, by altering the modulationscheme, additional information about the measurement as a function oftissue depth and distance between the light source and the detector ispossible, thereby resulting in a two-dimensional measurement. Extendingthis concept by employing a grid or array of emitters and detectors canprovide truly 3D measurements of tissue characteristics. Acquiringresults of the measurements at different points in time (e.g. throughthe heart cycle or as a function of an appropriate signal-modulationscheme), adds a temporal component to the measurement, in turn resultingin a 4D measurement (three spatial dimensions and time). Carrying themeasurement at multiple wavelengths provides an additional fifth,spectral dimension. The extension from one dimensional physiologicalmonitoring to five dimensions of monitoring greatly enhances thecapabilities and reach of embodiments of the invention presented in thisdisclosure. Applications of the disclosed embodiments include importantand still unresolved clinical problems such as, e.g., oral cancers thatcan be identified when stimulated with light having appropriately chosenspectrum, or abnormal protein buildup in utero that can be noninvasivelymonitored by means of optical techniques.

Embodiments of the present invention provide a method for transmitting,receiving, and quantifying two or more concurrent signals. The methodincludes transmitting a first signal in response to a first excitationsequence, at a first wavelength and from a first location whiletransmitting a second signal concurrent to the transmission of saidfirst signal, in response to a second excitation sequence, at a secondwavelength and from a second location. The method additionally includesreceiving the transmitted signals as an aggregated signal; andprocessing the aggregated signal by correlating it with atime-controlled representation of the first excitation sequence and,additionally or in the alternative, by correlating the aggregated signalwith a time-controlled representation of the second excitation sequence.The quantified result(s) of each of the first and second correlationsprovides measure(s) of, for example, attenuation and scattering of thefirst and second signals through the media of propagation. The media mayinclude, for example, at least one of turbid media, molecular dye, anorganism, and similar spectrally distinguishable material.

A system of the invention is then used to determine media-characterizingproperty values based on relative and absolute results of the magnitudeof the received signal components (such as variations over time). Themedia-characterizing parameters include, for example, hemoglobin analyteconcentrations, blood glucose level, blood protein level, blood lipidlevel, total body water, central blood volume, heart rate, heart ratevariation, concentration of a fluorescent biomolecule, and concentrationof a fluorescent dye.

Useful variations of the composition of the excitation sequences includetime periods of silence, possibly with the arrangement of the silentperiods such that a received signal can be calibrated in the absence ofany other system excitation, the repeating duration of the variousexcitation sequences may be different, the excitation sequences maybeneficially comprise codes substantially balanced in duration betweenmaximum and minimum energy, the excitation codes may be based onpseudo-random repeating sequences, the excitation sequences can be basedon actual random processes, the excitation sequences can be producedfrom stored sequences. The excitation sequences of the multiple channelsare beneficially relatively orthogonal, that is, the long-termcorrelation of one of the excitation signals with another excitationsignal will not be substantially different from correlation to thebackground noise of the system. If the excitation sequences are producedwith a repeating pattern, it may be beneficial to provide a longerbaseline for noise reduction by using different pattern lengths ordurations. Durations that are relatively prime one to another willproduce longer independent sequences, thereby reducing any cross-inducednoise from the received and correlated result. If an excitation sequenceis produced digitally by a clocked process, that is, the sequence can berepresented by a sequence comprising a fixed set of levels, with theoutput level changing at time periods which are a multiple of a clocksignal on a fixed time base, the changes in level of the excitationsequence can be separated in time by at least the nominal impulseresponse time of the sample being measured; however with the mixing andcorrelation techniques described herein, the time period separatingchanges in level of the excitation sequence can also be substantiallyshorter than the nominal impulse response time of the sample beingmeasured and still provide accurate quantification of properties of thesample being measured.

The transmitter or the receiver, or both may be juxtaposed to a sampleto be characterized; said sample altering the signal between thetransmitter and the receiver. The signal in transit may be altered by atleast one of turbid media, molecular dye, an organism, and similarspectrally distinguishable material. The time-controlled representationof the excitation sequence may be produced responsive to the excitationsequence passed through a variable delay element. This method may beused in a system where the receive processing includes scanning thedelay or adjusting the delay element over an interval—thereby providinga representation of the impulse response of the intervening material.The delay element can also be set at a determined delay to correspond tothe measurement of interest. The setting of the delay elements can beaccomplished in a manner that synchronizes the delay used for eachcorrelation process, or the time-controlled delays of the multiplecorrelation processes may be scanned or set independent from each other.Alternately, an array of similarly constructed multipliers andintegrators can be connected, each to a varied delay of the excitationsignal. In this case, multiple values of the impulse response can beproduced concurrently. The transmitting device may usefully be, forexample, an antenna, a laser, a VCSEL, an LED, or similar device. Thedetector of the receiver may usefully be, for example, an antenna, a PINdiode, an avalanche photo diode, a photo multiplier tube, or similardevice. The correlator may usefully include, for example, a mixerfollowed by a summing amplifier, followed by an analog to digitalconverter; a set of delay elements interconnected with a representationof the excitation sequences, a low-pass filter, an analog-to-digitalconverter, and a threshold detector and counter; analog delay linecorrelator; or similar circuits. The logical unit for processing theaggregated signal correlator may usefully include, for example, aprogrammable device, a device programmed with software, a fixed(non-programmable) logic unit, an analog circuit element, or combinationof such elements.

Embodiments of the present invention may include an apparatuscomprising: two or more transmitters, each transmitting its signal inresponse to an excitation sequence, at a unique wavelength or from adistinct location; said transmission being made at times concurrent tothe one or more other transmitters. Said apparatus also comprising areceiver for receiving some or all transmitted signals as an aggregatedsignal; and a logical unit for processing the received aggregatedsignal, whereas: a distinct correlator is used for quantifying thereceived portion responsive to each transmitted signal; said aggregatedsignal is correlated with a time-controlled copy of the excitationsequence to generate a data sequence; and the data sequence isquantified uniquely corresponding to the transmitted signal of interest.The time-controlled representation of the excitation sequence canalternately be produced by passing the original excitation sequencethrough a calibrated time delay element, such as tapped delay line orother similar technique, or it may be created independently at thereceiver, based on a calibration pulse received from the transmitterassembly.

Making Use of the Fine-Grained Impulse Response

Making use of the fine-grained impulse response provides increasedfidelity of measurement in several systems. In addition to themulti-transmitter, multi-receiver, and micro-mirror-based localizationcontrol techniques, an embodiment of the current invention includesusing information about the time delay (and hence radiation path length)contained in the fine-grained impulse response to increase the precisionof spatial localization of the sample's features that are beingcharacterized. Data processing based on the fine-grained impulseresponse can increase the contrast of the details of interest containedin the received data of the impulse response. This contrast increase, asdescribed below, uses the fine-grained detail of the impulse response tonarrow the analysis to the areas or regimes of interest in the sample.The improved fidelity provides greater precision of measurement, and,with appropriate calibration, greater accuracy. The greater accuracy, inturn, facilitates making finer distinctions between material componentsand/or their properties.

Making use of optimized choice of wavelengths of operation of thespectrometric device, as discussed above, improves the propagation ofvariance in the calculations that determine material properties Z (see2110 in FIG. 21) based on the results of the amplitude measurements (see2070 in FIG. 20B).

According to the idea of the invention, the analysis of datarepresenting the measurement of a parameter characterizing the tissue(for example, a concentration of a particular species such as oxygenatedhemoglobin, or glucose, or specific protein, or a physiologicalparameter such as aortic pressure) is performed in such a fashion as todynamically and spatially localize a portion of the tissue to which suchdata pertains. Such spatial localization is achieved by analyzingportions of a curve that represents the parameter of interest and thatis obtained by accounting of the impulse response of the tissue.Specifically, and in further reference to FIGS. 8C and 8D, whiledifferent portions of an impulse response curve of the tissue areassociated with different paths that the EMW is taking from thetransmitter through the tissue to reach the receiver, the impulseresponse curve integrated over the area of the tissue through which theEMW propagated is associated with the averaged over the tissuecross-section path length for the EMW (which is, of course, wavelengthspecific).

The tissue is excited by an impulse of the EMWs (at multiplewavelengths) by at least one transmitter and the impulse response of thetissue is measured, by at least one receiver, in terms of a chosenoptical characteristic (e.g., absorption, scattering, opticalanisotropy). Referring to Eq. (37c), such measurement produces thewavelength-dependent wavefront parameter denoted as [A(λ_(i))]. Usingthe pre-determined characteristics such as the molar extinctioncoefficients [ϵ(λ_(i),C_(i))] and estimating the path-lengths [δ^(Ci)]of the EMWs through the tissue, the concentrations [MC(λ_(i))] of thespecies (material components) are then calculated based on the measuredimpulse response curves. FIG. 24A illustrates an example of the methodof the invention. As shown, the use of the impulse response of thetissue sample (that is determined based on data representing attenuationof the signal transmitted through the sample) is made. The three graphson the left (2401, 2402, and 2405) represent impulse response curvescorresponding to time-dependent amplitudes of the signal as received atthe detector and associated with a particular wavelength and/or locationof a transmitter with respect to the sample. The amplitudes 2401, 2402,2405 are associated with the respectively corresponding uniqueexcitation sequence series of pulses used at the transmitter, and arerecorded as a function of time lapsed after the moment of transmissionof the corresponding signal. That is, the abscissa on each of the graphsrepresents the duration of time that a particular signal takes topropagate from the transmitter to the receiver. If the amplitude of theoutput generated by a correlator of the system is the result of a singlepulse, then the single pulse must have enough energy to distribute overthe time period and be measurable in a statistically repeatable manner.In an embodiment of the invention, many pulses are aggregated togetherby one of a variety of correlation circuits described above. However,other techniques could be used as long as the polarity of a pulse in aseries transmitted pulses is normalized (in the examples above bymultipliers, multiplexors, or exclusive-OR gates), and the time-delaycorresponding to each pulse is normalized based on the transmission timeto allow for correlation among different pulses from the series. Thistime-delay is directly related to the length of the path of propagationof a pulse and indirectly related to the transit speed. The value shownin the extinction curves of, for example, FIG. 3C, is in the units of(1/[concentration*distance]), meaning that to determine theconcentration takes a distance factor (nominally the delta delay) and anattenuation, or received power/transmitted power. In further referenceto curves 2401, 2402, and 2405, values of amplitude at a chosen“impulse-response delta time” for the respectively corresponding decodedexcitation sequences (λi) 2411, 2412, and 2415 are used in a dataprocessing step corresponding to the operation shown as 2420 (see Eq.37C).

The output from the system of the invention at each time step is aseries of values representing material component concentrations. Thesevalues, when aggregated together in a time-dependent fashion, are shownby the graphs 2421, 2422, and 2425 of FIG. 24A that represent,respectively, impulse responses for component concentration1, componentconcentration 2, and component concentration N. The ellipses are used inFIG. 24A to indicate that the data collection and correspondingmeasurements are performed at multiple wavelengths in a similar fashion.

When referring to λi (with respect to data), the cohort referencedincluded both those correlated impulse responses that are differentiatedby wavelength, and those correlated impulse responses that aredifferentiated by transmitter or receiver location. All of these variousλi's are each extracted from the received signal by correlating them atime-controlled representation of the source excitation sequence ofpulse.

A further improvement made possible by the fine-grained impulse responsedata is the potential to adjust the linear combination of the λi'simpulse responses based on the diffusion-varied speed of propagationthat is dependent on wavelength. This property of the sample materialwas discussed above in reference to the optimal selection of wavelengthsof operation of the spectrometric device. If this information isavailable, then it can be used by modifying the process described inFIG. 24A by appropriately time-shifting the samples taken from theimpulse response based on wavelength—for each particular analyte—whencreating the instantaneous material property of the impulse response.

For example, and referring to FIG. 24B, graph 2441 represents the timeshift appropriate for wavelength λ₁ specific to material component1. Thedesired time-parameter 2436, representing the “impulse response timedelta” in the concentration of material component 1 output, is used asindex 2451 on graph 2441, to produce an adjusted time delta 2461. Theadjusted time delta is used to index 2411 on graph 2401 to provide inputto the matrix calculation 2420.

In a similar manner, graphs 2442 . . . 2445 represent the time shiftappropriate for wavelengths λ_2 . . . λ_N specific to the materialcomponent1. The desired concentration output impulse response time delta2436 is used as indices 2452 . . . 2455 on graphs 2442 . . . 2445, toproduce the adjusted time deltas 2462 . . . 2465. These adjusted timedeltas are used to index 2412 . . . 2415 on graphs 2402 . . . 2405 toprovide inputs to the matrix calculation 2420.

Although this has been illustrated as a simple time shift, it will beappreciated that the input could also be a function of the wavelengthimpulse response, taking into account more points, and providing for thepossibility of weighted summation, interpolations, and other similarmethods of choosing the input to the matrix operation 2420. Additionallyit is noted that the term [δ^(BM)]⁻¹, labeled as 2430 in FIG. 24B, isused to account for the path length of each particular data sample. Asdiscussed above, this may be a scaling for the increased extinctionexperienced by the path; it may as well be augmented by gained knowledgeof the heterogeneous makeup of the sample. This knowledge may be gained,for example, from the sampled-measurement of the sample with the use ofthe micro-mirror array described above, from results ofexternal/independent analysis using the micro-mirrors in combination toprovide feedback, optimizing signal transmitted or received by modifyinglevel, frequency of chips, and excitation sequence characteristics,providing a larger field of view, increasing the effective density oftransmitters and receivers, and techniques described above. Allinformation is contained in the signal from a single transmitter andreceiver with or without a mirror or mirrors in place. Additional usefulinformation is gained by increasing the number of transmitters andreceivers or by increasing their effective density as can be done withthe micro-mirror array described earlier.

The Use of Embodiment(s) for Determination of Beat-to-Beat CardiacStroke Volume and Cardiac Output

The magnitude of cardiac output in cardiac-compromised patients (e.g.,those with damaged heart valves, or those who has a heart attack orcoronary artery disease, or who have been heavy smokers, etc.) is acritical measure of health and illness. The units of cardiac output are“liters per minute” of blood pumped by the heart. At a finergranularity, cardiac output is composed of the individual amounts ofblood pumped during each heartbeat, referred to as “beat-to-beat strokevolume,” where the “stroke” is one heartbeat or contraction of theheart's left ventricle. In the early days, before the onset ofhospital-based ultrasound systems, there was no reliable way ofmeasuring beat-to-beat stroke volume; an X-ray based technique couldmeasure stroke volume, but only a few medical centers had thecapability, and a cardiac catheterization laboratory and staff wererequired (performing a highly skilled procedure). Modern ultrasoundsystems can estimate beat-to-beat stroke volume if used by a highlytrained ultrasonographer. Although this is a costly test, the cost isless than the cost of a cardiac catheterization laboratory.

In Continuous Determination of Beat-to-Beat Stroke Volume From AorticPressure Pulses in the Dog, Circulation Research, 39(1):15-24, July1976, Bourgeois, M. J., B. K. Gilbert, G. von Bernuth, and E. H. Woodproposed a method for calculating a beat-to-beat stroke volume (andhence cardiac output) based on the beat-to-beat pulse pressure wave inthe central thoracic aorta. The proposed method raised at least twoquestions. First, according to the method, the aortic pressure pulsewave had to be measured with a catheter inserted through a femoralartery and advanced into the aorta. Such a procedure requires theavailability of specialized hospital environment. Second, to be ofclinical value, the calculations had to be carried out on a beat-by-beatbasis, in “real time,” which was beyond even the mainframe computers atthe time of the research reported in the manuscript. Accordingly, theproposed method was never implemented in general clinical practice.

Non-trivial extensions to the work of Bourgeouis et al. that make use ofthe oximetry aspects of the system herein described can be used toprovide accurate non-invasive measurement of the beat-to-beat strokevolume of a measurement subject. The described spectrometer providesmeasurement of the near-instantaneous changes in blood componentconcentration between the transmitter and receiver of the sensors; thewaveform of this information is used as a “surrogate” waveform for theaortic pressure pulse contour. The quality of this surrogate waveform isfurther improved where a subset of the material component concentrationcurve is selected from the spacing and depth region of the sample, whichprovides a maximum variation or the oxygen bearing components, ratherthan the variation of the oxygen-depleted components. This selection ofthe source for the saturation waveform provides for measurement of thearterial rather than venous variation. Choosing a source correspondingto an artery (not a vein) close to the skin surface provides a bettersurrogate for the aortic pressure pulse: the arterial pressure pulse isa “damped” version of the aortic behavior, since the peripheral arteryis fed from the aorta. Using this local resolution capability and theselection criteria, the system computes beat-to-beat cardiac strokevolume, and hence cardiac output, in real time, in the patient'sfree-living environment. The resulting information provides benefit tocardiologists treating patients whose cardiac function is compromised.Further, the use of multiple sites of interrogation (i.e., multiplephoto detectors) time-locked together, with or without the previouslydescribed excitation sequence encoding, would allow the detection ofvariations in blood flow and cardiac output associated withcardiovascular state, e.g. hypovolemia detection. Any changes in thevalues of these variables then may be predictive of deterioration orinstability of cardiovascular state.

What is provided by the present invention is the ability to identifyand/or localize the ideal measurement site, to measure instantaneouschanges in oxygen saturation and to create a surrogate aortic pressurepulse waveform, and thus calculate cardiac stroke volume and cardiacoutput.

Properties that are beneficially measured using the optimized selectionof wavelength, the introduction of challenges to the sample, and the useof selected possibly relatively orthogonal selected transmitterexcitation sequences include, but are not limited to the varioushemoglobin analyte concentrations and time based variation described, aswell as blood glucose level, blood protein level, blood lipid level,total body water, central blood volume, concentration of fluorescentbiomolecules or dyes, respiratory tidal volume, beat-to-beat cardiacstroke volume, and beat-to-beat cardiac stroke volume variation.

Based on the multiple relatively orthogonal coding sequences describedelsewhere, the sensor output can be processed simultaneously by multiplecorrelators to produce the measurement value for each input wavelength.Being able to run these measurements simultaneously can advantageouslyeliminate motion-based inaccuracies when comparing multiple wavelengths.

Challenge Based Measurements

Material properties of interest, characterizing the tissue, includephysiological parameters some of which (for example, central bloodvolume or total amount of water in the tissue) cannot be measureddirectly. Diagnostic techniques, developed to measure such physiologicalparameters indirectly are based on the dilution of “tracers” introducedinto the blood and tissue. For example, tracers are injected into anartery or vein, administered orally, or inhaled (in a form of, forexample, small quantities of a gas such as carbon monoxide), and theblood sample is subsequently analyzed to quantify the percent of thetracer therein. The tracers are selected to facilitate their detectionthrough various in vitro (laboratory) or in vivo (in the humans) means,but in this context, spectrometrically. Some tracers that have been usedinclude carbon monoxide, and fluorescent dyes, several of which areavailable from commercial sources and have been or are being tested forsafe clinical application in patients. Depending on the tracer used, andthe accuracy of measurements, various physiological variables can bemeasured.

Other material properties of interest are associated with variousphysiological analytes or metabolites. A subset of these metabolites arethose chemical derivatives that are formed in vivo through normalbiochemical processes in response to the administered drug (as describedabove) that has been metabolized (i.e., broken down or degraded throughin vivo biochemical processes). A chemical substance is a materialcomponent with uniform properties. Given the uniform properties, achemical substance may be analyzed spectrometrically to determine itsunique spectral signature. Using techniques described elsewhere in thisapplication, specific EM wavelength(s) are chosen to best identify themetabolite that is associated with the drug of interest. In this manner,the metabolite is used as a “tracer”.

For the purposes of this disclosure, the term “challenge” is used torepresent a substance introduced into a material (e.g., blood) to laterbe measured as a proxy for a parameter associated with the material(e.g., blood volume or drug uptake). Examples of challenges include butare not limited to trace gases and drugs. Some drugs could be detecteddirectly and/or through their associated metabolite(s). Given achallenge, the process used in analysis of the material parameter isillustrated in the flowchart of FIG. 22. The process starts at steps2200 and 2210, with a baseline measurement of the sample. Once thebaseline of a sample characteristic is established, the challenge isintroduced (e.g., through injection, inhalation, or oral introduction)at step 2220 and the so-modified sample is then measured to acquirerepresentative data at step 2230. This process of data collectioncontinues until sufficient (as measured by an externally defined rule)data are gathered, at step 2240, to allow for the calculation of theproperty value of interest at step 2250. At this point the process iscompleted (END 2260).

The techniques used to-date to quantify the properties, of a materialcomponent of a sample, that are caused by a challenge, have varyingdegrees of accuracy, invasiveness, and, at times, undesired sideeffects. What is required in clinical applications is a means to measuresuch quantities with sufficient accuracy for clinical requirements,non-invasively, and with minimal side effects to the patient. Thetechniques for selection of optimal EM wavelengths, transmitting thewaveform(s), receiving the waveform(s) after interacting with thesample, and quantifying the result that have been described elsewhere inthis application are as applicable to the challenge-based process asthey are to the detection of other properties.

As an instance of the challenge-based property measurement techniquethat is improved in substantial and non-obvious ways, the measurement oftotal blood volume can be augmented with the additional accuracyprovided by the invention to change the use case. Support of measurementof the change in blood volume over time would be of significant benefitin the triage of care. Carbon monoxide (CO) is used as a detectabletracer for the measurement of total blood volume, because it bindspreferentially over oxygen to hemoglobin.

Burge and Skinner, in “Determination of hemoglobin mass and blood volumewith CO: evaluation and application of a method,” Journal of AppliedPhysiology, vol. 79, pp. 623-631, Aug. 1, 1995, described measuringtotal blood volume by introducing a measured aliquot of CO to a subject,with carboxyhemoglobin concentration measures before and 10 minutesafter the administration of the aliquot. Their method used the apparatusillustrated in their report, and illustrated here in FIG. 23. In C. J.Gore, W. G. Hopkins, and C. M. Burge, “Errors of measurement for bloodvolume parameters: a meta-analysis,” Journal of Applied Physiology, vol.99, pp. 1745-1758, Nov. 1, 2005, the effectiveness of the techniquecompared to other techniques was demonstrated with a meta-analysis ofseveral studies. The accuracy of the technique is in the range of 0.8%,but achieving such accuracy requires the introduction of CO in theamount sufficient to raise the percentage of hemoglobin bound to CO by(or on the order of) 6.5 percentage points. As the dangerous limit of COconcentration in blood is ≤15%, and the half-life time of CO in bloodvaries from 50 minutes with pure oxygen to 2-8 hours in normal subjects,it is understood that the use of multiple detectable aliquots is notpossible within a short time period using the current technology.According to Gore 2005, the two reasons behind using such a largealiquot are, first, delivering a smaller amount may be subject todelivery quantity error, and second, that current external reading bloodgas monitors provide precision only to 0.1%.

With the increased precision offered by the methods and systems of thepresent invention, a much smaller aliquot can be used to provide usableand reliable results. Using aliquots low enough such that regularadministration of CO aliquots is safe over a period of time provides,for example, an effective monitor for central blood volume withoutcompromising the health of the subject Given the rate of rapid rate ofdispersion of CO through the system, it would be possible to presentadditional aliquots as often as every 10 minutes, provided theincremental detection is sufficiently sensitive.

Fluorescence

Local anesthesia is, in many cases, a superior means of pain control inavoiding the side effects of narcotic pain medication. Nerve blocks arecommonly used for postoperative pain control as well as the primaryanesthetic when the avoidance of general anesthesia is desired. The useof local anesthesia in this application renders an area insensitive topain without affecting consciousness, speeding the detection ofcomplications and patient recovery. However, when performing a nerveblock, relatively large drug doses can be required to inhibit the nerveconduction to the surgical site. By way of example, the requiredinjection volumes are often in the range of 20-40 cc. Unfortunately, aslittle as 1 cc of local anesthesia injected into an artery can causeseizure, and large intravascular local anesthesia injection will blockthe conduction system of the heart, possibly leading to cardiac arrestand death. Symptoms evidencing the onset of toxic anesthesiacontamination include numbness of the tongue, lightheadedness and visualand auditory disturbances, and can occur at concentrations below 7.5mcg/ml. Approximately 5 in 10,000 surgeries report severe complicationscaused by intravascular local anesthesia, but thousands of regionalblocks are done every day.

There are very few tools available to anesthesiologists to protectpatients from this problem. One known method involves an epinephrinemarked ‘test dose’ of local anesthesia. During injection, doctors watchfor an increase in heart rate or blood pressure indicating aninadvertent intravascular injection. Unfortunately, this method is veryunreliable. Patients taking beta-blockers or wearing a pacemaker may beunresponsive to epinephrine. Children often have blocks done undergeneral anesthesia that alters the body's response to adrenaline andmakes it difficult to determine if the local anesthesia isintravascular. An epinephrine marker during an epidural can causecomplications by decreasing blood flow to the fetus in pregnant women.There remains, therefore, a continuing need for improved systems andmethods for providing accurate and early warning of intravascular localanesthesia.

One embodiment of the invention is a method for operating a spectrometerto monitor a patient receiving local anesthetic marked with dye thatabsorbs infrared light. This method includes: (1) applying light pulsesto a patient that is receiving (or has received) anesthetic marked withfluorescing dye; (2) detecting the light emitted by the patient andproviding an output representative of the light detected; (3) processingthe output to derive information representative of the presence ofdye-marked anesthetic; and (4) displaying the information representativeof the presence of the dye-marked anesthetic as a function of theinformation representative of the fluorescence.

It is appreciated that results of applying such method that has beenenhanced by using the mechanisms described above in producing theimpulse response are advantageous. Due to the delaying nature of thefluorescing activity, the impulse response curve for any of thewavelengths that trigger the fluorescence will have extended tails. FIG.28 illustrates the impulse response of the tissue without thefluorescent dye (curve 2810) and with the presence of the fluorescentdye (curve 2820). The use of the method for operating a spectrometerthat utilizes the impulse response approach is advantageous over a peakor bulk correlation in that it increases the sensitivity of themeasurement.

Localization

Referring to FIG. 1C, an improved measurement of a cardiac stroke volumecan be enabled based on an embodiment of the invention. FIG. 1A shows anembodiment with the sensor 160 juxtaposed to the finger. This is oftenreferred to as a peripheral measurement and is appropriate for manyuses. However, when considering a measurement of the cardiac activity,the measurement of the pulse wave at the finger is strongly affected bythe impedance of the circulatory system. The waveform of the pulse wave(when measured at location 140) closely resembles, in shape, the curve2530 of FIG. 25 discussed below, with potentially undetermined noiseadditionally present in the heart-centric measurements. Therefore, it isadvantageous to ensure alternate placement of the spectrometric deviceonto the body. FIG. 1B shows the sensor placed closer to the heart, ascompared with FIG. 1A.

FIG. 1C includes a sketched side view of the human head, neck,shoulders, and a portion of the arterial circulation system 190. Theneck 193, just at the shoulder joint, provides proximity to at least onemajor artery, the carotid artery 191. A sensor 145 in this illustrationis an instance configured to monitor closely the concentrationcharacteristics due to the blood flow in the carotid artery 191. Whenelements of an implementation of the invention are taken together(including at least the transmitter(s) 110, and the receiver(s) 130,131, and 132) the spectrometric system is referred to, collectively, asthe sensor 145.

Referring to FIG. 25, the graphs 2500, 2510, 2520, and 2530 illustratevariation of blood pressure in a blood vessel over time. Systole is theportion of the complete cardiac cycle when the left ventricle isejecting blood into the ascending aorta through the aortic valve.Diastole is that portion of the cardiac cycle when the ventricle isbeing refilled by the left atrium, the aortic valve is closed, no bloodis flowing into the aorta from the left ventricle; however, duringdiastole, blood continues to flow from the aorta into the arterialvasculature and then into the arterioles and capillaries deep in thebody tissues. During the systole, blood is also flowing into theperipheral arteries as noted above; but the primary effect is thepartial refilling of the aorta by the left ventricle.

Referring to graph 2500, the point at (time=t_(D), pressure=P_(D))represents the minimum pressure at the end of diastole 2501, just as thesystolic ejection phase begins, as measure at a location point specifiedwith respect to the heart. The systole continues, achieving a peakpressure P_(S) at 2502. As the left ventricular contraction ends and theaortic valve closes, the aortic pressure begins to decrease to a localminimum at 2503, followed by a slight pressure increase at 2504(referred to as the incisura or dichrotic notch). The dichrotic notch iscaused by a pulse wave reflection traversing the ascending aorta, theaortic arch, and the descending thoracic aorta as the aortic valvecloses. At the time following the appearance of the dichrotic notch at2505 in FIG. 25, the aortic valve is completely closed, and the aorticpressure decreases as blood drains out of the aortic, which functions asan elastic reservoir to maintain as much pressure in the arterial treeas possible. It has been shown (in Continuous Determination ofBeat-to-Beat Stroke Volume From Aortic Pressure Pulses in the Dog,Circulation Research, 39(1): 15-24, July 1976, Bourgeois, M. J., B. K.Gilbert, G. von Bernuth, and E. H. Wood) that the phase 2505 can becharacterized by an exponential decay with time constant τ. This timeconstant, which can be determined during the diastolic portion of eachcardiac cycle, is used in the calculation of beat volume. The beatillustrated ends at 2506 with the beginning of the next cardiac systoleas the aortic valve reopens.

As the pressure wave is measured at additional points across the bodythat are located progressively farther from the heart than the point atwhich the pressure wave graph 2500 has been measured, the correspondingtime-dependent pulse waves (shown as 2510, 2520, and 2530) weaken due toimpedance of the circulatory system. The peaks of the correspondingpulse-waves are generally reduced, the steady-state pressure increases,and the exponential decay portions of the curves are also affected. Ingraph 2530 (acquired from the measurement taken at a point that is thefarthest from the heart among those discussed), for example, reflectionsof the pulse wave can be seen as 2531, and a smoothed delay at thebeginning of the systole can be seen as 2532. It is understood thatthese graphs represent decay and broadening of a pulse wavecorresponding to a healthy and unobstructed circulatory system, and areused to illustrate the operational advantages of placement of the sensoras close to the central aorta as feasible.

Bourgeois et al. proposed a method for calculating a beat-to-beat strokevolume (and hence the cardiac output) from measurements of andsubsequent calculations based on the beat-to-beat pulse pressure wave inthe central aorta. The chart of FIG. 26 (corresponding to FIG. 1 ofBourgeouis et al.) illustrates the “method for determination ofindividual beat-to-beat stroke volume from aortic pressure pulses”. Thesolid curve represents a normal pressure pulse from the central thoracicaorta. The dashed curve above the systolic portion of the pressure pulserepresents the hypothetical aortic pressure contour, resulting in apressure increment, ΔP_(ES), which would occur if all peripheral vesselswere to be closed during systole, thereby preventing systolic drainage.The pressure equivalent to the diastolic drainage in the steady statecase is (P_(ES)−P_(D)). The pressure equivalent to the total strokevolume. ΔP_(SV), is the sum of these values. K_(A) is a proportionalityfactor.” (quoted from the Bourgeouis et al.)

To calculate the stroke volume (SV), the following operations areperformed:

$\begin{matrix}{{SV} = {K_{A}\Delta\; P_{SV}}} & (50) \\{{SV} = {K_{A}\left\lbrack {\left( {P_{ES} - P_{D}} \right) + {\Delta\; P_{ES}}} \right\rbrack}} & (51) \\{{SV} = {K_{A}\left\lbrack {\left( {P_{ES} - P_{D}} \right) + \frac{SA}{\tau}} \right\rbrack}} & (52)\end{matrix}$

As discussed earlier, the time constant r can be calculated aftermeasurements are obtained. The other constant, K_(A), requires acalibration phase followed by the desired measurements if absolutereadings of cardiac stroke volume are required, or can remain as aconstant if only relative changes in stroke volume over time aredesired. The pressure values are then obtained from the measurement.

Regarding K_(A), Bourgeois 1976 used only a single calibration “based onthe first indicator-dilution measurement of cardiac output performed atthe beginning of an experiment” for a set of experiments. This “was usedfor all subsequent determinations of stroke volume.” The results of themeasurements demonstrated “the relative invariability of the K_(A)value” over the durations and under the conditions of the experiment.The invariability “is implied by the numerical similarity of thestandard deviations calculated from the entire set of data measuredduring short experiments as well as from those of extended duration.”Therefore, a single K_(A) derivation prior to measurement is adequate toproduce sufficient accuracy for many applications requiring absolutereadings of cardiac stroke volume.

Additionally, as the arterial pulse is propagating, at higher pressurethe arterial cross section increases slightly (depending on theelasticity of the aortic wall and the strength of the ventricularcontraction). The increase of the arterial cross section facilitates, inturn, the transmission of the EMW that traverses the sample, through thecharacteristic oxyhemoglobin. Because the impulse response is measuredwith such fine grain (see FIG. 8E), an attendant increase in attenuationoverall will be measured that corresponds to the vessel expansion. As aresult, there will be an increase in the width of that portion of theimpulse response that is caused by the material component thatcontributes a larger percentage to the optical path (between the opticalsource and the detector) during systole than that during diastole. Inthe present example, such material components if blood. Accordingly,both the peak concentration of the impulse response and the area underthe curve of the impulse response will increase. The increase of theaorta in width is, of course, dependent on the elasticity of thearterial walls, such elasticity is estimated from the rate of risebetween the points 2501 and 2502 of the graph 2500 of FIG. 25.Accentuating the contribution from the portions of the step functionthat vary with changes in effective optical path length throughout thecardiac cycle improves the measurement as a surrogate for pressure.

Accordingly, the oximetry system of the invention enables theidentification of the optimal site or location for measurement ofinstantaneous changes in oxygen saturation and creation of a surrogateaortic pressure pulse waveform, based on which the calculation of thecardiac stroke volume and cardiac output is carried out. Similarly,other pulse waveforms may be measured throughout the body.

Measurements are not restricted to the aforementioned physiologicalvariables. Other physiological variables can be measured using similartechniques. With the increased sampling frequencies, the transit timeand therefore also the velocity of blood may be estimated in real-time.

The shape of individual returned signals (see FIG. 8D) can be used tomeasure additional information regarding the central aortic vessels.e.g., their elasticity. This analysis of the shape and change of shapeof blood vessels is illustrated in FIG. 8E that shows a blood vesselboth constricted, 850, and dilated, 851. The constricted vessel 850 isoutlined with a fine dashed line and the dilated vessel 851 is outlinedwith a course dashed line. The time-dependent amplitude of the signalreturned by the constricted vessel 850 is shown in the graph as line 852and for the dilated vessel 851 as line 853. The primary observation fromthis figure is that the size of the vessel influences both the width ofthe curve representing the detected signal and the area under the curve.From this graph it can be seen that the determination of the change ofsize a blood-vessel and its elasticity is enabled, according to anembodiment of the invention, by measurement of the shape of the waveformof the impulse response of the blood-vessel.

As was alluded to before, and according to one embodiment of theinvention, the fundamental instrumentation block includes a “transceiverpair” formed by a transmitter and a single detector. The impulseresponse is derived for such a transceiver pair. At minimum, when thetransmitter operates at a single wavelength, the derivation of thisimpulse response enables the estimation of the average,wavelength-specific path length light travels through the tissue. (Thiscorresponds to the centroid of the true impulse response, i.e.calibration.) In addition, the fundamental optical properties of thetissue can be extracted including the scattering, absorption, andanisotropy of the tissue. These data represent data collected at a pointin time at a single wavelength. Once additional data is collected, in asimilar fashion, as a function of time, the changes of the physiologicalor optical parameter of the tissue with time can be determined (e.g., achange of the pulse pressure wave, see for example FIG. 26). Oneadditional result is the ability to characterize changes in theunderlying tissue as a function of time and physiology (see for example,FIG. 8E). According to an embodiment of the invention, the parametercharacterizing the tissue is spatially localized when such parameter,determined based on the measured impulse response of the tissue, iscoarsely spatially mapped to the tissue with the use of predetermined ortabulated frequency dependent attenuation characteristics (see FIG. 8D).

Extending the measurement further by employing multiple wavelengths oflight, determining the spectrally-dependent impulse response of thetissue, and determination of the spectral dependence of the parametercharacterizing the tissue enables, for example, the estimation ofconcentration of as many species as there are wavelengths in the EMWpropagating between the transmitter and the receiver of the pair, at onepoint in time. Extending the determination of the spectral dependence ofthe tissue parameter by carrying out multiple measurements separated intime, enables spatial localization of the parameter. (see FIG. 8D)

The next logical extension of the method providing for spatiallocalization of the tissue parameter is the use of a single transmitterand multiple receivers, spaced from the transmitter in different ways.For each new transceiver pair all previous steps occur. The localizationof the tissue parameter determined based on multiple overlapping orspatially shifted impulse responses results in increased density ofsampling in FIG. 8D, and enables the formation of an image of theunderlying tissue that references specific anatomical structures ofinterest. Yet the next fundamental extension of the method is thederivation of spatial distribution of the tissue parameter(s) based onimpulse responses determined at multiple transceiver pairs at multiplewavelengths of operation and multiple separations, enabling acomprehensive spatially dense map of tissue parameters estimated throughtime and space.

Referring again to FIGS. 8D and 24 A, it is appreciated that differentportions of a given material characteristic curve (such as curverepresenting concentration of a material component or a physiologicalparameter derived, in accordance with the embodiments of the invention,based on the measured or determined impulse response curve) areassociated with spatially-different portions of the tissue through whichthe EMWs propagated from the transmitter to the receiver. Putdifferently, the spatial location of a region of the tissue is coded inboth the measured impulse response curve and the material parametercurve calculated based on such impulse response curve. By analyzingdifferent portions of the calculated material parameter curve, suchparameter can be spatially mapped across the tissue (and, therefore,localized to a given portion of the tissue). In this case, theconclusion about the value of the sought-after parameter characterizingthe tissue is made on a dynamically localized spatial basis. Forexample, by considering the front portions of the concentration curves2421, 2422, through 2425 of FIG. 24A, the concentrations of respectivelycorresponding material species are defined in the regions of tissue thathas delayed the EMW by the least amount, while the tail portions of thecurves 2421, 2422 through 2425 contain information of the correspondingconcentration values related to the portions of the tissue that havedelayed the EMW the most. The integrated over time impulse responsecurve and the corresponding integrated concentration curve would beassociated with a spatially-averaged-over-the-tissue concentration ofthe given material parameter.

In reference to FIG. 27, an embodiment of the invention directed to themeasurements of the sample that are spatially localized to the region ofthe sample that with the use of the impulse responses is initiated atstep 2700; the impulse response (IR) (t) is measured for each ofmultiple wavelengths λi and at each time point of interest forgenerating a time-derived sequence, at step 2710. An example of atime-derived signal would be deriving the contour of a pressure wavesurrogate from a sequence of samples [S_(T1), . . . . , S_(TN)], witheach S_(T) is conceptually composed of the integral ∫IR_(T) _(X) (t)dtof a particular measurement IR_(T) _(X) (t) of the impulse responsebased at time Tx. (The resulting time-derived sequence would nominallyresemble the form represented by the graph 2500 of FIG. 25).

The next step is to combine the impulse response (t) with the sampledλ_(i) values to determine impulse response values (t) representing aproperty of the material component of interest.

From so-determined impulse response (t) sampled values, a valuerepresenting the bulk property is formed, at step 2730, by summing thevalues of the individual samples 2720. Then, the derived time-basedsequence is formed at 2740. This sequence, corresponding to, forexample, a pressure wave surrogate, is based on the bulk, not localized,averaged over the tissue cross-section property value 2730. The resultmay be similar to, for example, the graph 2500 of FIG. 25. Using thatgraph representation as an illustration, the process would proceed toidentify time points of interest (MVT_(1-N)) at step 2750. In this case,if the details of the local maximum of the measured curve are ofinterest, the data values corresponding to the points on the curve 2500of FIG. 25 and labeled as 2501, 2502, 2503, 2504, and 2506 are likelythe points of interest that could be identified as inflection points ofthe time-derived property value. If, on the other hand, achievingprecision in time is the goal, then the local cohort of data surroundingthese points of interest can be chosen.

Once the specific points in time are chosen, analysis of the data can beconducted to determine the impulse response delay time values thatexhibit the highest variance, or stated alternately, those that providethe greatest contrast between the areas of the derived time-basedproperty of interest 2760. This can also be described as choosing, atstep 2760, the regimes of larger cross variance from the sampled impulseresponses at MVT_(1-N). From so obtained information, produce aweighting function W(t) over the impulse response time range at step2770. Using the weighting function together with the impulse response(t) λ_(i) values, a refined property value can be determined as∫IR(t)*W(t)dt, at step 2780. Finally, a refined derived time-basedsignal (such as the refined pressure wave) can be determined based onthe refined property value, at step 2790. This completes the process ofproducing the weighting function.

The weighting function, once determined, can be reused for subsequentmeasurements. Many events will cause a need to recalculate the weightingfunction, including, but not limited to motion of the sample.

Implementations of the invention can employ a processor controlled byinstructions stored in a tangible memory to perform the steps ofoperation of the system described above. The memory may be random accessmemory (RAM), read-only memory (ROM), flash memory or any other memory,or combination thereof, suitable for storing control software or otherinstructions and data. In an alternative embodiment, the disclosedsystem and method may be implemented as a computer program product foruse with a computer system. Such implementation includes a series ofcomputer instructions fixed either on a tangible non-transitory medium,such as a computer readable medium (for example, a diskette, CD-ROM,ROM, or fixed disk) or transmittable to a computer system, via aninterface device (such as a communications adapter connected to anetwork over a medium). Some of the functions performed during theexecution of the method of the invention have been described withreference to flowcharts and/or block diagrams. Those skilled in the artshould readily appreciate that functions, operations, decisions, etc. ofall or a portion of each block, or a combination of blocks, of theflowcharts or block diagrams may be implemented as computer programinstructions, software, hardware, firmware or combinations thereof. Inaddition, while the invention may be embodied in software such asprogram code, the functions necessary to implement the invention mayoptionally or alternatively be embodied in part or in whole usingfirmware and/or hardware components, such as combinatorial logic.Application Specific Integrated Circuits (ASICs), Field-ProgrammableGate Arrays (FPGAs) or other hardware or some combination of hardware,software and/or firmware components.

While the invention has been described through the above-presentedexamples of embodiments, it will be understood by those of ordinaryskill in the art that modifications to, and variations of, theillustrated embodiments may be made without departing from the inventiveconcepts disclosed herein. Furthermore, disclosed aspects, or portionsof these aspects, may be combined in ways not listed above. Accordingly,the invention should not be viewed as being limited to the disclosedembodiment(s).

What is claimed is:
 1. An apparatus for collection of optical data, theapparatus comprising: a source of light including at least onetransmitter configured to emit a beam of first light at at least onecorresponding wavelength along a first axis; a device including a firstset of selectively and individually re-orientable reflecting elementsdisposed parallel to a first plane when each of said reflecting elementsfrom the first set is in a first reference orientation; a reflectingelement from the first set configured to redirect the first light,received from said beam, towards a predetermined point, wherein thepredetermined point is separate from the first set; an optical receiverdisposed in optical communication with said device to receive secondlight, the second light being the first light that has been scattered atthe predetermined point; and a computer-readable processor and tangible,non-transitory storage medium containing program code thereon which,when downloaded onto the processor, causes the processor to receivereference data representing an empirically-defined spectrally-dependentcharacteristic of an object present at said predetermined point; todetermine an impulse response of the object, the impulse response beingassociated with emission of the first light by the source; and todetermine, as a function of time, at least one of (i) awavelength-dependent material parameter characterizing the object and(ii) a position-dependent parameter characterizing the object, based onthe determined impulse response and the reference data.
 2. An apparatusaccording to claim 1, further comprising a second set of selectively andindividually re-orientable reflecting elements disposed parallel to asecond plane when each of said reflecting elements from the second setis in a second reference orientation, the second plane being transverseto the first plane and substantially parallel to the first axis.
 3. Anapparatus according to claim 2, wherein a reflecting element from thesecond set is configured to receive said second light and to redirectthe second light towards the optical receiver.
 4. An apparatus accordingto claim 2, configured to vary an optical path associated withpropagation of the first light from the reflecting element from thefirst set to the predetermined point and then to the reflecting elementfrom the second set as a result of spatially reorienting at least one of(i) the reflecting element from the first set and (ii) the reflectingelement from the second set.
 5. An apparatus according to claim 1,wherein the device comprises a baffle disposed to prevent said opticalreceiver from i) receiving said first light that has not been scatteredat the predetermined point and from ii) receiving said second light thathas interacted with the reflecting element from the first set.
 6. Anapparatus according to claim 1, configured to satisfy at least one ofthe following conditions: (i) different reflecting elements of the firstset are configured to be reorientable synchronously with one another;and (ii) different reflecting elements of the second set are configuredto be reorientable synchronously with one another.
 7. An apparatusaccording to claim 1, wherein said predetermined point is a point at asurface of an optical disk.
 8. An apparatus according to claim 1,configured to vary an optical path associated with propagation of thefirst light from the reflecting element from the first set through thepredetermined point to the optical detector as a result of spatiallyreorienting the reflecting element from the first set.
 9. An apparatusaccording to claim 2, wherein at least one of the first and second setsincludes only one, single reflecting element.
 10. An apparatus accordingto claim 2, wherein the first and second sets includes reflectingelements disposed in respectively-corresponding first and second rows.11. An apparatus according to claim 2, wherein first and second sets ofreflecting elements are portions of a two-dimensional array ofreflecting elements.
 12. An apparatus according to claim 2, wherein thedevice is configured to reposition all reflecting elements from both thefirst and second sets to assume the second reference orientation.
 13. Anapparatus according to claim 2, wherein the device is configured tospatially-continuously vary an angle formed by the first and secondplanes in response to changing at least one of the first and secondorientations.
 14. An apparatus according to claim 2, wherein the deviceis configured to change orientation, of at least one of the reflectingelements of the first set and the reflecting elements of the second set,in an angularly-discrete fashion.
 15. An apparatus for collection ofoptical data, the apparatus comprising: an optical detector; a source oflight including at least one transmitter configured to emit a beam offirst light at at least one corresponding wavelength along a first axis;an optical system containing an array of spatially-reorientablereflecting elements and configured to deliver the first light, emanatingfrom the source of light, to the optical detector while scanning anobject with the beam by reorienting first and second groups of saidreflecting elements between first and second reference orientations,respectively; and a data-processing system in operable communicationwith the optical detector and configured to acquire an output producedby the optical detector in response to being irradiated with said firstlight that has been scattered by the object, the data-processing systemincluding a computer-readable processor and tangible, non-transitorystorage medium containing program code thereon, wherein said programcode, when loaded onto the processor, causes the processor to receivereference data representing an empirically-defined spectrally-dependentcharacteristic of the object; to determine an impulse response of theobject, the impulse response being associated with emission of the beamof light by the source; and to determine, as a function of time, atleast one of (i) a wavelength-dependent material parametercharacterizing the object and (ii) a position-dependent parametercharacterizing the object, based on a determined impulse response andthe reference data when the object is positioned at an intermediatepoint along a direction of propagation of the beam between the first andsecond groups of the reflecting elements, wherein the first group ofreflecting elements is positioned to receive the beam from the source oflight, wherein the second group of reflecting elements is configured toredirect the first light, that has been scattered by the object, towardsthe optical detector.
 16. An apparatus according to claim 15, whereinthe device is configured to spatially-continuously reorient reflectingelements from at least one of the first and second groups between thefirst and second reference orientations, a surface of the reflectingelement in the first orientation being transverse to the first axis, asurface of the reflecting element in the second orientation beingparallel to the first axis.
 17. An apparatus according to claim 15,further comprising a baffle disposed in a first plane to block the firstlight from propagating directly between the source and the opticaldetector, wherein the first plane is perpendicular to a plane definedwhen all reflecting elements of the optical system are in the secondreference orientation.
 18. An apparatus according to claim 16,configured for use with the object that includes a biological tissue,and wherein data-processing system is configured to determine an opticalproperty of the object without the use of an arterial pulse in theobject, the optical property including at least one of a)wavelength-dependent scattering, b) absorption, and c) anisotropy of theobject, said optical property having been determined without the use ofan arterial pulse in the sample.
 19. An apparatus according to claim 17,wherein the data-processing system is configured to determine the atleast one wavelength-dependent material parameter characterizing theobject based on an averaged wavelength-dependent path length that hasbeen estimated with the use of the determined impulse response andscattering of light from said beam in the object, said averagedwavelength dependent path being different from a straight path between areflecting element from the first group and a reflecting element of thesecond group.
 20. An apparatus according to claim 19, wherein thedata-processing system is further configured to sample a curverepresenting the determined wavelength-dependent material parameter as afunction of time to obtain values of said material parameter at samplingpoints in time, and to form a weighting function by calculating variancevalues among the values of said material parameter at the samplingpoints in time, the weighting function being defined by said variancevalues as a function of time.
 21. An apparatus according to claim 20,wherein the data-processing system is further configured to modify thedetermined wavelength-dependent material parameter characterizing theobject based on (i) the computed impulse response that has been weightedby the weighting function and (ii) the reference data.
 22. A method forcollecting optical data, the method comprising: redirecting a beam oflight, that has been received from an optical source, with a firstspatially-reorientable reflecting element from an array of reflectingelements of an optical device, towards a first point of an object,wherein said beam of light is characterized by at least one wavelength;receiving first light from said beam, that has been scattered by theobject, from a second point of the object with a secondspatially-reorientable reflecting element from said array to redirectsaid first light towards an optical detector; and with a data-processingsystem including a computer-readable processor: acquiring an outputgenerated by the optical detector in response to being irradiated withsaid first light, the output representing an optical property of theobject determined along a path from said first reflecting element tosaid second reflecting element through the object, said path defined bya gradient of concentration of a material component in the object;receiving, with the data-processing unit of the system, reference datarepresenting an empirically-defined spectrally-dependent characteristicof the object; determining an impulse response of the object, theimpulse response being associated with emission of light from theoptical source; and determining, as a function of time, at least one of(i) a wavelength-dependent material parameter characterizing the objectand (ii) a position-dependent material parameter characterizing theobject, based on a determined impulse response and the reference data.23. A method according to claim 22, wherein the object includesbiological tissue, and wherein the acquiring the output includesacquiring the output representing the optical property including atleast one of a) wavelength-dependent scattering, b) absorption, and c)anisotropy of the object, said optical property having been determinedwithout the use of an arterial pulse in the object.
 24. A methodaccording to claim 22, wherein the determining, as a function of time,the wavelength-dependent material parameter characterizing the objectincludes determining said wavelength-dependent material parameter basedon an averaged wavelength-dependent path length that has been estimatedwith the use of the determined impulse response and scattering in theobject, wherein said averaged wavelength dependent path is differentfrom a straight path between the first and second points of the object.25. A method according to claim 24, further comprising: with thedata-processing system, sampling a curve representing the determinedwavelength-dependent material parameter as a function of time to obtainvalues of said parameter at sampling points in time, and forming aweighting function by calculating variance values among the values ofsaid determined wavelength-dependent material parameter at samplingpoints in time, the weighting function being defined by said variancevalues as a function of time.
 26. A method according to claim 25,further comprising: modifying the determined wavelength-dependentmaterial parameter characterizing the object based on (i) a computedimpulse response that has been weighted by the weighting function and(ii) the reference data.
 27. A method according to claim 22, wherein theobject includes biological tissue; and wherein said determining thewavelength-dependent material parameter includes determining one or moreof a level of oxygenation of the object and a concentration of a chosenspecies in the object, and further comprising: at different wavelengths,modulating said beam of light with first and second excitation sequencesapplied to the source of light, the first and second excitationsequences being orthogonal to one another.
 28. A method according toclaim 22, wherein the determining the impulse response includesdetermining the impulse response of the object based on the output, theoutput being a time-dependent characteristic of a pair formed by thesource of light and the optical receiver, the pair identified at thedata-processing unit of the system by a respectively correspondingexcitation sequence series of pulses produced by the source of light, afirst excitation sequence series of pulses at a first wavelength beingorthogonal to a second excitation sequence series of pulses at a secondwavelength.
 29. A method according to claim 22, wherein the objectincludes a biological tissue; and further comprising adding a chemicalsubstance to the object; and wherein said determining thewavelength-dependent material parameter includes determining a change inthe wavelength-dependent material parameter caused by said adding.
 30. Amethod according to claim 29, further comprising generating said beam oflight as electromagnetic pulses in a defined excitation sequence andtransmitting said pulses through a delay line with multiple delay tapsto enhance said determining the change in the at least onewavelength-dependent material parameter and to reduce a cost of thedetermining of said impulse response.
 31. A method according to claim22, wherein the object includes a biological sample, and furthercomprising repositioning at least one of the source of light and theoptical detector along a surface of the biological sample to enhancelocalization of said determining the wavelength-dependent materialparameter based on an image of the biological sample created as a resultof said repositioning, the image representing a structure of thebiological sample.
 32. A method according to claim 22, wherein theobject is homogeneous.
 33. A method according to claim 22, furthercomprising blocking light from said beam of light from propagatingdirectly from the source of light to the optical detector along asurface of said array.